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现代 社会 ， 数 学 理论 已 经 发 展 得 十 分 完善 ， 各 种 数学 工 
具 对 现实 的 生产 显示 出 很 强 的 指导 作用 。 然 而 ， 在 数学 的 理 
论 工具 与 实际 生产 的 结合 过 程 中 ， 效 率 低下 的 传统 手工 计算 
方法 成 为 最 大 的 障碍 。 随 着 计算 机 技术 的 进步 ， 各 种 数学 软 
件 也 如 雨后春笋 般 展 示 出 勃勃 生机 。 

MATLAB 是 一 款 简 单 易 用 的 数学 软件 ,由 美国 mathworks 
公司 发 布 ,主要 面 对 科 学 计算 、 可 视 化 以 及 交互 式 程序 设计 的 
高 科技 计算 环境 ,具有 强大 的 数据 处 理 能 力 和 出 色 的 图 形 处 
理 功能 。 学 习 并 掌握 它 , 有 助 于 人 们 从 繁重 的 数学 计算 中 解 
脱出 来 ,把 更 多 的 精力 投入 数学 理论 的 学 习 和 研究 。 

目前 ，MATLAB 在 我 国 的 应 用 越 来 越 广泛 ， 很 多 科学 领 
域 都 用 它 来 完成 数值 、 图 像 和 信号 的 分 析 与 处 理 ， 而 且 其 独 
特 的 、 功 能 丰富 的 应 用 工具 箱 为 用 户 提供 了 大 量 方便 实用 的 
处 理工 具 ， 受 到 广大 用 户 的 一 致 好 评 。 

本 书 主要 介绍 MATLAB 基础 与 应 用 ,目的 就 是 要 培养 读 
者 的 数值 计算 和 数据 处 理 能 力 ,加 强 读者 用 数学 工具 分 析 和 
解决 问题 的 意识 。 本 书 在 编写 过 程 中 ,力求 体现 以 下 特点 : 

1) 语言 简练 ， 趣 味 性 强 ， 易 于 阅读 与 接受 。 

2) 知识 难度 不 大 ， 适 合 初 学 者 使 用 。 

由 于 水 平 有 限 ， 时 间 人 仓促， 书 中 难免 有 不 足 之 处 ， 欢 迎 
读者 批评 指正 。 
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第 1 章 MATLAB 简介 





和 其 他 数学 软件 相 比 ，MATLAB 具有 简洁 直观 、 使 用 方 
便 、 符 合 人 们 的 习惯 思维 、 库 函数 丰富 等 优点 。 除 卓越 的 数 
值 计 算 功 能 外 ，MATLAB 还 具有 专业 水 平 的 符号 计算 、 文 字 
处 理 、 可 视 化 建 模仿 真 等 功能 ， 儿 乎 能 解决 所 有 的 工程 计算 
问题 。 在 国外 ，MATLAB 软件 已 经 受 了 许多 年 的 考验 。 在 欧 
美 等 高 校 ，MATLAB 已 经 成 为 线性 代数 、 自 动 控 制 理论 、 数 
理 统计 、 数 字 信号 处 理 、 时 间 序 列 分 析 、 动 态 系 统 仿 真 等 高 
级 课程 的 基本 教学 工具 。 

MATLAB 是 Matrix Laboratory 的 简写 ， 意 为 矩阵 实验 室 。 
它 产生 于 20 世纪 70 年 代 后 期 ， 是 美国 New Mexico 大 学 计 
算 机 系 主任 Cleve Moler 在 给 学 生 讲授 线性 代数 课程 时 ， 为 
学 生 使 用 EISPACK 和 LINPACK 而 编写 的 接口 程序 。1984 
年 ， 由 Маф Works 公司 正式 推 向 市 场 ， 并 不 断 更 新 完善 。 


1.1 MATLAB 的 启动 






























































双击 桌面 上 的 MATLAB 图 标 ， 显 示 如 图 1-1 所 示 MAT- 
LAB 工作 界面 ， 其 中 右边 窗口 为 命令 窗口 (Command Win- 
dow) ， 用 于 输入 操作 命令 ; 在 左下 方 窗口 为 历史 记录 窗口 
(Command History) ， 保 留 自 安装 时 起 所 有 命令 的 历史 记录 ， 
并 标明 使 用 时 间 ， 以 方便 使 用 者 查询 。 双 击 某 一 行 命 令 ， 即 
在 命令 窗口 中 执行 该 命令 ; 左上 方 窗口 为 工作 空间 管理 窗口 
( Workspace) ， 显 示 所 有 目前 保存 在 内 存 中 的 MATLAB 变量 




















.2. MATLAB 基础 与 应 用 





的 变量 名 、 数 据 结构 、 字 节 数 以 及 类 型 ， 不 同 的 变量 类 型 分 
别 对 应 不 同 的 变量 名 图 标 。 


工作 空间 
管理 窗口 















1.2 常用 操作 


1. 恢复 默认 的 工作 界面 
有 时 可 能 由 于 误 操 作 ， 更 改 了 MATLAB 的 工作 界面 ， 
要 恢复 到 默认 的 工作 界面 ， 操 作 如 图 1-2 所 示 。 





File Edit Debug Desktop Window Help 
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Move Command Window 
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Resize Command Window 
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图 1-2 
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如 果 只 是 工作 区 窗口 和 历史 区 窗口 丢失 ， 只 要 打开 界面 
菜单 栏 的 desktop, 分 别 单 击 Command Window, Command 
History, Workspace 即 可 。 

2. 清除 操作 





>> clc p 清除 命令 窗口 中 的 所 有 内 容 
>> clear % 清除 内 存 中 的 所 有 MATLAB 记忆 
>> elf 和 清除 图 形 窗 口内 容 


其 中 ,符号 “% ”表示 注释 ， 在 命令 窗口 中 不 运行 。 

3. 调用 以 前 运行 过 的 语句 

(D 按 键盘 上 的 上 下 方向 键 ; 

G@ 在 左下 角 的 历史 记录 窗口 双击 欲 调 用 语句 ; 

(3 在 右边 的 命令 窗口 或 左下 角 的 历史 记录 窗口 选中 语 
句 ， 顺 次 进行 如 下 操作 : 


而 在 命令 窗口 有 击 paste 


4. 其 他 
QO 请 求 帮助 文件 >> help 请求 内 容 ; 

如 : >>help efun J % 关 于 基本 水 数 的 帮助 信息 
>>help exp «| % 指数 函数 exp 的 详细 信息 

@) 超 文本 格式 的 帮助 文件 : >> doc 请求 内 容 ; 
1: >>doc еп] 乡 以 超 文 本 格式 显示 关于 基本 
函数 的 帮助 信息 
>> дос exp .J % 以 超 文 本 格式 显示 关于 指数 
函数 的 帮助 信息 

(ЗК В: >> lookfor 请 求 内 容 ; 
当 要 查找 具有 某 种 功能 但 又 不 知道 准确 名 字 的 指令 时 ， 
help 的 能 力 就 不 够 了 ， 此 时 可 使 用 ”lookfor ”命令 。lookfor 可 
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以 根据 用 户 提供 的 完整 或 不 完整 的 关键 词 ， 去 搜索 出 一 组 与 
之 相关 的 指令 。 
如 : >>lookfor integral -| ЖЕНА 
>> lookfor fourier J % 查找 能 进行 侵 里 叶 变 换 
的 指令 
DRRR: >> demo 请 求 内 容 。 
命令 help, doc, lookfor 都 有 各 自 的 特点 ， 其 中 help 与 
doc 只 是 显示 方式 不 同 。help、doc 的 请 求 内 容 必 须 完 整 准 
确 ， 而 lookfor 后 面 的 请 求 内 容 可 以 不 完整 。 


1.3 常量 与 变量 


МАТГАВ 语言 本 身 具有 一 些 预 定义 的 变量 值 ， 这 些 特 殊 
的 有 特定 值 变量 称 为 常量 。 表 1-1 给 出 了 MATLAB 语言 中 经 
常 使 用 的 一 些 常 量 。 



























































表 1-1 
常 量 表示 数值 

pi Ја т 

inf 正 无 穷 大 

NaN 表示 不 定 值 

i,j 虚数 单位 

eps 计算 机 的 最 小 数 

realmax 最 大 可 用 正 实数 

realmin 最 小 可 用 正 实数 



































MATLAB 语言 中 的 变量 是 由 字母 、 数 字 、 下 划 线 组 成 ， 
主要 命名 规则 为 : 
中 以 字母 开头 ; 
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@ 区 分 大 小 写 。 
MATLAB 语句 有 两 种 最 常见 形式 : 
Ф» 变量 = 表达 式 ; 
运行 结果 显示 为 “变量 ORD 
— >> 表达 式 ; 
行 条 结果 显示 为 “ans = =...” 


























站， “ans” 是 指 当 前 的 计算 结果 ， 奉 计算 时 用 户 没 有 对 表 












































达 式 设 定 变 量 ， 系 统 就 自动 赋 当 前 结果 给 “ans” 变量。 如: 
>> a=1 +2 >>] +2 
a= ans = 
3 3 
14 算术 运算 符 
MATLAB 的 加 、 减 、 乘 法 运算 符 的 输入 和 通常 的 电脑 输 
人 是 一 致 的 。 除 法 运算 分 左 除 (\ MERC), 2/3 是 2 除 








以 3， 而 2\ 3 实际 是 3 除 以 2。 为 了 避免 混 消 ， 对 一 般 除 法 





























运算 采取 前 者 。 乘 方 运算 符 为 电脑 键盘 上 的 “*”。 算 术 运 算 
符 输入 方式 见 表 1-2. 
表 12 
MATLAB 输入 
运算 符 
JE 阵 数 组 
加 十 十 
减 z = 
乘 
除 / ./ 
乘 方 
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MATLAB 的 运算 分 矩阵 运算 和 数组 运算 两 种 。 线 性 代数 
HE m ÍF n Яо ЕНЕЛ АЯН РЕКЕ о ЦП: 


а v Qin 





只 有 1 行 或 1 列 的 矩阵 叫做 向 量 或 数组 。MATLAB 的 基 
本 数据 单位 是 和 矩阵 ， 因 此， 正常 的 运算 是 矩阵 运算 ， 在 运算 
符 前 加 点 的 运算 是 数组 运算 。 

在 MATLAB 实际 操作 时 如 果 刻 意 区 分 运算 符 前 是 否 加 

可 能 会 使 问题 复杂 化 ， 制 造 出 不 必要 的 混乱 。 符 合 人 们 的 
и 大 优点 ， 因 此 ， 在 操作 时 ， 可 先 
按 不 加 点 的 方式 进行 输入 ， 如 果 输 入 没有 其 他 错 ， 而 命令 运 
行 不 了 ,那么 其 运算 就 可 能 是 数组 运算 ,在 相应 的 运算 符 前 

































































ЛИ“. "АА 
如 : >> х= -5: 0. 5:5; 
>> y =x2 
2??? Error using = = > mpower Inputs must be а scalar and 


a square matrix. То compute elementwise POWER, use POWER 
(. `) instead % 这 行 英文 红色 语句 是 错误 警告 ,警告 命令 
输入 有 误 , 并 提示 用 “. 7 РНЕ? 


>>у=х. ^2 
1.5 逻辑 运算 符号 


逻辑 运算 是 MATLAB 中 数组 运算 的 一 种 运算 形式 ， 也 
是 几乎 所 有 的 高 级 语言 普遍 使 用 的 一 种 运算 。 它 的 符号 运算 
符 、 功 能 及 函数 名 见 表 1-3。 
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表 13 
符号 运算 符 Ж: TE 
= 等 于 ed 
E 不 等 于 ne 
< 小 于 lt 
> RF gt 
<= 小 于 等 于 le 
>= 大 于 等 于 ge 
& 逻辑 与 and 
| 逻辑 或 or 
~ 逻辑 非 not 
说 明 : 在 算术 运算 、 比 较 运 算 和 逻辑 与 、 或 、 非 运算 


中 ， 它 们 的 优先 级 关系 先后 为 :比较 运算 、 算 术 运 算 、 逻 辑 












































与 或 非 运 算 。 
1.6 其 他 常用 符号 
其 他 常用 符号 见 表 1-4。 
表 1-4 
符号 “| MATLAB 输入 用 ж 
55 ; 分 隔 变 量 、 表 达 式 、 和 矩阵 的 列 
分 号 ; 分 隔 命令 行 而 不 显示 运行 结果 ， 分 隔 矩 阵 的 行 
单 引号 定义 字符 串 
x =a: b: c 
а 表示 x 从 a IEK b 取 值 至 
等 号 = 变量 赋值 
百 分 号 % 命令 注释 
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( 续 ) 
符号 “| MATLAB 输入 я я 
3 个 句点 Зе 续 行 
到 括号 () 区 分 运算 次 序 
方 括号 г) 构成 矩阵 或 向 量 
1 2 3 
例 1.1: 在 MATLAB 中 输入 矩阵 4=|4 5 6. 
7 8 9 


MATLAB 程序 如 下 : 
>> А=[1,2,3;4,5,6;7,8,9] % 行 与 行 之 间 用 分 号 
分 隔 , 每 行 的 ( 列 ) 元 素 间 用 逗号 分 隔 。 














А = 
3 
6 
7 9 
实 训 题 
1996 
1 计算 表达 式 ”的 结果 。 
2sin0. 3 ров 
2. 输入 表达 式 一 一 一 ， 并 运行 。 





25 
1 4 
3: кни 5 
3 6 
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RAAY MATLAB 计算 与 作 图 


2.1 基本 初等 函数 的 输入 


在 MATLAB 中 ， 函 数 输入 的 总 体 原则 是 将 变量 整体 用 
括号 括 起 来 。 如 cos2x; 的 MATLAB 输入 为 cos (2 x (х^3)), 
x 本 来 不 需要 用 插 号 括 起 来 ,但 括 起 来 后 ， 运 算 次 序 更 加 




















清晰 。 常 见 初等 函数 的 具体 输入 方式 见 表 2-1。 










































































表 2-1 
名 称 AT MATLAB 命令 # Ë 
х“ ха x 括 不 括 都 可 以 
i Зан(х) (102) | 1⁄2 必须 括 起 来 
а" а^х x 括 不 括 都 可 以 
指数 函数 
e” exp(x) 不 能 用 ex 
lnx log(x) 对 数 函 数 只 有 e、2、 
10 三 个 底 ,其 他 底 的 情 
对 数 函 数 log, x log2 (x) 况 需 用 换 底 公式 : 
log.b 
logiox log10(x) РЕР а 
sinx ,cosx sin(x) ,cos( x) pe es 
= ff K% tanx ,cotx tan(x) ,cot( x) B м. 的 是 
secx ,CSCX зес(х) ,csc(x) 
arcsinx , arccosx asin( x) ,acos( x ) 在 三 角 函数 输入 前 
М = Ж | arctanx ,arccotx atan(x) ‚асоЕ(х) 加 E = i Z< RU 
a 




















arcsecx ,arCCSCX 





азес(х) ,acsc( x) 
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2.2 系统 运算 与 操作 函数 的 输入 


在 MATLAB 中 ,通常 以 由 基本 初等 函数 扩展 的 数学 函 
数 作为 处 理 的 对 象 。 此 外 ，MATLAB 系统 还 设计 了 具有 运算 
和 操作 性 质 方面 的 函数 ， 它 们 常 作为 处 理 的 工具 。 这 类 函数 
常见 的 有 以 下 几 种 : 

QO 绝对 值 函 数 |x |: abs(x); 

(ОР РА sign x: sign(x) ; 

@ 求 和 函数 : sum; 

СОЖ: prod; 

@ 求 最 大 值 : max; 

@ 求 最 小 值 : min, 


2.3 国 数 值 的 计算 
(1) 数值 计算 方式 











>>x=…; ФМЛ x 的 数值 (不 能 为 字母) 
>>y =… PIRA y 的 表达 式 (表达 式 中 除 x 外 不 能 有 其 
他 字母 ) 


е. 
0] 2.1: Úr y = 3⁄2 о 用 МАТГАВ 计算 


y(1) 的 值 。 
MATLAB 程序 如 下 : 
>> x=1; 
>> у=3 x (x2) -2/(3^х) + (2°x)/3 -4 * exp(2 * x) 


– 26. 5562 
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xX +l, х<0 
0] 2.2: у= (2° -1, 0 <х<10 ,用 MATLAB 计算 
2х +3, х>10 
y(5) КИК. 
MATLAB 程序 如 下 : 
> х=5; 
> if х<0 
y=x2 +1 
elseif x>0& x<=10 % else 与 计 之 间 不 能 空格 ,否则 要 
用 两 个 end 
у=2^х – 1 
else 
у=2 *жх+3 


end 


31 
(2) 符号 计算 方式 
>> syms x 其 他 字母 % ë x 和 其 他 字母 为 符号 


>>y=Í(x); % 输 入 y 的 表达 式 
>> subs (y,x,a) % 计 算 x=a 时 y 的 值 


注意 : 用 符号 计算 时 ， 对 数 函 数 只 识别 以 e 为 底 的 对 数 。 如 
果 要 计算 在 多 个 点 *=al，…,*=& Жу 的 值 ， 则 可 用 
La, а, 0, а, а, ША ап = ”是 以 符号 形式 
给 出 时 ， 输 入 double(ans) 即 可 得 到 数值 型 结果 。 


3 1 
例 2.3: 设 y=3lnx” – logyx + logio | 十 ae (4x), 
x 


Я MATLAB 计算 y(1)、y(2) 的 值 。 
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MATLAB 程序 如 下 : 
>> syms x 
>> у= 3 ж log(x*2) - (log (x)/log (2)) * (log (3/x)/ 
log(10)) +log(4 *x)/(2 *log(3)); 
>> subs (y,x,[1,2] ) 
ans = 
0. 6309 4. 9292 
例 2.4: Ш у=л” -2ax, Л] MATLAB H у (а), y(b) 
的 值 。 
MATLAB 程序 如 下 : 
>> syms xab 
>> у=х2 -2 *a*x; 
>> subs (y,x,[a,b] ) 
ans = 
[ -a2, b2 - 2*a*b] 
这 个 例子 告诉 我 们 ， 用 符号 计算 方法 时 ， 相 关 式 子 中 的 所 有 
字母 都 要 先 定义 为 符号 。 


2.4 ЖИ 


MATLAB 有 很 强 的 图 形 功 能 ， 可 以 方便 地 实现 数据 的 视 
觉 化 。 下 面 着 重 介 绍 二 维 图 形 的 画 法 。 

1. 一 般 函 数 y =f(x) 的 作 图 (二 维 ) 

(1) 作 图 基本 形式 ”二 维 图 形 的 绘制 是 MATLAB 语言 
图 形 处 理 的 基础 ，MATLAB 最 常用 的 画 二 维 图 形 的 命令 是 
plot，MATLAB 命令 格式 : 
>> x=a:c:b Ф х А а, b], KX c 
>> y=f(x); % W у 的 表达 式 , 表达 式 中 的 运算 
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符 加 点 


>> plot(x，y) 
如 : 


的 光滑 的 正弦 曲线 。 


% l| Hi pK 2 BJ sl 


% 生成 一 组 线性 
等 距 的 数值 


>> x =linspace(0, 2 * pi, 30); 


>> y =sin(x); 
>> plot(x, y) 
生成 的 图 形 如 图 2-1 Рук, ЕГО, 2т] Е 30 个 点 连 成 





7! 














(2) EZER ”在 同一 个 画面 上 可 以 画 许多 条 曲线 ， 只 


需 多 给 出 几 个 数组 ，MATLAB 命令 程序 格式 : 


>> x=a:c:b 

у 

>> ріоі(х,у1,х,у2) % 在 同一 平面 画 出 两 个 函数 的 
图 像 

如 : 
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>> x =0:pi/l5:2 * pi; 

>> yl =sin(x) 

>> y2 =cos(x); 

>> plot(x,yl ,x,y2) 
则 可 以 画 出 图 2-2。 


и. 








图 22 


注 : 中 也 可 用 hold оп 语句 ， 达 到 作 多 重 图 的 效果 
MATLAB 命令 格式 : 

>> plot (表达 式 1,[a,b]) 

>> hold on 

>> plot (表达 式 2,[a,b]) 

@ 如 果 要 在 一 个 画布 上 作 开 个 小 图 , 则 可 用 subplot(m， 
n,k), MATLAB 命令 格式 : 

>> subplot(m,n,k); 

fplot (' 表 达 式 1,[a,b]") 

即 表 达 式 1 所 表示 的 曲线 画 在 mm 行列 第 个 位 置 上 
(МЕЖ, 再 从 上 至 下 计数 ) 。 





` 
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(3) 作 图 的 线 型 和 颜色 为 了 适应 各 种 绘图 需要 ， 
MATLAB 提供 了 用 于 控制 线 色 、 数 据点 和 线 型 的 3 组 基本 参 
数 。 它 的 使 用 格式 如 下 : plot(x, у, 'color_ point_ linestyle'), 
具体 参数 见 表 2-2。 











































































































表 22 

b 蓝 色 m 紫红 色 

c 青色 f 红色 

5 绿色 w 白色 

k 黑色 y 黄色 

= 实 线 ( 默认) : 点 连 线 

= 点 画 线 == 虚线 

点 s 正方 形 

+ 十 字号 d 菱形 

o( 字 母 ) 5 h 六 角形 

* 星 号 p 五 角 星 

x( 字 母 ) 又 号 > 右 三 角 
(4) 作 图 的 网 格 和 标记 ”在 一 个 图 形 上 可 以 加 网 格 、 标 


题 、x 轴 标记 、y 轴 标 记 ， 用 下 列 命令 完成 这 些 工 作 。 
>> х =linspace(0, 2 *pi, 30); y=sin(x); z=cos(x); 
>> plot(x, y, x, z) 





>> grid % 网 格 
>> хјаһе (АХ) % 横 坐标 标签 
>> ylabel( АҮ) % 纵 坐标 标签 
>> title('Sine 和 Cosine 图 像 ') % 标 题 


注 : 如 果 要 使 图 形变 得 更 加 美观 ,也 可 作 一 些 技巧 性 的 处 理 。 
如 果 限 制 画 布 , 则 需 在 输入 plot 语句 前 输入 >> axis([a,b,e， 
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d] ) ,这 个 命令 是 将 图 形 限 制 在 La,b] x[c,d] E, HF a,b, 
c,d 必须 是 数值 。 
由 此 生成 图 2-3 : 











Sine 和 Cosine 图 像 





E ! ! | Е: 


纵 坐 标 Y 和 z 

















图 23 

2. 特 殊 函 数 的 作 图 

人 作 参 数 方程 人 (的 图 形 ， 也 可 以 用 pot 命令 
其 MATLAB 命令 格式 ; 
>>t=a:c:b 


>>x=f(t); y=g(t); 
>> plot (x,y, 'S') % 单 引号 里 为 线 型 和 颜色 参数 ,参数 可 
选 , 默 认 蓝 色 
例 2.5: 作出 函数 人 оч 的 图 像 。 
y =3s1nt 
MATLAB 程序 如 下 : 
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>> t= -2 *pi:0.1:2 * pi; 
>> х=2. ж соѕ(1) ;у =3. * sin(t); 
>> plot(x,y, r —. ') 
函数 图 像 如 图 2-4 所 示 。 
3 $ Н sg ° ть 
| | 
k k: 
| | 
N aal 
-24 
s= 
图 24 
=t 
Н. 由 于 y=/(*) 可 转变 为 参数 方程 | АШ, 
y=f(x) 也 可 用 此 法 作 图 。 
@ 作 分 段 函 数 的 图 形 ， 也 可 以 用 plot 命令 ， 其 MATLAB 
命令 格式 : 
>> xl =al:cl :bl ;х2 =a2:c2:b2;……… 
>>yl = 352 =; 


>> plot (x1 ,yl,S1',x2,y2, 52',.….) 


—x, х<0 
5 的 图 像 。 
x x >0 


例 2. 6 : 作出 y= 


MATLAB 程序 如 下 : 
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>>xl = -2:0.1:0 ;x2 =0:0.1:2 ; 
>>yl = -xl ;y2 =x2. 2 ; 

>> plot(x1 ,yl ,x2,y2,r-') 
函数 图 像 如 图 2-5 所 示 。 


4 T T T T T 





35 j} 
引 аа 


25| Ў = 











图 2-5 


(3 作 一 些 数据 点 的 散 点 图 : 

>>х=[... | ;y=[... ] 

>> plot(x,y,'*') 

例 2.7: 已 知 表 2-3 组 数据 是 黄河 小 浪 底 在 24 个 不 同时 
间 的 调 沙 量 


























表 2-3 (单位 : t) 
时 点 1 2 3 4 5 6 
调 沙 量 57.6 114 157.5 187 207 235.2 
时 点 7 8 9 10 11 12 
调 沙 量 250 265. 2 286. 2 302. 4 312. 8 307. 4 
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( 续 ) 
时 点 13 14 15 16 17 18 
调 沙 量 306. 8 300 271.4 231 160 111 
时 点 19 20 21 22 23 24 
调 沙 量 91 54 45.5 30 8 4.5 























根据 试验 数据 建立 数学 模型 ， 用 拟 合 的 方法 得 出 任意 时 
刻 排 沙 量 的 变化 关系 。 

> х=[123456789 10 11 12 13 14 15 16 17 18 19 
20 21 22 23 24]; 

>> у= [57.6 114 157. 5 187 207 235. 2 250 265.2 286. 2 
302. 4 312. 8 307. 4 306. 8 300 271.4 231 160 111 91 54 45. 5 














308 4.5]; 

>> plot (х, у, '*') 

散 点 图 如 图 2-6 所 示 。 
350 
300[- К. "еж а 
250! Е i 
200} 
a а Š 
100} F 
50t ^ * 
K: 5 10 15 wawa 
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在 2.5 节 将 讲 到 多 项 式 回 归 ， 如 果 想 将 此 曲线 拟 合并 在 
同一 画面 作 图 ， 则 需 加 如 下 程序 : 
>> polyfit(x, y, 3) % 拟 合 数据 


ans 二 





0.0798 -5.0661 75. 6607 -3l.5583 
>> hold оп 
>> plot (' 0. 0798 x х^3 - 5. 0661 * х^2 + 75. 6607 * х - 
31. 5583 ',[0, 25], ' =) 
则 得 如 下 散 点 图 和 拟 合 图 像 (图 2-7) : 


350 





300} p S S 
250 


f 
200} / 











3. 其 他 绘图 命令 

QD 符号 函数 的 简易 绘图 函数 ezplot : 

MATLAB 命令 格式 : 

>> ezplot('f(x)) Ф 默认 横 坐 标 范围 [ -2 * pi,2 * pi ] 
或 
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>> syms x y 


>> y=f(x); 
>> ezplot( у) 90 КУЛА МАУС -2 = pi ,2 * pi] 


当然 ,ezplot(f,[xmin,xmax]) 可 以 使 用 输入 参数 来 代替 默认 
横 坐 标 范 围 [ -2* pi,2* pi], 
例 2.8 画 出 函数 y=tanx 的 图 形 。 
MATLAB 程序 如 下 : 
>> ezplot ('tan (x) ') 


函数 的 图 形 如 图 2-8 所 示 。 























tan(x) 
| | 
4r | | | 
| | | | 
a | j J / 
S P 4 / 
Ta Pa Pa y 
: ⁄ ; z и A 
2 { / / | 
| | | | 
м Ка 
$ 4 2 0 2 4 6 
图 2-8 


H х“ + y“ — 8< -107° +16 =0 的 图 形 。 





例 2.9: (ЕН 
MATLAB 程序 如 下 : 


>> syms x y 
>> Е=х4 +y4 -8 ж х^ -10 *y2 +16; 





>> ezplot( F) 
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函数 的 图 形 如 图 2-9 所 示 。 


x*+y*-8 х2-10 у2+16 = 0 














ФРАЗ E pq 2 fplot: 

fplot 与 plot 命令 相似 ， 其 中 fplot(fun, lims) 绘制 由 字符 
HE fun 指定 函数 名 的 函数 在 x 轴 区 间 为 lims = | хи, xmax | 
的 函数 图 。 知 lims =[xmin，xmax，ymin，ymax] M) y 轴 也 
被 限制 。fun 必须 为 М 文件 的 函数 名 或 对 变量 x 的 可 执行 字 
符 串 ， 此 字符 串 被 送 入 函数 eval 后 被 执行 。 函 数 fan(x) 必 
须要 返回 针对 向 量 x 的 每 一 元 素 结果 的 向 量 。 

x+], x<1 
812.10: ш/а) =) ，1 Я. 
х 
ВЕ: (1) 首先 用 M 文件 funl. m ЕХ) F: 


function у =funl (х); 














Их <1 
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у=х+1; 
else 
у=1 +1. Zx; 
епа 
(2) fE matlab 命令 窗口 输入 : 
fplot('fun1',[ -3,3]) 
就 可 画 出 函数 fx) 的 图 形 。 
这 里 也 可 以 使 用 匿名 函数 ， 编 写 程序 如 下 : 
funl =@ (x) (х+1) ж (х<1) +(1+1ИХ) ж (x> =1); 
fplot(funl,[ -3,3]) 


2.5 多 项 式 回 归 


多 项 式 的 回归 就 是 将 两 组 数据 用 多 项 式 曲 线 进行 拟 合 ， 
其 程序 为 : 
>>х= [... |;у=[... 1]; % Л х ТУ 的 数据 
>> polyfit(x,y,n) % n 为 多 项 式 的 次 数 

例 2.11: 某 公司 2009 年 前 6 个 月 的 销售 收入 R( 万 元 ) 
( 见 表 24) ,根据 该 表 的 数据 ,在 Excel 中 发 现 , 该 数据 的 走势 
大 体 接近 三 次 函数 ,请 用 MATLAB 拟 合 得 出 该 曲线 ,其 中 ,t= 
0 对 应 1 月 。 假 车 这 个 趋势 继续 下 去 , 问 2009 年 11 HBA 
司 的 销售 收入 是 多 少 ? 


















































销售 收入 32.7 36.5 41.1 51.1 66.5 78 





解 : 曲 线 的 拟 合 程序 如 下 : 
>>1=[0,1,2,3,4,5], В = [32.7,36.5,41.4,51. 5, 
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66.5,78 ] 
>> polyfit(t,R,3) 
Ans = 
-0.2065 2.9861 -0.6963 33. 1444 
所 以 ,R(t) = -0. 2065г +2. 9861 -0. 6963t +33. 1444 
И 月 销售 收入 实际 上 对 应 的 是 上 = 10 时 函数 的 值 。 
R(10) = - 0. 2065 x 10° +2. 9861 x 102 — 








(T; 





0. 6963 x 10 +33. 1444 
=118( 万 元 ) 
例 2.12: 表 2-5 所 列 数据 是 某 公 司 网 上 电子 商务 每 年 收 
入 (单位 : 百 万 元 ) 。 试 利用 MATLAB 建立 一 个 3 次 多 项 式 
回归 模型 y=y(t) ,使 1=0 对 应 1997, 在 [0,， 4] х[0, 20] 
范围 内 作出 y =y(b 的 图 形 ， 并 计算 上 =0，1，2，3，4 时 收 
A y=y(t) ИВ 






































表 2-5 
年 份 1997 1998 1999 2000 2001 
收入 2.4 5 8 12 17. 4 


解 : (1) 首先 对 数据 进行 拟 合 ，MATLAB 程序 为 : 
>> t=[0,1,2,3,4],R=[2.4,5,8,12,17.4] 
>> polyfit(t,R,3) 
Ans = 
0.0833 -0.0286 2.5310 2. 4029 
y(t) =0. 0833г — 0. 0286 +2. 5310: +2. 4029 
(2) 在 指定 区 域 作 出 图 像 ， 在 MATLAB 窗口 输入 如 下 


命令 : 
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>> axis(0 ,4,0,20 ) 

>>t=0 :0.2 :4 

>> у =0.0833. * (t.3) -0.0286. * (2) +2. 5310. * t +2. 4029 
>> plot(t,y) 

于 是 ,得 到 如 图 2-10 所 示 的 回归 模型 图 像 。 

















д 
16| " 
74 
РА 
14 Р 
we 
12 РА 
i 
Pa 
10| g 
P 

8 пай? 
外 
2 А 

0 05 1 15 2 25 3 35 4 

图 2-10 


(3) 回归 验证 模型 ,在 MATLAB 窗口 输入 如 下 命令 : 
>> syms t 
>> у=0. 0833. ж (1 3) -0. 0286. * (t2) +2. 5310. * t +2. 4029 
>> subs (y,t,[0,1,2,3,4]) 
>> ans = 


2.4029 4.9886 8.0169 11.9876 17.4005 
Z 训 w 


1. 已 知 函数 /(x) =4x -5 5+1, RSO) WE 
2. АЯ (х) = x° — x +3. 45sinx -2cosx， 求 f(3) 的 值 。 
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3. 画 出 函数 y = tanx 的 图 形 。 






































4. 1] 0 <х<8 绘制 y=2e "sinx 的 图 形 。 
х+1, х<1 
5. Ш /(х) = 1 的 图 形 。 
1+—, х21 
x 











6. 画 出 y=sinx， y=sin(x+ Э +2, у= соѕх 的 对 比 图 。 
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3.1 极限 的 计算 


在 MATLAB 符号 工具 中 求 极限 的 指令 是 limit， 其 MAT- 
LAB 命令 格式 : 

>>syms x 

>>y=f(x); % 输 入 y 的 表达 式 

>>limit(y, x, a) % RINK АЖ f = x—a 时 极限 

说 明 : 

(Dlimit(f，a) : 表示 求 了 中 的 自 变量 (系统 默认 的 自 变 
EH x) Та 时 的 极限 。 

(limit(f) : 表示 求 /中 的 自 变 量 趋 于 0 时 的 极限 。 

(Зай (Е, x, a,，'left') : RIR S ха 时 的 左 极限 。 

Olimit(f, x, a 'right'): 表示 求 f 当 ха 时 的 右 极限 。 

(limit(f，x，inf) : 表示 求 当 x 一 % 时 的 极限 。 

ПЕ 

"Эхе “ln(l1 -2x) 

解 : MATLAB 程序 如 下 : 

>> syms X 

>> y=sin(2 xx) *tan(x°2)/((x”) x exp( -2 *x) * 
log(1 -2 *x)); 

>> limit(y) 














ans = 
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x a2 
= ° 
解 : MATLAB 程序 如 下 : 
>> syms x 
>> у= (1 +х/(х^2 +1) )^(2 жх); 
>> limit (у,х, — inf) 





例 3.2: 计算 lim [1+ 


ans 二 





ехр(2) 
>> double (ans) % 双 精度 型 数据 类 型 (精确 到 小 
数 点 后 第 四 位 ) 
ans = 
7. 3891 


3.2 导数 的 计算 


在 MATLAB 符号 工具 中 求 导 数 的 指令 是 diff( y), H 
MATLAB 命令 格式 : 





>>syms x 

>>y=f(x) ; 

>> diff(y, x, n) % 表示 求 函数 y 的 n 阶 导数 (n 必 
须 为 正 整 数 ) 

说 明 . 


(Ddiff( y) : 表示 diff(y， X, 1) ‚ 系统 默认 值 。 

@ 上 述 方法 得 到 的 结果 是 一 个 MATLAB 符号 表达 式 ， 
如 果 想 让 结果 能 如 通常 的 数学 表达 式 一 样 好 看 一 点 ， 则 可 增 
加 一 句 : 


>> pretty (ans) 
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若 与 计算 相关 的 式 子 中 有 其 他 字母 ， 则 必须 在 第 一 句 的 
x 之 后 全 部 列 出 ， 且 用 去 号 或 空格 分 隔 。 

OME yO Æ x =a 的 值 ， 只 需 修 改 第 三 句 为 : 

>>d=diff (y,x,n) ,subs (d,x,a) 

例 3.3: y=xsin2x-e*, 3K y'o 

解 : MATLAB 程序 如 下 : 


>> syms X 











>> y =x *sin(2 xx) —exp(3 *x) ; 
>> diff( y) 

ans = 

sin(2*x) — 3 жехр(3 жх) + 2 *x *cos(2 жх) 
例 3.4:y = x lnx, фу", 

解 :MATLAB 程序 如 下 : 

>> syms x 

>> y=(x3)*log(x) ; 

>> а (у, х,3) 

ans = 

бх Іов(х) + 11 

>> ргейу (апѕ) 

6 log(x) + 11 

01] 3.5: у = е*соѕ3х, Жу"(0). 
Е: MATLAB 程序 如 下 : 


>>syms x 





>>у=ехр(х) * соѕ(3 *x) ; 
>> A =аШ(у,х,3) ; 
>> subs (А,х,0) 


ans = 
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例 3.6:  y=xarctanx, > y"( a). 
解 : MATLAB 程序 如 下 : 
>>syms x 
>>y =x * atan( x) ; 
> b =АШ(у,х,2) ; 
>> subs (b,x,a) 
ans = 
2/(1 +a?) -2 ж а^2/(1 +a 2)2 
>> pretty (ans) 


3.3 积分 的 计算 


在 MATLAB 符号 工具 中 求 积分 (intergration ) 的 指令 是 








int， 其 MATLAB 命令 格式 : 
>>syms x 
>>y=f(x) ; 
>> int(y,x) % 求 y 的 不 定 积分 


>>int(y,x,a,b) ”% 求 y 在 [a,b] 上 的 定 积分 

其 中 ， 当 积分 上 下 限 或 被 积 函 数 中 还 含有 除 x 以 外 的 其 
他 字母 时 ， 必 须 在 第 一 句 中 的 x 之 后 列 出 ， 且 用 逗号 或 空格 
分 隔 。 

如 果 想 让 不 定 积 分 的 计算 结果 好 看 一 点 ， 则 可 修改 第 三 
№]: 

>>l=int (у,х) ,ргейу(1) 

定 积分 的 计算 结果 是 一 个 符号 ， 要 转换 为 数值 需 再 输入 
double (ans) ， 或 修改 第 四 句 为 : 

> J=in (y,x,a,b) ,double (J) 
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913.7: у=—^ НЯ [yde 和 [ yde 的 MAT- 
= Гра 和 | 


1 +x 





LAB 计算 程序 。 





解 : MATLAB 程序 如 下 : 

>>syms x 

>>y= (x3)/sqrt(1 +x°2); 

>> int(y,x) 

ans = 

1/3 ж х^2 ж (1 +x°2)^(1/2) -2/3 x (1 + х^2)^(1/2) 
>> шКу,х,0,1) 

>> апѕ = 

1/3 ж 2^(1/2) +2/3 

>> double ( ans) 


ans = 


0. 1953 








2 1 
例 3.8: y= 一 ， 写 出 计算 | ydx 的 MATLAB 计算 程 
1 +x ° 


2? 


f: MATLAB 程序 如 下 : 

>> syms x 

>>y= (2x*x)/ (1 +x2); 
>>] =шКу,х,0,1), double(J) 
J =log(2) ,ans =0. 6931 


913.9; у= (1 +27), ЗНАЯ | yde 的 MATLAB 计 


算 程 序 。 


解 : MATLAB 程序 如 下 : 
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>>syms ха b 
>>y=log(1 +x2) 
>> int(y,x,a,b) 
ans =b * log(1 +b°2) -2 xb +2 * atan(b) -а* log(1 + 
x2) +2 * a —2 ж atan( a) 
3.4 求 函 数 极 值 、 最 值 
灵活 地 运用 MATLAB 的 计算 功能 ， 再 根据 极 值 理论 可 
以 很 容易 求 得 函数 的 极 值 。 
例 3.10: SRŽ y =x +252 -5х +1 ВИН. 
解 题 思路 一 : 用 diff 命令 求 函 数 导 数 ; 
再 用 solve 命令 求 驻 点 ; 


再 用 plot 命令 绘制 函数 图 像 ， 根 据 图 像 判 断 驻 点 是 
否 为 极 值 点 。 


MATLAB 程序 如 下 : 
>> syms xy 





>> y=x3+2*x2-5*x+1; 
>> dy = (у) % >K 
dy = 
3 жх2 +4 жх -5 
>> solve( dy) % 求 驻 点 
ans = 
19^(1/2)/3 -2⁄3 
_19°(1/2)/3 -2⁄3 
>> x = double( ans) 
x= 


0. 7863 





% 驻 点 近似 双 精 度 取 值 


第 3 章 微 积分 的 MATLAB 计算 .33. 





—2. 1196 
>> ezplot(y,[ -4,2]) % 画图 
图 像 显 示 如 图 3-1 所 示 。 


e ERTA 

















根据 图 像 判 断 ，x = -2.1196 是 极 大 值 点 ，x = 0. 7863 
是 极 小 值 点 ， 可 以 代入 原 式 求 极 大 极 小 值 。 
>> subs(y,x, 0. 78633) 








ans = 
– 1. 2088 % 极 小 值 
>> subs(y,x, –2. 1196) 
ans = 
11. 0607 % 极 大 值 


解 题 思 路 二 : 用 df 命令 求 函数 导数 ; 
再 用 solve 命令 求 驻 点 ; 
再 用 df 命令 求 函数 二 阶 导数 ; 
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再 用 subs 将 驻 点 代入 二 阶 导数 ， 根 据 二 阶 导 数 的 
符号 判别 是 否 为 极 值 点 。 
MATLAB 程序 如 下 : 
>> syms x y 





>> y=x3+2*x2-5*x+1; 
>> dy = (у) % >K 
dy = 
3*x2+4*x-5 
>> solve( dy) 9% RIER 
ans = 
19^(1/2)/3 -2/3 
— 19^(1/2)/3 -2/3 
>> х = double( ans ) % 驻 点 近似 双 精 度 取 值 
x= 
0. 7863 
–2. 1196 
>> ddy = diff( dy) % 求 二 阶 导 
ddy = 
6*x+4 
>> subs(ddy, x, 0. 78633) 
ans = 
8.7180 
>> subs(ddy, х, -2. 1196) 
ans = 


-8.7176 


因为 x = -2.119 代入 二 阶 导 的 值 为 负 ， 所 以 x = 
-2. 1196 是 极 大 值 点 ，x = 0.7863 代入 二 阶 导 的 值 为 正 ， 所 
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以 x=0.7863 是 极 小 值 点 。 
Зх? +4х +4 
例 3.11: у= урин. 
x +x+] 
解 : MATLAB 程序 如 下 : 
>>syms х y 


>>y= (3 ж х^2 +4 жх+4)/( x2+x+1); 


>> ду =аШ(у); % >+ 

>> xz = solve ( ду) % 求 驻 点 
xz = 

[0] [ -2] 

>> d2y = diff(y,2); % 求 二 阶 导 


>>zl =subs( d2y,x,0) % 驻 点 x=0 代入 二 阶 导 

21 三 

-2 

>> 22 =subs (42у,х, -2) EA x= -2 代入 二 阶 导 

72 = 

2⁄9 

可 知 在 x =0 处 二 阶 导 数 的 值 为 z, = -2， 小 于 0， 函 数 
有 极 大 值 ; 在 = -2 处 二 阶 导数 的 值 为 2 =2/9, КТО, 
函数 有 极 小 值 。 如 果 需 要 ， 可 顺便 求 出 极 值 点 处 的 函数 值 : 

>>y, =subs (у,х,0) 





yq = 

4 

>> у, =subs (y,x, -2) 
JYS 

















事实 上 ， 如 果 知 道 了 一 个 函数 的 图 形 ， 则 它 的 极 值 情况 
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和 许多 其 他 特性 是 一 目 了 然 的 。 而 借助 MATLAB 的 作 图 功 
能 ， 我 们 很 容易 做 到 这 一 点 (图 3-2) 。 
>> ezplot(y) 


(3 X2+4 X+4)MX2+X+1) 














图 32 


在 MATLAB 的 语言 中 ， 求 函数 在 给 定 区 间 上 的 最 小 值 
命令 是 : fminbnd， 调用 格式 如 下 : 

>> х =fminbnd(y,x1 ,x2) 

说 明 : у 是 函数 的 符号 表达 式 ; fminbnd 仅 用 于 求 函数 
的 最 小 值 点 ， 如 果 要 求 最 大 值 点 ， 可 先 将 函数 变 号 ， 求 得 最 
小 值 点 ， 即 得 到 所 求 函数 的 最 大 值 点 ;x1，x2 是 变量 x 的 
取 值 范围 。 

例 3.12: 求 函数 y=e“+ (x+1) 在 区 间 [ -3, 3] 内 
的 最 小 值 。 

解 : MATLAB 程序 如 下 : 

>> x=fminbnd('exp( -х) +(x+1)2',-3,3) 
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х = 
– 0. 3149 
>> y=subs('exp( =x) + (х+1)^2 ',х, —0. 3149) 
у= 
1. 8395 
所 以 最 小 值 为 1. 8395 
实 训 题 
1. 求 极限 : 
(1)1 sinx (2) li Еа 2 
о а | bn ( ст) 
(2х +1. _ 1 
(3) lim( ) (4) lim 一 
== \ 2х — 5-х 
1 
x 0sinx 
а [39+ гек 
з) ЕЎ Sh 





3. 求 函 数 的 导数 : (1) y=xe” 


(3) y =e ™tan2x 
4. 求 不 定 积分 : (1) [e'sin?wdx 


G) | 于 


г COSX 








(2) у= cos ух 
Inx 
(4) y= = 


2 
x 


(2) [+ х) хах 


+ siny _ 


e dx 


5. 求 定 积分 : (1) [G< — 5) агссозхах (2) | МТ - sin2xdx 


6. 讨论 函数 Kx) = e “ВИН. 
7. ЖЖ ==" -2x? +5 在 区 间 [ -2，2] 内 的 最 小 值 。 





第 4 章 用 MATLAB 解 常 
微分 方程 ( 组 ) 


用 MATLAB 软件 解 常 微分 方程 (组 ) 有 两 种 方式 ， 一 种 
是 符号 求解 ， 这 样 求 得 的 通常 是 方程 的 解析 解 ， 男 一 种 是 数 
值 求解 ， 这 是 一 种 通过 近似 计算 而 得 到 近似 解 的 方式 。 能 采 
用 符号 求解 的 常 微分 方程 (组 ) 的 范围 十 分 有 限 ， 通常 仅 限 
于 常 系数 的 线性 常 微分 方程 (组 )， 以 及 少数 的 变 系数 微分 
方程 (组 )。 对 于 更 广泛 的 非 线 性 微分 方程 (组 )， 一 般 很 难 
得 到 其 解析 解 ， 通 常 采 用 的 是 数值 求解 方式 。 微 积分 学 习 过 
程 中 ， 所 遇 到 的 微分 方程 一 般 都 能 求 得 其 解析 解 ， 因 此 ， 下 
面 仅 考 察 党 微分 方程 (组 ) 的 符号 求解 。 

常 微分 方程 (组 ) 常 用 的 符号 求解 命令 是 dsolve， 求 解 的 
MATLAB 命令 格式 如 下 : 

>> [yl, ..., yu] =dsolve( 方 程 1，:…… ， 方程 n'，' 条 
件 1，…， 条件 m'，'var') 

程序 输入 注意 事项 如 下 : 

Qm. n 必须 是 有 限 数 ， 不 能 为 字母 ; 

Фу, ... , yn 是 待 求 的 未 知 函数 变量 , n ТН, af 
用 [yj 或 y НВ, ..., т]; 

@ 方 程 中 的 y" 输 入 为 Dny， 且 DD 必须 大 写 ; 

@ 条 件 中 y(a) =b My” (с) =d 的 输入 分 别 为 yY(a) =b 
和 Dky(c) =d(D 必须 大 写 ); 

OREL, oos ， 条 件 m 默认 其 程序 是 在 求 微分 方程 





















































k=’ 
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(组 ) 的 通 解 ; 
Ovar 为 指定 的 自 变 量 字 母 ， 默 认 所 求 得 的 解 的 表达 式 
中 的 自 变量 为 i。 
W, Ry" -y'= sinx 的 通 解 ， 在 MATLAB 命令 窗口 输入 
>> dsolve('D2y - Dy = sin(x)') 
ans = 
Cl —sin(x) + C2 * exp(t) -t * sin(x) 
不 难看 出 ,在 求解 过 程 中 ,sinx 被 当做 了 常数 ,这 显然 不 
是 我 们 所 要 求 的 结果 。 指 定 var 为 x, 在 MATLAB 命令 窗口 
输入 





Ш 








>> dsolve('D2y – Dy =sin(x)',' x ') 

ans = 

СЗ +cos(x)⁄/2 – зт(х)/2 + C4 * ехр(х) 
才 是 正确 的 结 
例 4.1: 求 方程 yy + (у)? =0 的 通 解 。 
解 : MATLAB 程序 如 下 : 
>> y=dsolve('y * D2y – (Оу)^2 =0','x') 
y= 








C3 

Cl x ехр( C2 * x) 

即 方程 的 通 解 为 :y = Ce” y =C. 

注意 :y = C, 只 含有 一 个 任意 常数 , 称 为 奇 解 。 因 此 , 方 
ERNER y = Сџе ЕА 

例 4.2: 求 方程 2xyy' +y =2 的 通 解 。 

解 :MATLAB 程序 如 下 : 

>> у =dsolve('2 x (x2) *y*Dy+y2 =2','х') 


Ira 








у = 
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2^(1/2) 
-2^(1/2) 
(exp( Cl +1⁄x) +2)^(1/2) 
— (exp( C1 +1⁄x) +2) (1⁄2) 
即 方程 的 通 解 为 : У’ =2 + Се“ 
例 4.3: REE y +у -2y =x WERF yO) =4， 
y (0) =1 的 特 解 。 
解 : MATLAB 程序 如 下 : 
>> y=dsolve('D2y + Dy -2 жу=х','у(0) =4,Dy(0) = 
Lao 
y= 
11/(12 * ехр(2 *x)) – х/2 + (10 * ехр(х) )/3 – 1/4 
БОРТ 1 1 10, п 
即 方程 的 特 解 为 : у 4 = и 3 e tihe 
им 


1. 求 微分 方程 y =1+у’ 的 通 解 。 
2. 求 常 微分 方程 v +y+ (x —2y)y' =0 的 通 解 。 


一 2x 




















第 $ 章 和 矩阵 运算 的 MATLAB 实现 


5.1 ERRA 


和 矩阵 的 输入 程序 格式 为 : 

>> 和 矩阵 名 =[ 第 1 行 ;…; 第 n 行 ] 

具体 规则 如 下 : 

@ 和 矩阵 的 元 素 置 于 方 括号 之 中 ; 

@) 行 与 行 之 间 用 分 号 或 回 车 符 分 隔 ; 

多 同行 元 素 之 间 用 逗号 或 空格 分 隔 。 

MATLAB 软件 关于 和 矩阵 的 输入 有 数字 和 矩阵 和 符号 矩阵 之 
分 ， 数 字 和 矩阵 按 规则 直接 输入 ， 符 号 矩阵 在 按 规则 输入 前 ， 
必须 先 将 矩阵 中 的 字母 定义 为 符号 。 

















1 2 3 
例 5.1: 写 出 输入 矩阵 4 =|4 5 6 | 的 MATLAB fy 
7 8 9 








令 。 


解 : 矩阵 A 的 MATLAB 程序 可 以 有 如 下 三 种 输入 法 : 
Ф>>А=[1,2,3;4,5,6;7,8,9 | 
(2) >>А=[123;456;789] 
(3) >> А=[123 
456 
789] 








运行 结 


. 42. MATLAB 基础 与 应 用 





5.2 











A=1 2 3 
4 5 6 
7 8 9 
1 2 3 
例 5.2: 写 出 输入 矩阵 4=| a 0 i MATLAB 命 
-1 се 2 
解 : MATLAB 程序 如 下 : 
>>syms a b с 
>A=|123;a0b;-1 c2] 
A=|1 2 3] 
[a,0 ,bj 
[-1,е ,2] 
FE JA: JJ H: 
КА, ВАНИЕ 
>>А+В 多 矩阵 的 加 法 
>A-B 儿 和 矩阵 的 减法 
>>А*В % ҖЕ ЕЕ 
>> А. * В % 数 组 的 乘法 一 一 对 应 行列 位 置 
上 的 元 素 相 乘 
>> A\B % A'B 
>> А/В % AB `' 
>>A' % A 的 转 置 
>> inv( A) % А HJ A! 
>> det( A,) %n WIERE А, 的 行列 式 


>> rref( A) % 将 A 化 为 行 阶梯 形 
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T 


>>[P,Q] =eig(A) 9% n ИННЫ. Q 和 特征 














向 量 了 
1 2 3 

例 $.3: 求 矩阵 4=|4 2 1 的 行列 式 值 .特征 值 和 特 
2 1 3 

征 向 量 。 

解 :MATLAB 程序 如 下 : 

中 求 行列 式 的 值 : 

>> А = [1,2,3;4,2,1;2,1,3]; 

>> det( A) 

Ans = 


-15 
@) 求 特征 向 量 和 特征 值 : 
> A=[1,2,3;4,2,1;2,1,3] 
>> [P,Q] =eig (А) 
P= -0. 5495 -0.6894 -0.0186 
-0.6420 0.7102 -0. 8331 
-0.5347 0.1428 0. 5528 








Q =6. 2561 0 0 
0 -1.6818 0 
0 0 1.4257 





即 特征 值 Q1 = 6. 2561, 02 = -1. 6818, 03 = 1. 4257, 
其 对 应 的 特征 向 量 分 别 为 : 
-0. 5495 -0. 6894 -0.0186 
-0.6420 |, | 0.7102 |, | -0. 8331 


– 0. 5347 0. 1428 0. 5528 
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#1 5.4: КЕША = 





l 2 3 1 2 

42153) 
阶梯 矩阵 。 

2 1 3 1 2 

2 4 6 2 4 


解 : MATLAB 程序 如 下 : 
>> A=[1,2,3,1,2;4,2,1,5,3;2,1,3,1,2 ;2,4,6,2,4] 


>> rref( A) 


ans = 1 


1. ЕИ А 





0 


0 
1 


0.9333 0.4667 
0. 9333 0.4667 
–0. 6000 0. 2000 
0 0 
Z 训 题 
123 
45 6 
7 8 9 
3 4 2 5 
‚ В= ‚ Ж(1) A+B; (2) А- 
2 2 -3 4 


0 Ë waa 
) B=| 2 l| 3ROO1)2A; (2) АВ, 
-2 -5 0 


-3 | Ж (1) АП; (2) А’; (3) A; (4) Ж 
2 


第 5 章 ЖӘЕ MATLAB 实现 


. 45. 





JERE A ЙК, 
5. ЖА MB 是 满足 关系 式 AB = A +2B BJ AE EE, 











4 2 3 
1 1 | ЖИ: B, 
3 


2 





` 


ГЕНА = 





第 6 章 用 MATLAB 求解 
线性 方程 组 


6.1 求解 单个 代数 方程 

格式 : solve (f,t) 

功能 : 对 变量 1 解 方 程 f=0, t 默认 为 * 或 最 接近 字母 x 
的 符号 变量 。 

如 : 求解 一 元 二 次 方程 y=ax + bx +c 的 实 根 。 

> symsa b c x 





>> f=a*x2+b*x+c; 
>> solve (f,x) 
ans = 
[1/2/а* ( -b+ (2 -4 жаже)^ (1/2))] 
[1/2/а ж ( -b - (Ь2 -4 xa x c)^ (1/2))] 
如 果 想 对 非 默认 x 变量 求解 ，solve 必须 指定 变量 。 
>> solve( 'axx2+bxx+c','b') % 求解 6 
ans = 
— (аж х^2 +е)/х 
带 有 等 号 的 符号 方程 也 可 以 求解 ,如 : 


>> f=solve( 'cos(x) =sin(x) ') 








А1 
>> t=solve( 'tan(2 ж х) =sin(x)') 
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一 
lI 


[ 0] 
[acos(1⁄2 +1⁄2 x 3^(1/2)) | 
[асоз( 1/2 – 1/2 ж 3^( 1/2) ) ] 

并 得 到 数值 解 。 

>> numerie(f)  & 数量 化 


0.7854 








>> numeric(t) 
ans = 
0 
0 +0. 83141 
1. 9455 


6.2 求解 线性 方程 组 


线性 方程 组 的 求解 分 为 两 类 :一 类 是 方程 组 求 唯一 解 或 
求 特 解 ， 另 一 类 是 方程 组 求 无 穷 解 ( 即 通 解 ) 。 可 以 通过 系 
数 和 矩阵 的 秩 来 判断 : 

若 系 数 和 矩阵 的 秩 "=n(n 为 方程 组 中 未 知 变量 的 个 数 ) , 
则 有 唯一 解 ; 

若 系数 矩阵 的 秩 r<n， 则 可 能 有 无 穷 解 ; 

线性 方程 组 的 无 穷 解 = 对 应 齐 次 方程 组 的 通 解 + 非 齐 次 
方程 组 的 一 个 特 解 ， 其 特 解 的 求法 属于 解 的 第 一 类 问题 ， 通 
解 部 分 属 第 二 类 问题 。 

І. 求 线性 方程 组 的 唯一 解 或 特 解 (第 一 类 问题 ) 

利用 矩阵 除法 求 线性 方程 组 的 特 解 ( 或 一 个 解 ) 
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方程 : AX =b 
解法 : X=A\b 
5х; +бх, =] 
х, +5х, +6x; =0 
例 6.1， 求 方程 组 к +5, +64 =0 的 解 。 
хз +5x, +6x; =Ü 
Я: 三 | 
解 : MATLAB 程序 如 下 : 
>>A=[5 6 0 0 0 
15600 
01560 
00156 
0 0 0 1 5]; 
>>B=[10001]; 


>> КА =rank(A) % 求 秩 


>> X = 


А\В 


多 求解 


运行 后 结果 如 下 : 


R_A= 


5 


X= 


2. 2662 
—1. 7218 
1. 0571 
– 0. 5940 
0. 3188 





这 就 是 


方程 





组 的 解 。 
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或 用 函数 rref 求解 : 
> С= [А,В]  % 由 系数 矩阵 和 常数 列 构成 增 广 矩 阵 С 
>> R=rref(C) %#С 化 成 行 最 简 形 























R= 
1. 0000 0 0 0 0 2. 2662 
0 1.0000 0 0 О -1.7218 
0 0 1.0000 0 0 1. 0571 
0 0 0 1.0000 0 -0. 5940 
0 0 0 0 1.0000 0. 3188 





W R 的 最 后 一 列 元 素 就 是 所 求解 。 
м + м) 一 3x3 - x, =1 
例 6.2: 求 方程 组 ,3x, – х, - 3х; +4x, =4 的 一 个 特 解 。 
xi +5x, —9x, — 8х, =0 
解 : МАТГАВ 程序 如 下 : 
>>A=[11-3-1;3-1-34;15-9-8]; 
>>В=[1 4 0J; 
>>Х=А\В % 由 于 系数 矩阵 不 满 秩 ,该 解法 可 能 存在 误 
差 。 
X=[0 0 -0.5333 0.6000] (一 个 特 解 近 似 值 ) 。 
若 用 reef 求解 , 则 比较 精确 : 
>> A=[11-3-1;3-1-34;15-9-8]; 
B=[1 4 01; 
>> С=[А,В]; ЖАЛИ AE 
>> В =rref( С) 
В = 
1. 0000 О -1.5000 0. 7500 1. 2500 
0 1.0000 -1.5000 -1.7500 -0.2500 
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0 0 0 0 0 
由 此 得 解 向 量 导 = [1.2500 -0.2500 0 0]'( 一 个 特 
解 ) 。 
2. 求 线性 齐 次 方程 组 的 通 解 
TE MATLAB HF, РЖ null 用 来 求解 零 空间 ， 即 满足 4 
“站 =0 的 解 空间 ， 实 际 上 是 求 出 解 空间 的 一 组 基 ( 基础 解 
系 )。 
格式 = пи % Z 的 列 向 量 为 方程 组 的 正 交 规范 
Ж, MEZ xZ=I 
Z=null(A, 'r') %Z 的 列 向 量 是 方程 AX =0 的 有 理 
基 





х +25, +253 + х =0 
016.3: 求 方程 组 的 通 解 : (25 + x, -2x3 -2x4 =0 
х= х 一 4x3 - Зх, =0 
解 : MATLAB 程序 如 下 : 
>A=|1 2 2 1;2 1 -2 -2;1 -1 -4 -3]; 
>> format rat % 指定 有 理 式 格式 输出 
>>В=пи (А ,'г') 2 求解 空间 的 有 理 基 
运行 后 显示 结果 如 下 : 
B= 





1 
或 通过 行 最 简 形 得 到 基 : 
>> B=rref( A) 
B= 


第 6 竟 M MATLAB 求解 线性 方程 组 . S1 






























































1. 0000 О -2.0000 -1.6667 
0 1. 0000 2. 0000 1. 3333 
0 0 0 0 
即 可 写 出 其 基础 解 系 ( 与 上 面 结果 一 致 ) 。 
写 出 通 解 : 
syms kl k2 
Х=К! *В(:,1) +К2 *В(:,2) % 写 出 方程 组 的 通 解 
pretty(X ) % 让 通 解 表达 式 更 加 精美 
运行 后 结果 如 下 : 
X= 


[ 2xkl+5/3 * К] 

[ -2 = kl —4⁄3 * k2 ] 

[ к] 

[ k2] 

% 下 面 是 其 简化 形式 : 
[2kl +5/3k2 | 


[ -2kl -4/3k2] 


k1 








k2 


з. 求 非 齐 次 线性 方程 组 的 通 解 





非 齐 次 线性 方程 组 需要 先 判 断 方程 组 是 否 有 解 ， 若 有 


解 ， 再 去 求 通 解 。 
因此 ， 步 又 为 : 




















第 一 步 : HE AX =b ESA, FARETE. 
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第 二 步 : RAX =b H 

第 三 步 : 求 4X=0 的 通 解 。 

第 四 步 : AX = b 的 通 解 = 
个 特 解 。 


K = 285 + 33%; – 


一 个 特 解 。 





АХ =0 的 通 解 + AX = b 的 一 














x, =1 


例 6.4: 求解 方程 组 (3xi - м +5x3 一 3x4 =2 


2 + х, +2^. -2х. =3 


解 : MATLAB 程序 如 下 : 





在 MATLAB 中 建立 М 文件 如 下 : 


А=[1 
Ь= [1 

В=[ АБ]; 

п=4; 

В_А =rank( A) 

В_В =rank( B) 

format rat 

ИВ_А = =К_В&В_А = =n 
X=A\b 

elseif В_А = =В_В&В_А <n 
X=A\b 
C=null(A,'r') 


else X =' equition no solve ' 


-2 3 
2. 3]: 


end 


运行 后 结果 显示 : 





-1;3 -1 


5 -3;2 1 2 -2]; 


9 判断 有 唯一 解 


和 判断 有 无 穷 解 

% 求 特 解 

% K AX =0 的 基础 解 系 
Yo 判断 无 解 
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Х = 
equition по solve 
说 明 该 方程 组 无 角 
ee - м. =1 
例 6.5: 求 方程 组 的 通 解 : (3xi - x 一 3x3 +4x4 =4 
xi +5x, -9х. — 8х, =0 
解法 一 : 在 MATLAB 编辑 器 中 建立 M 文件 如 下 : 
А=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8]; 
Ь=[140}'; 
B=[Ab]; 
n=4; 
R_A =rank(A) 
R_B =rank(B) 
format rat 
И В_А = =В_В&В_А = =n 
X=A\b 
elseif КА = =В_В&В_А <n 
X=A\b 
C=null(A,'r') 


else X =' Equation has no solves ' 





end 
运行 后 结果 显示 为 : 
R_A = 
2 
R_B = 
2 
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Warning: Rank deficient, rank =2 tol = 8. 8373е - 
015. 
> In D:\Matlab\pujun\lx0723. т at line 11 
Х = 
0 
0 
-8/15 
3⁄5 
C= 
3⁄2 -3/4 
3⁄2 7/4 
1 0 
0 1 
3⁄2 -3/4 
所 以 原 方程 组 的 通 解 为 X= 厂 | ee) а 
0 1 
0 
0 
-8/15 
3⁄5 


解法 二 : 用 rref 求解 














A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8]; 
b=[140]; 

B=[Ab]; 

С =rref( B) 多 求 增 广 和 矩阵 的 行 最 简 形 ， 可 得 最 简 同 解 方 
程 组 。 
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运行 后 结果 显示 为 : 
C= 
1 0 -3/2 3⁄4 5/4 
0 1 -3/2 – 7/4 – 1/4 
0 0 0 0 0 
3/2 – 3/4 
3/2 7/4 
对 应 齐 次 方程 组 的 基础 解 系 为 : £, = | £ = g 
0 1 
5/4 
N | З – 1/4 е 
非 齐 次 方程 组 的 特 解 为 : 1” = Ó F 所 以 ， 原 方程 组 的 
0 
通 解 为 : Х= Кё, + k,ë, +n” o 
хим 


1. 求 方程 组 的 通 解 : 
x, +2x,) +2853 + x, =Ü 
2x, +X, – 2x, — 2x, =0 
x, — x, —4x, — 3x, = Ü 
2. 求 非 齐 次 方程 组 的 解 : 
2x; +4x, =11 
За := 5х5, = 3 
ху +22 =6 
25; +% =7 


3. 用 最 小 二 乘法 解 方程 组 : 

















. 56- 


MATLAB 基础 与 应 用 





. 求 方程 组 的 解 : 


X d 
хү 
X11 


хү 


求解 方程 组 : 


求解 方程 组 : 
2% 





д — Xs + x, = Ü 
— x, +x, —3x, = 1 


— x, —2x, +. = — 1/2 


+x, – 5x; + x, = 8 


xı — 3x, — 6x, = 9 


2х, 


Xl 


21 


3x, 





| 


线性 方程 组 : 


21 

















+ х, 


(1+À)x, +x, 


=y +2, = =5 


+4x — 7x, + 6x, = Ü 


-3 = =: 1 


— x, = 3х; +4x, =4 


+5x, —9x, — 8x, = 0 


+ x, =Ü 





(1 +À)x, 





х1 


| 





x, =3 


xi tx, +(1+À)x; ЕЛ 
Ил 取 何 值 时 ， 此 方程 组 有 








解 ? 











第 7 章 用 MATLAB 解 线性 
规划 问题 


MATLAB 只 能 解 minZ 类 的 线性 规划 问题 ， 对 于 maxZ 
类 的 线性 规划 问题 ， 可 令 f= -Z， 化 为 求 mnf 问题 。 此 外 ， 
不 等 式 约束 中 的 “三 ” 需 在 不 等 式 两 边 同 乘 以 -1， 将 其 化 为 
“和 "后 再 确定 和 矩阵 A 和 hb。 
通常 当 x; 为 非 负 连续 且 变 量 个 数 较 多 时 ， 可 用 MAT- 
ГАВ 解 线性 规划 问题 ， 其 程序 如 下 : 

Эре ое І % 输 入 目标 函数 的 系数 
和 矩阵 
эх А енем ] ,bp=[ ] % 输 入 不 等 式 约束 的 系 
数 和 矩阵 和 常数 项 列 向 


ЕН 


>> Aeq =[ ] ,peq=[ ] ШЛА АЖ 
AE RE FI B£ ЗС A) [n] 


E. 


里 
ьар] 。% 输 入 决策 变量 的 下 、 

上 界 列 向 量 
>>[x,fval | =linprog (f,A,b,Aeq ,beq,lb,ub) 

其 中 ,(f,A,b,Aeq ,peq,lb,ub) 中 元 素 的 次 序 不 能 颠 
倒 , 若 某 个 元 素 省 略 , 则 该 元 素 用 [ ] 人 代替 。 但 当 它 后 面 的 所 
有 元 素 也 省 略 时 该 元 素 也 可 省 略 。 

例 7.1: 解 线 性 规划 问题 : 
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maxZ =0. бх, +0. 4x, +0. 32х, +0. 48x, 
хм = - м, == ü 
X, + X, — x, 20 
Ке ху +, AS + = 1 
ху "м. * д, 20 
解 : 令 /= -Z， 并 将 约束 条 件 中 的 “之 "化 为 <"， 从 
而 MATLAB 程序 如 下 : 
>>f= -[0.6,0.4,0.32,0.48] 
>>A=[1,-1,-1,-1;0,-1,-1,1 ],b=[0;0] 
>> Аеа = [1,1,1,1 ],beq=[ 1] 
>>12=[0;0;0;0] 
>> [x,fva | = linprog (Ё, А,Ь, Аед, Бед, №) 








х = 0. 5000 
0. 2500 
0. 0000 
0. 2500 
fval = — 0. 5200 
注意 : ub 省 略 ， 但 后 面 没 有 元 素 了 ， 因 此 ， 不 必用 
“ГРА. 


817.2: 某 工 厂 生产 A、B 两 种 产品 ， 见 表 7-1。 已 知 生 
产 1kg 的 A 产品 ， 耗 煤 9t、 耗 电 4kW . h、 用 劳力 3 个 工 、 
利润 为 500 元 ; 生产 1kg 的 B 产品 ， 耗 煤 4t、 耗 电 5kW +h, 
用 劳力 10 个 工 、 利 润 为 900 元 。 若 可 提供 的 资源 数量 为 : 
炬 360t、 电 200kW .h、 劳 力 300 个 工 。 问 应 如 何 安排 生产 
量 ， 能 使 生产 A、B 产品 的 总 利润 为 最 大 ? 

解 : 设 A、B 两 种 产品 的 生产 量 分 别 为 x, 和 *,kg， 则 
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其 耗 煤 量 为 9x, +40, (1); 耗 电量 为 4x, +5х, (kW : h); 用 
劳力 为 3x, +10x,( 个 工 ) ， 总 利润 为 Z =500x, +900x, 元 。 









































表 7-1 
产品 
A B 可 供 资 源 
消耗 
煤 /t 9 4 360 

B/kW - h 4 5 200 
劳力 /个 工 3 10 300 
利润 /元 500 900 














因为 问题 的 目标 是 使 总 利润 为 最 大 ,约束 条 件 是 只 能 提 
供 煤 360t、 电 力 200kW - h 和 劳动 力 300 个 工 ， 转 换 成 数学 
式 子 即 为 线性 规划 模型 : 

тахй = 500х, + 900x; 
9х, +4x, <360 
4x, +5x, <200 
3x, + 10x, <300 


























Wa 三 0 
用 矩阵 表示 为 : 
maxZ = (500, 900) X 
AX<b 
д АЯ x, 20 
9 4 360 X, 
其 中 4=|4 5 |, В=| 200 |, X= 2 
3 10 300 x 


MATLAB 程序 如 下 : 
>>f= - [500,900 ] 
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>>A=[9,4;4,5 ;3,10 ],[360 ; 200 ; 300 | 
>>Ъ=[0;0 ] 
>> [х, Ва! ] =linprog ( f£,A,b,[ ],[ ],lb ) 
х =20. 0000 
24. 0000 
fval = -3. 1600е +0. 04 
注意 : aeq 和 beq 2786, (H lb FE, Kk, MALI” 


和 “[ ] "代替 。 





实 训 题 
1. 求解 线性 规划 问题 : 
minz = —2x, —%, + Xs , 


X, +X, +2x, = 6 





ху +4x, — x, <4 
2x, —2x, + x, < 12 
х 20, х,20, x, <5 
2. 求 下 面 的 规划 问题 : 
min —Sx—4x,— бх; 
x, — x; + x, <20 


sub. to 
3x, +2x, +4x, <42 





3x, +2x, <30 


Ол, 0<x,, 0<x, 


第 8 章 概率 问题 的 MATLAB 计算 


8.1 常见 分 布 的 MATLAB 各 称 ( 表 8-1) 











表 8-1 
分 布 名 称 MATLAB 名 称 
二 项 分 布 bino 
泊 松 分 布 poiss 
离散 型 
几何 分 布 geo 
超 几 何 分 布 hygo 
均匀 分 布 unif 
Т1 КЗ 
正 态 分 布 norm 
2 тя 
сы 了 了 分布 f 
гл ват 
В 分 布 beta 
x 分 布 chi2 




















注 : 表 中 的 分 布 名 称 取 自 英文 单词 的 前 几 个 字母 。 
8.2 离散 随机 变量 概率 的 计算 
概率 P(X = k) TEFIE: 
>> пате pdf(k, Al, А2, АЗ) 
zk >> pdf('name', К, Al, A2, АЗ) 
Hp, name HK 8-1 的 MATLAB 名 称 ， 参 数 (k，Al,， 
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А2, АЗ) 8 128-2, 




















= 8-2 
分 布 名 语句 参数 设置 
二 项 分 布 binopdf (Е) 
泊 松 分 布 poisspdf (x, À) 
几何 分 布 geoodf (w, р) 
超 几 何 分 布 hygepdf (x, N, M, k) 








例 8.1: #£ X< B(3, 0.05), Ye P(0.1), ЖР(Х= 
2), Р(У=2). 
#1: MATLAB 程序 如 下 
>> Px =binopdf(2 ,3 ,0. 05) 
Px = 
0. 0071 
>> Py =poisspdf(2,0. 1) 
Py = 
0. 0045 
注 : 如 果 name 是 连续 随机 变量 的 MATLAB 名 称 ， 则 
name pdf(k，al，a2，a3) 所 计算 的 是 相应 的 密度 函数 f(x) 
ЖЕ x = 上 时 的 值 ， 是 函数 值 而 不 是 概率 值 。 如 : 
>> normpdf (0,0,1) 
ans = 
0. 3989 
>> 1/sqrt(2 * pi) 
Ans = 
0. 3989 


不 难看 出 ，normmpdf(0，0，1) 计算 的 是 标准 正 态 分 布 的 
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ЗЕ Е RZS) ЖЕ x =0 处 的 函数 值 


1 
JO) = 一 一 
т 
8.3 分 布 函 数 的 计算 一 一 连续 函数 概率 的 计算 
ЛАР Р(х) =P(X<x) 的 计算 程序 格式 如 下 : 
>> пате cdf(x,Al ,A2 ,A3 ) 
或 >>cdf('name',x,Al,A2,A3) 
其 中 ，name 见 表 8-1 中 MATLAB 名 称 ， 离 散 型 分 布 的 
参数 设置 见 表 8-3 ， 其 他 连续 型 分 布 的 参数 设置 见 表 8-3。 




































































表 8-3 
分 布 名 称 记 5 W “J 参数 设置 
均匀 分 布 Xc U(a, b) unifedf (x, a, b) 
指数 分 布 XvE(1/A) expcdf (x, À) 
正 态 分 布 Х‹оМ(ш, о?) normcdf (х, u, о) 
t 分 布 Xcoti(P) tcdf (t, n) 
РА XOF(V, V) fedf (x, vl, v2) 
Гг Э Хо (a, b) gamcdf (x, a, b) 
Вв 分 布 Xc B(a, b) betacdf (x, a, b) 
x 分 布 Xay (п) chi2cdf (x, n) 
例 8.2: ë Xo N( - 1, 4), RP(-5<X<1) MP 
(|X|>2)。 


解 : 由 XoN( -1, 4) 知 : 人 = -1, o=2 


(1) P( -5<X<1) =F(1) -F( -5) 得 程序 : 


>> РІ = normcdf(1, 


Ans= 


-1, 2) -normcdf( -5, 


г’ 
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0. 8186 
(2) P(|X|>2) =1-Р(|Х|=2) =1-Р(Х=2) +Р(Х 
< -2) 得 : 
>> P2 =1 -normcdf(2, -1,2) +normcdf( -2, -1,2) 
Ans = 
0. 3753 
8.4 期 望 和 方差 的 计算 


1. 常见 分 布 的 期 望 与 方差 
常见 分 布 的 期 望 和 方差 的 计算 程序 格式 为 : 
[М,У] =namestat(Al ,A2 ,A3 ) 

其 中 ，name 见 表 8-1 中 MATLAB 名 称 ， 泊 松 分 布 的 调 
HAA poisstat( 少 一 个 s) ，(A1，A2，A3 ) HK 8-2 和 表 
8-3 中 除 掉 x 后 的 部 分 。 

例 8.3: AXAN, 9), Yok(5), 求 X 和 YY 的 期 望 
和 方差 。 

解 : (1) X 的 期 望 和 方差 的 计算 程序 为 : 

>> [M1 ,V1 ] =normstat(1,3) 

МІ =1 

V1 =9 

















1 
(2) 由 YoE(5) 知 ,一 =5, 和 =0.2。Y 的 期 望 和 方差 


的 计算 程序 为 : 
>> [М2,У2 | = expstat( 0. 2) 
М2 =0.2 
У2 =0. 4 
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2. 一 般 的 数学 期 望 和 方差 
离散 随机 变量 的 数学 期 望 通常 表示 为 E(X) = 
ХУ хр, Жр, 是 对 应 x; 的 概率 。 在 MATLAB 中 用 函数 sum 


求 这 种 简单 数学 期 望 。 
sum(x) AIRRA 





+оо 


连续 随机 变量 Xe f(x) 的 数学 期 望 EX = ЕС Jdr, 


其 中 f(x) 是 随机 变量 X 的 概率 密度 函数 。 在 MATLAB 中 用 
函数 int(x* fx, a, b) 计算 这 种 数学 期 望 。 

int(x * fx,a,b) 积分 函数 

样本 均值 给 定 一 组 样本 x = (x,，x,，…，%,)， 样 本 均 


fx = 一 六 = 。 此 时 期 望 值 等 于 各 元 素 的 算术 平均 值 。 


在 MATLAB 中 用 函数 Mean 求 样本 均值 。 

Mean 平均 值 函数 

对 于 样本 x = (x ，x,，…，%,)，mean (x) 得 到 样本 均 
值 x， 对 于 向 量 x，mean(x) 得 到 它 的 元 素 的 平均 值 ， 对 于 
MEX, теав( ХО 得 到 一 列 向 量 ， 其 每 一 行 值 为 矩阵 行 元 素 
的 平均 值 。 

对 随机 变量 ,为 表征 它 与 其 期 望 的 偏离 程度 ,引入 方差 
和 标准 差 。 

方差 D(X) =E|[X-E(X) | =E(X) - (ЕХ)? 

标准 差 oo(X) = /D(X) 

对 一 组 样本 x = (x ,x,,… ,x ) ,样本 方差 为 
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样本 标准 差 为 
l < - 
e J У (x; — x)? 


п- 1% 
ТЕ MATLAB 中 由 方差 函数 var 和 标准 差 函 数 std 计算 , PK 
数 形式 : 


var(x,1) 为 方差 | 管 前 因子 为 二 ) ;var(x) 为 样本 方差 











(午前 因子 为 一 )。 


1 
sd (x, D 为 标准 差 ( РТА f ) зма Сх) 为 样本 标准 








ә (язе). 


8.5 逆 累 加 分 布 


逆 累 加 分 布 函数 用 于 对 给 定 的 概率 P, 求 *, 使 得 随机 事 
件 P(X<x) =p, 其 计算 程序 为 : 
>> пате inv(p,Al,A2,A3) 
其 中 ,name 见 表 8-1 中 MATLAB 名 称 ,参数 (AL,A2 ,A3 ) 见 表 
8-2 R 8-3 中 除 掉 x 后 的 部 分 。 
例 8.4: 若 XeB(5,0.1) ,YoP(5) ,P(X<a) =0.9914, 
Р(У>Ь) =0.01, 求 a、b。 
解 :(1) >>a=bino inv(0. 9914,5,0.1) 
a= 
2 
(2) HT P(Y<b)=1 -0.01 =0. 99 
>> b = poiss inv(0. 99,5) 
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b= 
11 

41 8. 5:1 Хо N(5.9) ,Р(Х<а) =0. 9525, a. 
解 :由 XoN(5,9) 知 ,w=5,o =3。 因 此 

>> а = norm inv(0. 9525 ,5 ,3 ) 
a= 
10. 0088 
J 8. 6:21 X< E(0.0001) ,P(X<b) =0. 633 Б. 
Е: ХФ E(0. 0001): Л = 10000 ,因此 

>> b = exp шу (0. 633 ,10000) 

b= 

1. 0024е +0. 04 


Z 训 w 


1. 计算 正 态 分 布 YI0,1) 的 随机 变量 工 在 点 0.6578 的 密度 函数 
值 。 




















2. 自由 度 为 8 的 卡 方 分 布 , 在 点 2. 18 处 的 密度 函数 值 。 
3. 求 标准 正 态 分 布 随机 变量 工 落 在 区 间 ( - = ,0.4) 内 的 概率 (该 
值 就 是 概率 统计 教材 中 的 附 表 一 一 标准 正 态 数值 表 ) 。 

4. 求 自由 度 为 16 的 卡 方 分 布 随机 变量 落 在 [0,6. 91] 内 的 概率 。 

5. %Х-М3,2?)ЖР|2<Х<5],Р|-4<х<10},Р!|Х| > 
2] ,PIX>3| 

6. 设 随机 变量 X 的 分 布 律 为 : 


X -2 -1 0 1 2 









































P 0.3 0.1 0.2 0.1 0.3 


求 E(X) EX -1)。 
7. 设 随 机 变量 XX 的 分 布 律 为 : 
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X =2 0 2 





Р. 0.4 0.3 0. 3 


ЖЕ(Х),Е(3Х +5), 
8 设 随 机 变量 X 的 分 布 密度 函数 为 : f (x) = 
3 6, 
АЕ я (09541) Ех, 
0, 其 他 
9. 规定 某 型 电子 元 件 的 使 用 寿命 超过 15008 为 一 等 品 ,已 知 一 批 
1, 一 等 品 率 为 0.2 , 问 这 批 样品 中 一 等 品 元 件 的 期 望 与 方差 为 
一 等 品 元 件 的 个 数 的 最 大 值 )? 
10. 求 参 数 入 =6 的 泊 松 分 布 的 期 望 与 方差 。 
11. 随机 抽取 6 个 滚珠 测 得 直径 (单位 :mm) 如 下 : 
14.70 15.21 14.90 14.91 15.32 15.32 
试 求 样本 平均 值 。 
12. 求 下 列 样 本 的 样本 方差 和 样本 标准 差 .方差 和 标准 差 ; 
14.70 15.21 14.90 15.32 15.32 
13. 随机 的 取 8 只 活塞 环 , 测 得 它们 的 直径 (单位 :mm ) 为 : 
74.001 74.005 74.003 74.001 74.000 73.998 74.006 
74. 002 


求 样本 的 均值 方差 值 样本 方差 值 标准 差 值 样本 标准 差 值 。 



































У 





样品 20 
多 少 (有 





多 少 
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9.1 样本 的 数字 特征 


在 MATLAB 中 ， 样 本 车 为 向 量 ， 则 其 数字 特征 是 对 





中 的 数字 作 相 应 的 运算 ; 





统计 问题 的 MATLAB 计算 





若 样本 为 矩阵 ， 则 其 数字 特征 是 对 





关中 的 数字 作 相 应 的 运算 。 样 本 常见 的 数字 特征 见 表 9-1。 


















































Ж 94 
名 ж 数学 式 描述 MATLAB 程序 格式 
和 Ум о 
最 大 值 max x, [m,n] =max(X) 
最 小 值 min x, [m,n] =min( X) 
中 位 数 M =median(X) 
算术 平均 =Ë D Hent 
方差 У и -9 var(X) 
标准 差 Vvar(X) std( X) 
协 方差 Se -5 (w = 7) Cov(X,Y) 
相关 系数 Eo СЫР Corrcoef( X , Y) 
VD(X) VD(Y) 
нунде D larez] mad(X) 





.70 . 
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(2) 
名 Ж 数学 式 描述 MATLAB 程序 格式 
k FGE Ly - x)" moment( X ,k ) 
П kal 

PR Х-Е(Х) ү 
За 度 skewness ( X 
ыы = D( X) ) Aa 

х X-E(X) ү 
ВЕКЕ iz kurtosis( X 
样本 峰 a | rtosis( X) 

















Е: 样本 中 缺失 的 值 可 用 пап 补充 ， 在 求 和 、 求 最 大 
(小 ) 值 、 求 算术 平均 值 、 求 中 值 和 求 标准 差 时 只 需 在 相应 
的 MATLAB 语句 前 加 上 nan 即 可 。 

例 9.1: 将 某 校 19 名 同学 随机 分 成 四 组 进行 体能 测试 ， 

















测 得 的 结果 见 表 9-2。 试 计算 各 组 的 和 、 最 大 值 、 中 值 、 均 




















值 和 标准 差 。 
表 9-2 
序号 组 1 组 2 组 3 组 4 

1 24 29 30 27 
2 26 25 28 31 
3 20 21 32 32 
4 28 27 30 33 
5 28 26 

6 30 











解 : МАТГАВ 程序 如 下 : 
>> X = [24 26 20 28 пап nan 
29 25 21 27 28 30 
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30 28 32 30 26 пап 

27 31 32 33 nan nan | 
>> nansum( X) 
Ans = 

98 160 146 123 % 计 算 各 组 数据 之 和 
>> [М,М] = паптах(Х) 


M=28 30 32 33 % 各 组 中 的 最 大 值 
N=4 6 3 4 % 各 组 中 的 最 大 值 所 在 的 序号 


>> М = nanmedian ( X) 
M=25 27.5 30 3l.5 % 各 组 的 中 值 
>> nanmean( X) 
Ans =24. 5000 26.6667 29.2000 30.7500 
% 各 组 的 算术 平均 值 
>> nanstd( X) 
Ans =3.4157 3.2660 2.2804 2.6300 % 各 组 的 标准 差 


9.2 参数 估计 


参数 估计 包括 点 佑 计 和 区 间 佑 计 , 常 用 方法 有 和 矩 法 和 极 
大 似 然 法 。 由 于 极 大 似 然 佑 计 的 佑 计量 不 仅 满足 无 偏 性 有 
效 性 等 基本 条 件 ,还 能 保证 其 为 充分 统计 量 ,所 以 在 点 估计 和 
区 间 估 计 中 ,一 般 推 荐 使 用 极 大 似 然 法 。 参 数 估计 的 程序 格 
式 为 : 

[参数 估计 值 ] = namefit(X, 参数 设置 ) 
其 中 , 下 为 向 量 或 矩阵 ， 为 矩阵 时 是 对 互 的 列 进 行 运算 。 

1. 常用 分 布 的 参数 估计 

(1) 正 态 分 布 XAN (4,o ), 但 参数 ос Ж 
知 ， 则 其 点 估计 和 区 间 估 计 的 MATLAB 程序 为 : 
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>> [ ши, sigma, muci, sigmaci| =normfit( X, Alpha) 
HP, 1) mu, sigma 为 人 和 的 点 估计 。 
2) muci, sigmaci Я ш Жо 的 100 (1 - Alpha) % 置信 
度 的 置信 区 间 估 计 ， 其 行 数 为 2， 列 数 与 下 的 列 数 相同 ， 
上 、 下 值 分 别 为 置信 区 间 的 上 、 下 限 。 
3) Alpha 为 给 定 的 显著 水 平 a， 表 示 置 信和 度 为 (1-а) 
100% 。 默 认 值 为 0. 05, ， 即 置信 度 为 95% 。 
例 9.2: 已 知 某 种 灯泡 的 寿命 服从 正 态 分 布 ， 在 某 周 内 
从 所 生产 的 该 种 灯泡 中 随机 抽取 10 只 ， 测 得 其 寿命 〈 单 
位 : h) X: 
1067, 919, 1196, 785, 1126, 936, 918, 1156, 920, 
948 
设 总 体 参 数 均 未 知 ， 求 该 周 中 生产 的 灯泡 寿命 超过 1300h 概 
率 的 极 大 似 然 估计 值 。 
解 : Xo N(u,o2), Мр(Х>1300) =1 -p(X=1300) = 
1300 -u 
. 


G 


Eu Mo, 分 别 是 1 和 的 极 大 似 然 估计 ， 则 由 极 天 似 然 估 
1300 — u, 
计 的 不 变性 ， 所 求 极 大 似 然 估计 值 为 : p =1 - | | 


计算 程序 如 下 : 
>>X = [ 1067 ,919 ,1196 ,785 ,1126 ,936 ,918 ,1156 ,920 ,948 | 
>> [mu,sigma,muci,sigmaci] = normfit ( X ) 
mu = 997. 1000 sigma = 131. 5476 
muci = 1. Qe + 003 * 
0. 9030 
1. 0912 




















С 
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Sigmaci =90. 4831 






































240. 1546 
>>р=1 -normcdf( (1300 — mu)/sigma) 
Р = 
0.0107 
(2) 其 他 分 布 〈 表 9-3 ) 
表 9-3 
类 H 概率 分 布 MATLAB 程序 格式 
二 项 分 布 [ phat,pci ] = binofit( X , Alpha) 
mou 
泊 松 分 布 [lambda ,lambdaci | = poissfit( X , Alpha) 
均匀 分 布 [a,b,aci,bci] =unifit( X , Alpha) 
ВЖЕ [ lambda ,lambdaci | = expfit( X , Alpha) 
连续 型 
Г 2518 [phat,pci] = gamfit( X , Alpha ) 
В 分 布 [ phat,pci ] = betafit( X , Alpha ) 














例 9.3: 为 了 估计 制造 某 种 产品 所 需要 的 单 件 平 均 工时 
(单位 : p) ， 现 制造 10 件 ， 记 录 每 件 所 需 的 工时 如 下 : 

9.8, 10.2, 10.6, 11.2, 11.4, 12.5, 12.8, 9.9, 
10.4, 9.9 
设 制造 单 件 产 品 所 需 工 时 服从 指数 分 布 ， 求 平均 工时 的 极 大 
似 然 估计 值 和 95% 的 置信 区 间 。 

解 : >> X = (9.8, 10.2, 10.6, 11.2, 11.4, 12.5, 
12.8, 9.9, 10.4, 9.9) 

>> [mu, muci] =expfit( X , 0. 05) 

mu = 10. 9100 

muci = 5. 2318 

18. 6395 
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即 平均 工时 的 极 大 似 然 估 计 值 为 10. 91h，95% 置信 度 的 置 
2318, 18. 6395 ] 。 
参数 估计 的 通用 命 RD 令 

АР Gomphus 
有 一 个 通用 程序 : 
>> | phat,pci] =mle('name',X,Alpha) 
其 中 name 见 表 8-1 中 MATLAB 名 称 。 

例 9.4: 从 一 批 产 品 中 抽取 100 个 ， 经 检验 有 60 个 一 
级 品 ， 试 在 = 0. 05 的 显著 水 平 下 确定 这 批 产品 的 一 级 品 


о 











м 
я 


#: >>п =100, X =60, Alpha =0. 05 

>> | рва ‚ре! ] =mle('bino',X,Alpha,n) 

phat = 0. 6000 

pci =0. 4972 

0. 6967 

Е а =0. 05 的 显著 水 平 下 这 批 产品 的 一 级 品 率 为 0.6， 置 
信 度 为 95% 的 置信 区 间 为 [0. 4972, 0. 6967 ] 。 

Е: mle( ) 语 句 中 的 n 是 相应 二 项 分 布 的 额外 参数 。 


9.3 假设 检验 


假设 检验 是 对 总 体 X 的 分 布 或 分 布 参数 作 某 种 假设 ， 
然后 根据 抽样 的 样本 观察 值 ， 运 用 数理 统计 的 分 析 方 法 ， 检 
验 这 种 假设 是 否 正确 ， 从 而 决定 接受 假设 或 拒绝 假设 。 假 设 




















检验 分 以 下 两 大 类 : 
QD 参数 检验 : 观测 的 分 布 函数 类 型 已 知 ， 对 总 体 的 参数 
及 有 关 性 质 作 出 明确 的 判断 ; 


@) 非 参数 检验 : 要 求 判断 总 体 分 布 类 型 的 检验 。 
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例 9.5: 某 车 间 用 一 台 包 装机 包装 糖 。 当 机 器 正常 运转 
时 ， 每 袋 糖 的 质 ( 重 ) 量 的 均值 为 0. 5kg， 标 准 差 为 0.015， 
某 日 开工 后 检验 包装 机 是 否 正 常 ， 随 机 抽取 9 袋 ， 称 得 净重 
(单位 : kg) 为 : 

0.497 0.506 0.518 0.524 0.498 0.511 0.52 
0.515 0.512 问 机 器 运转 是 否 正 党 ? 

ВЕ: 假设 袋 装 糖 的 质量 是 一 个 随机 变量 ， 它 服从 正 态 分 
布 。 

已 知 总 体 均值 人 w =0.5， 标 准 差 o =0. 015, ， 该 问题 是 当 
o° 为 已 知 时 ， 在 显著 性 水 平 a =0.05 下 ， 根 据 样 本 值 判 断 
均值 人 =0.5 还 是 凡 关 0. 5。 

由 于 标准 差 已 知 ， 验 证 总 体 的 均值 ， 采 用 z 检验 。 

原 假设 : 有 :pn =o =0.5; 备 择 假设 : А: 70.5. 

用 MATLAB 求解 如 下 : 
>>X = [0.497 0.506 0.518 0.524 0.498 0.511 
0.52 0.515 0.512]; 
>> h,sig,ci,zval] = ztest( X ,0. 5 ‚0. 015 ‚0. 05) 

Н= 

1 























Sig = 

0. 248 % 样本 观察 值 的 概率 
Ci = 

0. 5014 0.5210 % 置 信 区 间 , 均值 0.5 在 此 区 间 之 外 
Zval = 

2. 2444 % 统计 量 的 值 


结果 表明 : h=1, 说 明 在 显著 性 水 平 a=0.05 下 ， 可 拒 
绝 原 假设 ， 即 认为 包装 机 工作 不 正常 。 
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Е: 在 已 知 总 体 服从 正 态 分 布 的 条 件 下 ， 若 总 体 方差 
о? 已 知 ， 则 总 体 均 值 的 检验 使 用 z 检验 ，z 检验 的 命令 为 : 

[h,sig,ci] = ztest(x,m ,sigma alpha ,tail) 

其 中 ，m 为 均值 ，sigma 为 已 知 方差 ，alpha 为 显著 性 水 平 ， 
alpha 的 默认 项 为 0.05，tail 的 取 值 决定 检验 内 容 。 

Tail =0， 检 验 假设 “X 的 均值 =m” 

Тай =1， 检 验 假设 “X 的 均值 >m” 

Тай = -1， 检 验 假设 “X 的 均值 <m” 

Tail 的 默认 值 为 0，h = 1 表示 拒绝 原 假设 ，h =0 表示 不 
拒绝 原 假设 ，sig 为 假设 成 立 的 概率 ，ci 为 均值 的 1 - alpha 
置信 区 间 。 

例 9.6: 某 公 司 为 测 某 电子 元 件 的 寿命 ， 进 行 抽样 检 
查 ， 现 测 得 16 只 元 件 的 寿命 (单位 : h) 为 : 

159 280 101 212 224 379 197 264 222 362 
168 250 149 260 485 170 
是 否 有 理由 认为 元 件 的 平均 寿命 大 于 jo (w =225h)? 

解 : 假设 电子 元 件 的 寿命 X (单位 : h) 服从 正 态 分 布 。 

假设 所 得 数据 为 抽样 数据 ， 这 里 总 体 方差 o 未知， 故 
可 以 采用 t+ 检验 。 在 显著 性 水 平 =0. 05 下 检验 假设 。 

原 假 设 : Нм <p =225 

备 择 假设 : Ни: u >225 

用 MATLAB 求解 如 下 : 
>>X=|159 280 101 212 224 379 197 264 222 
362 168 250 149 260 485 170]; 
>> [h,sig,ci] = певі(х,225,0. 05,1) % JEX 225, W 

性 水 平 alpha = 
0. 05 
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0 
516 = 
0. 2570 
Сі = 
198.2321 Inf ”% 均 值 225 在 该 置信 区 间 内 

结果 表明 : В =0 表示 在 水 平 a = 0.05 下 ， 应 该 接受 原 
假设 有 ， 即 认为 元 件 的 平均 寿命 不 大 于 225h。 

注 : 在 已 知 总 体 服从 正 态 分 布 的 条 件 下 ， 若 总 体 方差 
о? 未 知 ， 则 总 体 均 值 的 检验 可 使 用 t 检验 , t 检验 的 命令 
为 : 

[h,sig,ci] = ttest(x,m,alpha,tail) 

其 中 ，m，alpha，tail 的 意义 同 z 检验 命令 中 相应 的 参数 的 
意义 。h =1 表示 拒绝 原 假设 ，h = 0 表示 不 拒绝 原 假设 ，sig 
为 假设 成 立 的 概率 ，ci 为 均值 的 1 -alpha 置信 区 间 。 

例 9.7: 在 平 炉 上 进行 一 项 试验 以 确定 新 操作 方法 是 否 
会 增加 钢 的 产量 。 试 验 在 已 知 平 护 上 进行 ， 每 炼 一 炉 钢 除 操 
作 方 法 外 ， 其 他 的 条 件 做 到 尽 可 能 相同 。 先 用 标准 方法 炼 一 
炉 ， 然 后 用 其 他 方法 炼 一 炉 ， 以 后 交替 进行 ， 各 炼 10 炉 其 
兰 量 分 别 为 : 

(1) 标准 方法 : 78.1 72.4 76.2 74.3 77.4 78.4 
76.0 75.5 76.7 77.3 

(2) 新 方法 : 79.1 81.0 77.3 79.1 80.0 79.1 79.1 
77.3 80.2 82.1 

新 操作 方法 能 和 否 提高 产量 ? (На =0. 05) 

解 : 假设 这 两 个 样本 分 别 来 自 正 态 总 体 WCmw o°) 和 
№, а?) , pi, ш, о? 未 知 ， 且 假设 两 个 样本 均 相 互 独立 。 
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两 个 总 体 方差 不 变 时 ， 因 涉及 新 旧 两 种 方法 的 比较 ， 又 
由 于 两 种 方法 按 相 近 的 原则 可 配 成 对 ， 以 消除 混杂 因素 的 影 
响 。 下 面 采用 配对 的 t 检 验 。 在 显著 性 水 平 a =0. 05 下 检验 
假设 。 

原 假 设 : Нуш = о; 

备 择 假设 : Hiu <р о 

用 MATLAB 求解 如 下 : 

>>X=[78.1 72.4 76.2 74.3 77.4 78.4 76.0 
75.5 76.7 77.3]; 

>> Ү= [79.1 81.0 77.3 79.1 80.0 79.1 79.1 
77.3 80.2 82.1]; 

>> [В, 9, с! | = Нез (Х ‚У ,0. 05, -1) 





h= 

1 
515 = 

2. 1759e -004 % 说 明 两 个 总 体 均 值 相等 的 概率 很 小 
Ci = 


-Inf -1.9083 
结果 表明 : h =1 表示 在 显著 性 水 平 a =0.05 F, MIE 
绝 原 假设 ， 即 认为 新 操作 法 提高 了 产 率 。 

Е: 两 总 体 均 值 的 假设 检验 使 用 + 检验 的 命令 为 [h sig, 
ci] =ttest(x,y,alpha,tal) ， 用 于 检验 数据 x，y 关于 均值 的 
某 一 假设 是 否 成 立 。 其 中 ，alpha 为 显 车 性 水 平 ，tail 意义 同 
前 。 

例 9.8: 在 自动 化 车 床 连续 加 工 某 零件 的 一 道 工序 中 ， 
由 于 刀具 损坏 等 原因 会 导致 机 器 出 现 故 障 ， 而 故障 的 出 现 是 
随机 的 。 现 有 100 次 故障 记录 ， 故 障 出 现时 该 刀具 完成 的 零 
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件数 〈 件 ) 如 下 : 

459 362 624 542 509 584 
612 452 434 982 640 742 
926 653 164 487 734 608 
527 552 513 781 474 388 
775 859 755 49 697 515 
402 960 885 610 292 837 
699 634 555 570 84 416 
447 654 564 339 280 246 
621 724 531 512 577 496 
764 558 378 765 666 763 


433 
565 
428 
824 
628 
473 
606 
687 
468 
217 


748 
706 


815 505 
593 680 


1153 593 844 


538 
954 
677 


862 659 
771 609 
358 638 


1062 484 120 


539 
499 
715 


790 581 
544 645 
310 851 


试 分 析 该 刀具 出 现 故 障 时 完成 的 零件 数 服从 哪 种 分 布 。 
解 : 假设 工作 人 员 常 通过 检查 零件 来 确定 工序 是 否 出 现 
им; ВЕН 8 服从 正 态 分 布 。 
为 分 析 该 刀具 出 现 故 障 时 完成 的 零件 数 服 从 哪 种 分 布 ， 
先 作 频率 分 布 直方 图 ， 若 观察 出 样本 服从 正 态 分 布 ， 则 用 





MATLAB 中 的 函数 normplot( ) IH H 








一 条 直线 上 ， 则 表明 样本 服从 正 态 














态 分布 了 。 先 作 频 率 分 布 直方 图 。 




















出 样本 。 如 果 样 本 都 分 布 在 
人 分布， 否则 不 服从 正 态 分 
布 。 接 着 用 normfit( ) 进行 分 布 的 正 态 性 检验 ， 最 后 用 函数 
ztest( ) 进行 参数 检验 。 这 样 就 可 以 确定 一 组 数据 是 否 服从 正 


>>х= [459 362 624 542 509 584 
612 452 434 982 640 742 
926 653 164 487 734 608 
527 552 513 781 474 388 
775 859 755 49 697 515 
402 960 885 610 292 837 


433 
565 
428 
824 
628 
473 








748 815 505 
706 593 680 
1153 593 844 
538 862 659 
954 771 609 
677 358 638 
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699 634 555 570 84 416 606 1062484 120 
447 654 564 339 280 246 687 539 790 581 
621 724 531 512 577 496 468 499 544 645 
764 558 378 765 666 763 217 715 310 851]; 
>>hist(x(:),10) % 频 率 分 布 直方 图 
得 到 图 9-1 所 示 的 频率 分 布 直 方 图 。 








25 





20r 











800 1000 1200 


从 图 中 可 以 初步 认定 ， 该 刀具 出 现 故障 时 完成 的 零件 数 
服从 正 态 分 布 ， 但 必须 进行 参数 估计 和 假设 检验 。 
1 (x -1)2 
正 态 分 布 的 概率 密度 为 : f(x) = e 22, Hp 
分 布 的 概率 密度 为 : f(x) т ш 


TO 








是 平均 数 ，o 是 标准 差 。 
(1) 分 布 的 正 态 性 检验 MATLAB 的 求解 如 下 : 
>>normplot(x(:)) % 显 示 x 中 数据 的 一 个 正 态 分 布 图 ， 


若 x 中 数据 基本 分 布 在 一 条 直线 上 ， 则 x 服从 正 态 分 布 。 由 
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图 9-2 可 知 ， 刀 具 出 现 故障 时 完成 的 零件 数 近似 服从 正 态 分 
布 。 





Probability 














п H H Н i р i i 
0.003L-—4 css ss 800 ЗЕЕ pa p= =-= === == === Ea 
100 200 300 400 500 600 700 800 900 1000 1100 
Data 
图 9-2 


(2) 参数 估计 

>> [muhat ,sigmahat ,muci,sigmaci] = normfit( x( : ) ) 
% 正 态 分布 的 参数 估计 函数 
muhat = 
594 
sigmahat = 
204. 1301 

muci = 

553. 4963 

634. 5038 
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sigmaci = 
179. 2276 
237. 1329 
运行 结果 估计 出 该 刀具 完成 零件 数量 的 均值 人 = 594, 
7 25 о? = 204, 均值 的 0.95 置信 区 间 为 [553.4962 ， 
634.5038], МЭИ 0.95 的 置信 区 间 为 179.2276， 
237. 1329 ] 。 
(3) 假设 检验 
已 知 刀具 的 寿命 服从 正 态 分 布 ， 现 在 在 方差 未 知 的 情况 
下 ,检验 其 均值 m 是 否 等 于 594， 使 用 t 检验。 
>> [h,sig,ci] = ztest(x( :) ,594,0. 05) 
h = 
0 
Slg = 
1 
cl = 
[ 553. 4962 ,634. 5038 | 
结果 表示 : h=0 表示 不 拒绝 原 假设 , 说明 提出 的 假设 
寿命 均值 为 594 是 合理 的 。 
综 上 所 述 ， 可 以 认为 刀具 出 现 故 障 时 完成 的 零件 数 服 从 
正 态 分 布 ， 刀 具 的 平均 寿命 为 594。 
注 : 设 总 体 服从 正 态 分 布 ， 则 其 点 估计 可 同时 用 以 下 命 























[ muhat , sigmahat , muci , sigmaci | = normfit( X ,alpha) 
表示 在 显著 性 水 平 alpha 下 估计 数据 X 的 参数 ， 返 回 值 mu- 
hat 是 X 的 均值 的 点 估计 值 ，sigmahat 是 标准 差 的 点 估计 值 ， 
muci 是 均值 的 区 间 估 计 ，sigmaci 是 标准 差 的 区 间 估 计 ， 
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alpha 默 认为 0.05. 
9.4 回归 分 析 


回归 分 析 是 在 观察 数据 的 基础 上 ， 利 用 数学 方法 确定 变 
量 间 函数 关系 的 一 种 统计 分 析 方 法 。 在 回归 分 析 中 ， 如 果 自 
变量 和 因 变 量 的 关系 可 用 Y=b。 +bX 表示， 那么 这 种 回归 称 
为 线性 回归 ， 否 则 为 非 线 性 回归 ; 耕 自 变量 的 个 数 为 一 个 ， 
则 称 为 一 元 回归 ; 竺 自 变量 的 个 数 为 两 个 或 两 个 以 上 ， 则 称 
为 多 元 回归 。 

(1) 线性 回归 ”线性 回归 的 一 般 数学 模型 为 : Y = b, + 
bX +: +b X... 

我 们 希望 m 组 测量 数据 (Y(j) , G)... ,Х,() (< 

















ism) MENTE 
MATLAB 的 计算 程序 如 下 : 
>>X1 = [ ее 5 е ‚Хп = [ Е P N a| а 1 
Yo 按 列 的 方式 输入 各 变量 的 数据 
>> Х = [опеѕ(т,1) ХІ, + ,Xn | 
% m 为 各 变量 的 行 数 
>> [b,bint,r,rint,stats] =regress(Y,X,a) % 作 线性 回归 
>> rcoplot (r,rint) % 作 残 差 图 








HF, b 是 回归 系数 的 最 小 二 乘 佑 计 ，bint 是 置信 和 度 为 100 
(1-а) % 的 各 回归 系数 的 置信 区 间 ，r 为 观测 值 Y 与 估计 值 
的 残 差 向 量 ，rint 是 置信 度 为 100(1 — a.) % 的 各 残 量 的 区 间 
估计 ，stats 返回 的 是 R 统计 量 (相关 系数 的 平方 值 ) Е 2 
计量 和 显著 性 概率 P 值 。 

注意 : 当 所 有 数据 的 残 差 图 都 过 0 线 ， 且 显著 性 概率 P 
<0.01 时 ， 回 归 效 果 显 著 。 奉 某 个 数据 的 残 差 图 不 过 0 线 
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(此 时 数据 对 应 的 残 差 线 为 红线 ) ， 则 视 为 异常 值 ， 将 其 吻 
除 ， 再 重新 进行 回归 。 

例 9.9: 某 寝室 5 人 的 身高 和 体重 数据 如 下 : 
身高 (сш): 175, 165, 185, 178, 195 
体重 (kg): 73, 68, 81, 75, 90 
试 确定 这 5 人 的 身高 和 体重 的 关系 。 

解 : 

>>XL =[175, 165, 185, 178, 195 ]' 

>> Ү= [73, 68, 81, 75, 90 ]' 

>> = [опе (5, 1), ХІ] 

>> [Ь, bint, г, rint, stats] =regress( Y, X, 0. 05) 
>> гсорю (г, rint) 


其 残 差 图 如 图 9-3 所 示 。 














Residuals 





0.5 1 15 2 25 3 35 4 45 5 5.5 
Case Number 


图 9-3 
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wam 


注意 : 第 二 个 数据 对 应 的 残 差 图 不 过 0 26, ВЯ 
二 个 数据 后 继续 作 回归 : 
>> ХІ = [ 175, 185, 178, 195 | 
>> Ү= [73, 81, 75, 90]' 
>> Х = [опе (4, 1), ХІ] 
>> [Ь, bint, г, rint, stats] =regress( Y, X, 0. 05) 
Ь = -77. 5702 
0. 8585 
stats = 1. Ое + 003 * 
0. 0010 1. 3434 0. 0000 
>> rcoplot(r, гіпі) 
此 时 所 有 的 残 差 图 均 过 0 2 〈 见 图 94) ， 且 显著 性 概率 几 
乎 为 0， 因 此 ， 回 归 效 果 显著 。 





ut 














Residual Case Order Plot 


Residuals 
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回归 方程 为 : Y= -77.5702 +0. 8585X, 
其 中 , х 为 身高 ， 了 为 体重 。 

作 线性 回归 前 ， 也 可 通过 散 点 图 看 其 线性 效果 ， 如 果 这 
些 散 点 基本 在 一 条 直线 上 ， 也 可 以 作 线 性 回归 。 

$1 9. 10: 公司 要 在 13 个 地 区 推销 一 种 新 产品 ， 预 计 的 
使 用 人 数 、 销 量 以 及 当地 的 人 均 收 入 见 表 9-4， 若 其 中 某 地 
使 用 人 数 增 加 到 500 万 ， 人 均 收 入 增加 到 5000 元 ， 问 应 如 
何 安排 供应 量 ? 







































































表 9-4 
Лил | 收入 /元 销量 / 千 箱 人 数 /万 | 收入 /元 销量 / 千 箱 
274 2450 162 330 2450 192 
180 3254 120 53 2560 55 
375 3802 223 430 4020 252 
205 2838 131 236 2660 144 
86 2347 67 157 2088 103 
265 3782 169 370 2605 212 
98 3008 81 
解 : 


>> ХІ = [274, 180, 375, 205, 86, 265, 98, 330, 53, 
430, 236, 157, 370 ]' 

>> X2 = [ 2450, 3254, 3802, 2838, 2347, 3782, 3008, 
2450, 2560, 4020, 2660, 2088, 2605 | 

>> Y = [162, 120, 223, 131, 67, 169, 81, 192, 55, 
252, 144, 103, 212] 

>>Х = [ опез(13,1),Х1,Х2 ] 
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>>[b, bint, r, rint, stats] =regress( Y, X, 0.05) 
b =8. 5905 
0. 4973 
0. 0073 
Stats = 1. де + 004 * 
0. 0001 1. 3337 0 
>> rcoplot(r, rint) 


日 残 差 图 〈 见 图 9-5) 可 以 看 出 ， 所 有 数据 的 残 差 图 都 过 0 



































线 ， 且 显著 性 概率 P 近似 为 0， 故 回归 效果 显著 。 


Residual Case Order Plot 





Residuals 





图 9-5 


因此 有 : Y=8. 5905 +0. 4973Х, +0. 0073Х, 
`4 X, =500, X, =5000, Y=294 (Е). 

说 明 : 因 前 面 没 有 介绍 分 布 拟 合 检验 ， 所 以 这 里 省 去 了 
残 差 服从 均值 为 0 的 正 态 分 布 检验 。 
>>[Н,Р, СУ] = Lillietest(r, 0. 05) 
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H=0 
>> шеап (г) 
Ans = 

– 1. 3118е -014 

(2) 多 项 式 回 归 如果 回归 函数 是 多 项 式 ， 则 称 之 为 
ZME. H MATLAB 程序 格式 为 : 
>>[Р, S] =polyfit(X, Y, п) 
其 中 ，X 为 自 变 量 观测 值 ，Y 为 因 变量 观测 值 ，n 为 回归 多 
项 式 的 次 数 ; P 为 降 需 多 项 式 的 系数 向 量 ; S 是 一 个 矩阵 ， 
可 用 于 polyval( ) 或 polyconf( ) 作 误差 估计 。 





(Dpolyval 进行 多 项 式 评价 
>> Ү = polyval(P, X) % 在 给 定 值 X 处 求 了 
的 预测 值 
>> [Ү, Delta] = polyval(P, X, S) % 得 到 误差 估计 站 + 
Delta 
(Opolycon{ 一 一 进行 多 项 式 评价 和 置信 区 间 估 计 








>>[Y, Delta] = polyconf(P, X, $) % 831 Y 的 95% 置信 
区 间 Y + Delta 
>>[Y, Delta] =polyval(P, X, $, а) 019521] Ү 的 100(1 - 
а) % 置信 区 间 
例 9.11: 表 9-5 所 列 数据 是 某 公 司 销售 数字 信号 处 理 器 
的 收入 (单位 : 千 万 元 ) 。 
表 9-5 
年 份 | 2001 2002 | 2003 | 2004 | 2005 | 2006 | 2007 | 2008 | 2009 




















收入 | 3.2 29 2.8 3 3: 3 3.6 4 4.6 9:2 
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处 理 器 的 收入 是 多 少 ? 

解 : 先 做 散 点 图 以 确定 回归 函数 的 类 型 : 
>>1=[0,1,2,3,4,5,6,7,8] 
>>у=[3.2, 2.9, 2.8, 3, 3.3, 3.6, 4, 4.6, 5.2] 
>> ро (Е, у, '*') 
从 而 得 到 如 图 9-6 所 示 的 散 点 图 ， 其 图 像 像 一 条 二 次 抛物 
线 。 

















因此 考虑 作 二 次 多 项 式 回归 : 
>>[Р, $] =polyfit(t, у, 2) 
>>Р= 
0. 0596 – 0. 2087 3. 1055 

$ = 

В: [3 x3 double ] 

Df. 6 

Normr: 0. 2004 
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>> Ү =polyval (P, 9) 
Y= 
6. 0571 
故 2010 年 该 公司 销售 数字 信和 号 处 理 器 的 收入 是 6. 06 FH 
љо 
(3) 非 线 性 回归 在 许多 问题 中 ， 回 归 参 数 都 不 是 线 
性 的 ， 也 不 能 通过 转换 的 方式 将 其 变 为 线性 参数 ， 这 类 模型 
称 为 非 线性 回归 模型 。 其 MATLAB 程序 格式 为 : 
>> [beta, г, J] =nlinfit(X У, 'fun', beta0 ) 
其 中 ，X 为 自 变 量 失 了 泗 ， 每 一 列 对 应 一 个 变量 ，Y 为 因 变 量 
观测 向 量 ， 它 与 X 的 行 数 相同 ; fun 为 一 自 定 义 的 M 函数 文 
件 ，beta0 为 beta 的 迭代 初 值 ; r+ 是 观测 值 与 拟 合 值 的 残 量 ; 
ТАИТ, EZF nlparci, nlpredci 的 输入 ， 
生成 模型 参数 和 模型 预测 值 的 置信 区 间 。 
QO 输出 拟 合 曲线 图 
>>nlintool(x, у, 'fun', beta0, alpha, 'xname ', ' yname ') 
@ 计 处 预测 值 的 置信 区 间 
>> [ypred, delta] =nlpredci ('fun', х0, beta, г, J, alpha) 
其 中 ，ypred 是 x0 处 的 预测 点 估计 ，delta 是 预测 半径 ， 置 
ТЕЛЕ М 100 (1 — alpha) % 的 预测 区 间 为 
[ypred - delta, ypred + delta ] 
人 计算 参数 估计 值 的 置信 区 间 
>> ci = nlparci( beta, т, J) 
例 9.12: 体重 70kg 的 某 人 在 短 时 间 内 喝 下 2 瓶 啤酒 后 ， 
隅 一 段 时 间 测 他 血液 中 的 酒精 含量 (mg/100mL)， 结 果 见 
表 9-6。( 数 据 来 源 : 2004 年 全 国 大 学 生 数 学 建 模 竞 赛 题 ) 
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表 9-6 (单位 : mg/100mL) 
间隔 时 间 /h 10. 2510. 5 0.75| 1 11.5] 2 [2.5 3 |3.5 4 44.5 5 
酒精 含量 30 | 68 | 75 | 82 | 82 | 77 | 68 | 68 | 58 51 50 | 41 


























间隔 时 间 /h 6 7 8 9 | 10 | 11 | 12 |13 |14 |15 16 
酒精 含量 38 35 | 28 |25 | 18 | 15 | 12 | 10 | 7 7 4 

































































假设 该 人 现在 在 很 短 的 时 间 内 喝 了 3 瓶 啤酒 ， 通 过 微分 

方程 模型 得 知 他 血液 中 的 酒精 含量 ， 
C(t) =164. 877 (e° — et) 

已 知 国家 最 新 的 酒 后 驾车 标准 是 血液 酒精 含量 大 于 或 等 于 
20mg/100mL， 问 该 人 酒 后 需 经 过 多 长 的 时 间 才 可 以 开车 。 

ВЕ: (1) 编制 M XPF fun. m, 
在 MATLAB 界面 依次 单 击 Pile 一 New 一 M -fle， 然 后 在 其 窗 
口 输入 下 面 语句 : 
Function y = fun (beta0, х) 
K =beta0 (1) 
у = 164. 877 * (exp( – 0. 183 жх) -exp( -k *x)) 
输 完 后 以 默认 的 文件 名 fun 保存 到 MATLAB 的 Work 文件 夹 
(不 能 保存 到 其 他 文件 夹 ) 。 

(2) 打开 fun. m ( 需 单 击 Window МЕ, ), 在 
MATLAB 命令 窗口 运行 : 
>>1=3,1= [ 0.25, 0.5, 0.75, 1] 
>>у= [30, 68, 75, 82] * n/2 
>> beta0 = [0.2] 
>> beta, г, J] =nlinfit(t, у, 'fun', beta0) 
k= 

2. 1360 
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从 而 有 C(t) =164. 877(е7°'®% -e7 8%) 

在 Excel 空白 表格 的 Al 单元 格 输入 0.01, А2 单元 格 输 
和 人 0.02， 选 定 Al 和 A2 单元 格 , 将 鼠标 放 至 右 下 角 变 黑 
“十 ” 字 后 ， 按 住 左 键 下 拉 至 A20; 在 ВІ 单元 格 中 输入 
C(t) 的 表达 式 后 按 回 车 键 ， 再 将 鼠标 放 至 ВІ 的 右 下 角 变 黑 
“十 ” 字 ， 按 住 左 键 下 拉 。 看 到 1 二 9. 32h 后 , C(t) <20。 

因此 ,该 人 酒 后 需 经 过 9. 32h 后 才 可 以 开车 。 

хим 

1. 假设 某 种 清漆 的 9 个 样品 ， 其 干燥 时 间 (以 小 时 计 ) 分 别 为 
6.0, 5.7, 5.8, 6.5, 7.0, 6.3, 5.6, 6.1, 5.0, 设 干燥 时 间 总 体 服 
MESSA №(ш,о?) 。 

(1) R u A o 的 点 估计 值 ; 

(2) 求人 和 ez 的 置信 度 为 0. 95 的 置信 区 间 (o RA) o 

2. 分 别 使 用 金 球 和 铂 球 测定 引力 常数 : 

(1) 用 金 球 测定 观察 值 为 : 6.683, 6.681, 6.676, 6.678, 
6.679, 6.672; 

(2) 用 铂 球 测 定 观 察 值 为 : 6.661, 6.661, 6.667, 6.667, 
6. 664。 

设 测定 值 总 体 为 Wu,o) ,人 和 为 未 知 ， 对 (1), (2) 两 种 
情况 分 别 求 % 和 o 的 置信 度 为 0.9 的 置信 区 间 。 

3. 有 两 组 (每 组 100 个 元 素 ) 正 态 随机 数据 ， 其 均值 为 10， 均 
52202, 0 95% 的 置信 区 间 和 参数 估计 值 。 

4. 某 种 钠 头 在 正常 情况 下 ， 按 规格 平均 净重 379g， 标 准 差 11g， 
现 抽查 10 个 测 得 数据 (单位 : g): 370.74 372. 80 386. 43 398. 14 
369.21 381. 67 367. 90 371. 93 386. 22 393. 08。 试 说 明 平 均 净重 是 否 符 
合 规 定 要 求 (а =0. 05). 

5. 某 制 药 厂 生产 一 种 抗生素 , 已 知 在 正常 情况 下 , 每 瓶 抗生素 的 
某 项 指标 服从 均值 为 22. 3 的 正 态 分 布 。 某 天 开工 后 , 测 得 10 ЖМЖ 


ПЫ = 
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据 为 : 22.3 21.5 21.7 23.4 21.8 21. 4 23.4 19.8 24. 4 21.2, ЖЯ 





抗生素 每 瓶 的 该 项 指标 的 均值 是 否 正常 ? 

6. 美国 某 加 油 站 有 2008 年 1 月 和 2 月 二 组 汽油 价格 数据 样 
119 117 115 116 112 121 115 122 116 118 109 112 119 112 117 113 
109 109 118 和 118 115 115 122 118 121 120 122 120 113 120 123 



























































E 产 


本 : 
114 
121 


109 117 117 120 116 118 125 ， 设 该 加 油 站 汽油 价格 的 标准 差 为 每 加 仑 























4 美 分 。 

(1) 判断 1 月 份 每 加 仑 汽油 平均 价格 是 否 为 1. 15 美元 ; (2) 设 
价格 标准 差 未 知 ， 比 较 1 月 、2 月 的 两 个 样本 均值 。 

7. 从 一 批 化 工 产品 中 取 8 份 样品 ， 每 份 一 分 为 二 ， 分 别 用 甲 、 乙 




















两 种 方法 检测 某 成 分 得 数据 (单位 : mg): 57 56 61 60 47 49 63 61 和 
65 69 54 60 52 62 57 60。 设 含量 测 值 N(u,o2) ， 问 两 法 测定 的 平均 值 














ESAERA? (а =0. 05) 
8. 近期 在 某 城市 一 中 学 收集 到 24 名 17 岁 男性 中 学 生 身 高 
位 : ст) 数据 如 下 : 
170.1 179.0 171.5 173.1 174.1 177.2 170.3 176.2 
163.7 175.4 163.3 179.0 176.5 178.4 165.1 179.4 
176.3 179.0 173.9 173.7 173.2 172.3 169.3 172.8 


























( 单 


又 查 到 20 年 前 同一 所 学 校 同龄 男生 的 平均 身高 为 168cem， 根 据 上 








面 的 数据 回答 ，20 年 来 城市 17 岁 男性 中 学 生 的 身高 是 否 发 生 了 变 




















响 ， 测 得 数据 如 下 : 
温度 X 100 10 120 130 140 150 160 170 180 
生成 率 Y 45 51 54 61 66 70 74 78 85 
试 做 直线 y=a + bx 型 回归 。 
10. 用 ployfit KAA y =a; +a,x +a 型 回归 。 
x 0.5 1.0 1.5 2.0 2:5 3.0 


化 ? 


9. 为 了 研究 某 一 化 学 反应 过 程 中 ， 温 度 工 对 产品 生成 率 了 的 影 


190 
89 





y 1.75 2.45 3.81 4. 80 7. 00 8. 60 
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10.1 线性 最 小 二 乘 拟 合 


多 项 式 拟 合 的 MATLAB 命令 实现 : 

ФЗ f(x) = ах" +ах" +… +a x tan WE, М 
定 多 项 式 系数 的 命令 : 
P =polyfit(x,y,m) 
НР, х= (x,t) Y = (yi, ) p = (4,9,7, 
an Е (к) =ax"”+ax" | + +a, +а, ША 
数 ; т 是 拟 合 后 多 项 式 最 高 次 。 

预测 值 计算 命令 : 

Z =polyval(p,x) 

其 中 , р= (aanii) 是 多 项 式 Fxz) = ах" +a," + 
о +а,х +а, ВЖ, х= (on ) 是 某 一 取 值 , Ж 
polyfit 所 得 的 回归 多 项 式 在 * 处 的 预测 值 z。 

例 10.1: 对 表 10-1 中 数据 作 二 次 多 项 式 拟 合 。 

表 10-1 


























解 : (1) MATLAB 命令 
>>х =0:0.1:1; 
>>y=[ -0.447 1.978 3.28 6.16 7.08 7.34 7.66 
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9.58 9.48 9.30 11.2]; 
>> А =polyfit(x,y,2); 
>> Z = polyval( A,x); 

>> Plot(x,y,'k +',x,z,'r 
A=| -9. 8108 20. 1293 


ре 


– 0. 0317 ] 


即 拟 合 后 的 函数 为 : f(x) = – 9. 8108? +20. 1293x – 0. 0317 


拟 合 图 如 图 10-1 所 示 。 














图 10-1 


例 10.2: 通过 实验 测 得 一 次 
物 后 的 血 药 浓度 数据 ， 见 表 10-2, 


性 快速 静脉 注射 300mg 药 





表 10-2 
t/h 0.25 | 0.5 1 1.5 2 3 4 6 8 
y(ug/mL) 19.21 |18. 15 |15. 36 | 14. 10 |12. 89 | 9.32 |7. 45 | 5.24 | 3. 01 


求 血 药 浓度 随时 间 的 变化 规律 y(t 


) 。 


解 : 首先 利用 MATLAB 程序 描 出 血 药 浓 度 与 时 间 的 散 
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点 图 ， 如 图 10-2 所 示 。 

从 图 像 可 以 看 出 ， 散 点 图 大 致 呈 负 指 数 函 数 形态 , 邻 у = 
а ег", 其 中 , t 和 y 分 别 为 时 间 和 血 药 浓度 , a ,b(a,b >0) 为 
待定 系数 。 


20 














图 10-2 


将 其 进行 线性 化 处 理 ， 等 式 y =а e“ 两边 取 对 数 ， 
lny = lna — bt 
&lny=Y, Ina =8,, -b=B,, MZDE: 了 =B +В 
用 MATLAB 程序 为 : 
>>t=[0.250.511.523468]; 
>>y= [ 19.21 18.15 15.36 14. 10 12.89 9.32 7.45 5.24 
3.01]; 
>> Y =log(y); 
>> р = polyfit(t,y,1) 
P= 
-0. 2347 2. 9943 
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即 В, = 2.9943, В, = - 0. 2347. 从 而 a = e® = e” = 
19. 971, b =0. 2347 
则 血 药 浓度 与 时 间 的 关系 为 y = 19. 9714e…”“。 变 化 规 得 
曲线 如 图 10-3 所 示 。 
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图 10-3 
10.2 非 线 性 最 小 二 乘 拟 合 
序 命令 如 下 : 


>> К = lsqcurvefit(' fun ', х0 ,xdate ,ydate ) 
Hh, fun 是 一 个 事先 建立 的 定义 函数 F(x, xdate) 的 M- 文 
件 ， 自 变量 x，x0 是 迭代 初 值 ，xdate ，ydate 为 已 知 数据 点 。 
M 文件 的 编制 方法 : 在 工具 栏 点 击 File 一 一 单 击 New 

一 一 单 击 M-File， 在 Editor-Unitled 中 编辑 M 文件 。 
例 10.3: ЕЮ (ЕЛ А Ш, АИ 
ВИНЫ) PAE: 芍 驶 员 血 液 中 的 酒精 含量 大 于 或 
等 于 20mg/100mL， 小 于 80 mg/100mL 为 饮酒 驾车 ， 血 液 中 
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的 酒精 含量 大 于 或 等 于 80mg/100mL 为 醉酒 加 车。 假设 一 体 
重 约 70kg 的 人 第 一 次 饮酒 后 血液 中 的 酒精 含量 满足 以 下 模 
Ж: х=у(е“ -е), Нину, а, В 为 待定 常数 ,1 为 距 饮 
酒 的 时 间 。 当 此 人 在 短 时间 内 喝 下 2 瓶 啤酒 后 ， 隔 一 段 时 间 
测量 他 血液 中 的 酒精 含量 (mg/100mL) ， 所 得 数据 见 表 10- 
3 。 试 确定 驾驶 员 第 一 次 喝 两 瓶 啤酒 后 血液 中 的 酒精 含量 的 
模型 。 



































表 10-3 (单位 : mg/100mL) 
ИВ |0.25| 0.5 10.75| 1 151 2 125| 3 139514 14.5 | 5 
酒精 含量 | 30 | 68 | 75 | 82 | 82 | 77 |68 | 68 | 58 | 51 50 | 41 


























时 间 /h | 6 7 8 9 10 | 11 | 12 | 13 | 14 | 15 | 16 
酒精 含量 | 38 | 35 28 25 |18 |15 |12 |10 7 7 4 
























































解 : MATLAB 程序 为 : 

首先 在 File 一 New 一 M-File 中 定义 非 线性 函数 fun9 _3. m 
文件 : 
Function y=fun9 _3(a,x) 
у=а(3) * (ехр(1а(1) *х) 一 exp( -a(2) *х)); 
在 命令 窗口 中 运行 : 
>>x=[0.250.50.7511.522.533.544.55678910 
11 12 13 14 15 16] ; 
>> y = [30 68 75 82 82 77 68 68 58 51 50 41 38 35 28 25 18 
15 1210774]; 
>> a = lsqcurvefit(' fun9 _3',[0.2;2;100] ,x,y) 
а =0. 1855 

2. 0080 

114. 4329 
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ВР а =0. 1855, В =2. 0080, у=114. 4329 
由 此 可 见 ， 驾 驶 员 第 一 次 喝酒 后 ， 体 内 血液 中 酒精 含量 与 时 
间 的 关系 为 : 
x =114. 4329 (e "193 — e 200) 
拟 合 图 形 如 图 10-4 所 示 。 
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图 10-4 


对 于 表 10-2 中 的 数据 ， 也 可 以 采用 非 线性 最 小 二 乘 拟 
合 进行 处 理 。 

解 : 使 用 非 线性 拟 合 的 函数 nlinfit( ) ,首先 建立 fun9 _ 
1. m 文件 如 下 : 


function result =fun9 1(a,x) 





result =а(1) * exp( —a(2) *x); 

观察 散 点 图 , 不 妨 取 a=20, H £= 1 В, y=15.36 Ж Б = 
0. 264 , 取 beta0 = [20 0. 264]; 

МАТГАВ 的 程序 如 下 : 

>>t=[0.250.511.523468]; 
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>> у = [19.21 18. 15 15.36 14.10 12.89 9.32 7.45 5.24 
3.01]; 
>> beta0 = [20 0.264]; 
>> beta =nlinfit(t,y,'fun9 _ 1 ',beta0 ) 
beta = 
20. 2413 0.2420 
则 血 药 浓度 与 时 间 的 关系 为 : у =20. 2413е 0 
拟 合 效果 如 图 10-5 所 示 。 














12 `< 
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+ 
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图 10-5 





说 明 : 此 法 与 前 法 虽然 方法 不 同 ， 但 拟 合 的 效果 同样 较 





好 。 

例 10.4: 化 工 生产 中 获得 氯气 的 等 级 y 随 着 生产 时 间 x 
的 增加 而 下 降 ， 假 设 当 x*>=>8 PF, y Ej x 之 间 满 足 如 下 非 线 
性 关系 y=a+(0.49 -a)je "9, Е, а, b 为 待定 系数 ， 
现 收集 了 44 组 数据 见 表 104， 试 确定 氯气 与 生产 时 间 的 关 
系 。 
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表 10-4 























у |0.42 0.4 0.4 0.41 0.4 0.41 0.41 0.4 0.4 10.4 10.38 0.41 











у |0.4 0.4 0.41 0.38 0.4 0.4 0.39 0. 39 






































解 : 本 题 为 非 线性 拟 合 模型 。 要 求 a、5b 的 值 ， 可 调用 
非 线性 拟 合 函数 nlinfit( ) 。 

首先 定义 非 线性 函数 fan9 _ 2. m 文件 : 
function result =fun9 2(beta0 ,x) 
а = Беіа0(1) ; 
b = beta0 (2); 
result =a + (0. 49 —a) * exp( -bx (х-8)); 

然后 在 窗口 命令 中 运行 以 下 程序 : 
>>х= [8 8 10 10 10 10 12 12 12 12 14 14 14 16 16 16 18 18 
20 20 20 20 22 22 24 24 24 26 26 26 28 28 30 30 30 32 32 34 
36 36 38 38 40 42]; 
>>y = [0. 49 0. 49 0. 48 0. 47 0. 48 0. 470. 460. 460. 45 0. 43 
0. 45 0. 43 0. 43 0. 44 0. 43 0.43 0. 46 0. 45 0.42 0. 42 0. 43 
0.41 0.410.400. 420.40. 40. 410.40.410. 41 0.40.40. 4 
0. 38 0. 41 0. 4 0. 4 0. 41 0. 38 0.4 0.4 0. 39 0.39]; 
>> beta0 = [0.3 ,0. 02]; 
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>> betafit =nlinfit(x,y,'fun9 _ 2. m ',beta0 ) ; 
Betafit = 
0. 3904 0.1028 
Вр: а =0. 3904, 6 =0. 1028 
所 以 模型 为 : у = 0. 3904 + (0. 49 - 0. 3904) e 010808) = 
0. 3904 +0. 0996e 7° 08-9 
拟 合 效果 如 图 10-6 所 示 。 
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图 10-6 











9110.5: 我 国 是 一 个 人 口 大 国 ， 控 制 人 口 增长 是 我 国 
的 一 个 基本 国策 ， 有 效 控制 人 口 增长 的 前 提 是 要 认识 人 口 数 
量 的 变化 规律 ， 通 过 建立 人 口 模型 ， 可 以 作出 较为 准确 的 预 
报 。 表 10-5 给 出 了 1949—1994 年 我 国人 口 数 据 资料 。 
表 10-5 



































Ф х; | 1949 | 1954 | 1959 | 1964 | 1969 | 1974 | 1979 | 1984 | 1989 | 1994 





АНЖУ; | 5.4 | 6.0 167 17.0 8.1 91 98 | 10.3 | 11.3 11.8 
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试 分 析 我 国人 口 增长 规律 ， 建 立 我 国人 口 模型 ， 并 计算 
1999 年 我 国人 口 数 量 。 
ВЕ: 先 利 用 MATLAB 作出 已 知 数据 的 散 点 图 ， 如 图 10- 


7 所 示 。 
>> x = [ 1949 1954 1959 1964 1969 1974 1979 1984 1989 


1994 ] ; 
>>у= [5.4 6.0 6.7 7.0 8.1 9.1 9.8 10.3 11.3 


11.8]; 
>> ро (х,у,'*!) 











+ 











1 П 1 1 П L 1 L L 
145 1950 1955 1960 1965 1970 1975 1980 1985 1990 1995 


图 10-7 








模型 一 : 用 线性 函数 拟 合 。 设 我 国人 口 数量 满足 以 下 线 
y=a+bx 

Я, a, b 为 待定 系数 。 

用 MATLAB 算出 c = -283.2320， =0. 1480， 从 而 得 到 我 


国 的 人 口 模型 为 
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y = -283. 232 +0. 148x 
模型 二 : 用 指数 函数 拟 合 。 НАНИЖЕ а: 
у=ае" 
ЖФ, а, b 为 参数 。 模 型 两 边 同时 取 对 数 得 : lny = lna + bx 
下 面 用 简单 的 线性 最 小 二 乘法 拟 合 出 a、5， 用 MAT- 
LAB 算得 a =2. 33, Ь =0. 0179 
即 у=2. 33e 0 
用 MATLAB 求解 如 下 : 
>> x = [ 1949 1954 1959 1964 1969 1974 1979 1984 1989 
1994]; 
>>у=[5.4 6.0 6.7 7.0 8.1 9.1 9.8 10.3 11.3 
11.8]; 
>>а= polyfit(x,y,1); 
>>х| = [1949:5 ;1994 ]; 
>>yl =а(2) +а(1) «х1; 
>> Ь =polyfit(x,log(y) ,1); 
>> у2 = ехр(Ь(2)) жехр(Ъ(1) *х1) 
>> plot(x,y,'* ') 








>> hold on 
>> plot(xl ,yl ,'--r') 
>> hold on 
>> plot(x1 ,y2 ,'-k ' 
>> legend(' 原 离散 点 ', 模 型 一 曲线 ,模型 二 曲线 ') 
拟 合 后 得 到 两 个 图 形 ， 如 图 10-8 所 示 。 用 两 个 模型 分 
别 计 算 相 应 年 度 的 人 口 数 见 表 10-6。 计 算 模 型 一 的 误差 平方 
Ж 0, = У (у, -а – Бх,)?, 模型 二 的 误差 平方 和 0, = X (y, - 


bx; N 2 
ae 
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10 上 4 
gr J 
8 上 J 
E- J 
Er 模型 一 曲 
g€ 一 一 模型 二 曲线 
т 1950 1955 1960 1965 1970 1975 1980 1985 1990 1995 
图 10-8 
表 10-6 
Е) х, | 1949 | 1954 | 1959 | 1964 | 1969 | 1974 | 1979 | 1984 | 1989 | 1994 
АНЖУ, | 5.4 | 6.0 | 6.7 | 7.0 | 8.1 | 9.1 | 9.8 | 10.3 | 11.3 | 11.8 
模型 一 | 5.24 | 5.97 | 6.70 | 7.43 | 8.16 | 8.90 | 9.62 | 10.36 11. 09 |11. 82 
误差 | 0.16 | 0.03 | 0.00 |-0.43|-0.06| 0.20 | 0.18 |-0. 06] 0.01 |-0. 02 
模型 二 | 5.55 | 6.06 | 6.62 | 7.23 | 7.90 | 8. 64 | 9.44 | 10.31 | 11. 26 |12. 31 
误差 |-0. 15 -0.06 0.08 |-0.23 0.20 | 0.46 | 0.36 |0. 01-0. 13| =0. 51 
结果 分 析 : 


(РО, =0. 2915 <0. 7437 = 0,， 由 最 小 二 乘 指标 判断 可 
知 ， 线 性 模型 更 适合 
G@ 用 以 上 两 种 模型 预测 我 国 1999 年 的 人 口 分 别 为 
12.55 亿 ，13. 43 亿 。 
@ 我 国人 口 白皮书 披露 2005 年 我 国有 13.3 亿 人 口 ， 











国人 口 增长 。 
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2010 年 有 14 亿 人 口 。 
模型 一 得 到 的 预测 数据 为 : 2005 年 13. 43 亿 ，2010 年 14. 16 
亿 。 
模型 二 得 到 的 预测 数据 为 : 2005 年 14. 94 亿 ，2010 年 16. 33 
亿 。 

通过 实际 比较 可 以 看 出 ， 模 型 一 更 接近 我 国人 口 实际 数 
据 。 

人 口 增长 受 诸多 因素 影响 ， 为 了 使 人 口 预报 〈 特 别 是 
长 期 预报 ) 更 好 地 符合 实际 情况 ， 必 须 综合 考虑 诸多 方面 
的 因素 ， 不 断 修改 模型 。 


Z 训 w 


1. 用 最 小 二 乘法 求 一 个 形 如 y =a+ bx 的 经 验 公 式 ， 使 它 与 下 表 
所 示 数 据 拟 合 。 

















拟 合 数据 表 
x 19 25 31 38 44 
у 19.0 32.3 49 73.3 97.8 


2. 某 乡镇 企业 1990—1996 年 的 生产 利润 见 下 表 : 








年 份 1990 1991 1992 1993 1994 1995 1996 
利润 /万 元 70 122 144 152 174 196 202 





试 预测 1997 年 和 1998 年 的 利润 。 


第 11 章 ”层次 分 析 法 的 
MATLAB 程序 


11.1 层次 分 析 法 简介 


1. 建立 层次 结构 模型 

通过 调查 人 研究， 将 目标 准则 体系 所 包含 的 因素 划分 为 不 
同 层次 ， 如 目标 层 、 准 则 层 、 方 案 层 等 ， 并 采用 图 表 形 式 ， 
标明 递 阶 层次 结构 模型 。 

2. 构造 判断 矩阵 

按照 层次 结构 模型 ， 将 每 一 层 元 素 以 相 邻 上 一 层 某 元 素 为 
准则 ， 进 行 成 对 比较 ， 并 按 1 ~9 标 度 方法 ， 构 造 判断 矩阵 。 

具体 操作 : 设 B 层 有 三 个 元 素 B., В,, В,, # W JH 
上 一 层 的 元 素 4;， 则 其 判断 矩阵 为 : 

A, B, B, B, 




















用 左边 的 元 素 比 上 方 的 元 素 ， 按 1 ~9 标 度 方法 : 同等 
取 1， 稍 重要 取 3， 重 要 取 5， 很 重要 取 7， 非 常 重要 取 9; 
重要 性 介 于 他 们 之 间 时 取 2、4、6、8。 

3. 求 出 每 一 个 成 对 比较 矩阵 的 最 大 特征 值 ， 作 单项 一 
致 性 检验 
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Aa 
具体 操作 : 一 致 性 指标 公式 为 : CI = 
其 中 , 和 A 分 别 为 矩阵 的 阶 数 和 最 大 特征 值 。 
随机 一 致 性 指标 RI 由 表 11-1 确定 。 
表 11-1 
п 112|3 ЕЕЗ ЕЕ ЕЕ ЕТ: 





КІ 0 0 0.58 | 0.9 11.12 | 1.24 | 1.32 | 1.41 | 1.45 | 1.49 





CI 
从 而 得 到 一 致 性 比率 СЕ = 7 о #1 СК <0.1， 则 通过 一 


致 性 检验 ;否则 ， 需 重新 调整 比较 矩阵 。 

MATLAB 的 计算 程序 为 : 
>>P=[……] ”% 输 入 成 对 比较 矩阵 
>>[5,Т] =eig(P),Pmax =T(1,1),CI = (Pmax -n)/(n-— 
ВТ ,CR = СИНІ 

% 求 P 的 特征 值 T(1,1), 作 一 致 性 检验 

4. 作 总 体 一 致 性 检验 

具体 操作 : 

9 求 出 准则 层 对 目标 层 的 成 对 比较 矩阵 的 最 大 特征 值 对 
应 的 归 一 化 特征 向 量 W, ; 

@ 提 取 方 案 层 对 准则 层 的 所 有 CI, 和 RI, 分 别 构成 Ze 
和 Zn; 








Ze: W. 
@@ 总 体 一 致 性 比率 为 CR= 7y 
# CR<0.1， 则 通过 一 致 性 检验 ; 否则 需 重新 调整 比较 
和 矩阵。 
MATLAB 计算 程序 为 : 
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>>Wl=S(1,l)/sum(S(1,1)),ZC=[CH,CD ，……],ZR 
= [RI ,RI2 ,……… ] 
>> СВ =ZC * \1/ (78 * W1) 

5. 计算 层次 总 排序 权 值 

具体 操作 : 提取 方案 层 对 准则 层 的 所 有 成 对 比较 和 矩阵 的 
最 大 特征 值 对 应 的 归 一 化 特征 向 量 ， 构 成 矩阵 或 块 对 角 和 矩阵 
W,， 则 总 排序 权 值 为 : W =W, - W, , 

MATLAB 计算 程序 举例 : 
>>W2=[S1(1,1)/sum(S1(1,1)),S2(1,1)/sum($2(1, 





>>W=W2 * Wl 


11.2 应 用 


例 11.1: (ATER) 菜单 位 拟 从 3 名 干部 中 选拔 领 
T, 标准 有 : 健康 状况 、 业 务 知识 、 写 作 能 力 、 口 才 、 政 策 
水 平 、 工 作 作风 。 试 用 层次 分 析 法 对 3 人 作 综 合 评 佑 、 量 化 
排序 。 

解 : 1. 建立 层次 结构 模型 


目标 层 
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2. 构成 比较 矩阵 











L.E i 1 1 
2 
1 
1 1 2 Jo а 
2 
1 1 
j eS] 8: == 
> 2 
1 1 1 1 1 
4 4 5 3 3 
1 
Je 3; < 1 1 
3 
оо 1 1 
健康 状况 业务 知识 
гг тт 
4 2 5 8 
4 1 3 В, = 1 
Е 
2 工 1 
3 8 2 1 
写作 能 力 口才 
1 
1 1 1 — 5 
1 — 一 3 
3 5 
B,=|3 1 7 
1 1 
1 1 
5 1 1 Е, $ 
5 7 
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政策 水 平 工作 作风 
1 
1 1 1 — 1 
7 
1 1 
2 Bs=|I7 1 5 
— — 1 l 
6 7 1 — 1 
5 


3. 单项 一 致 性 检验 (用 MATLAB 程序 检验 ) 
>>Р=[1,1,1,4,1,1/2;1,1,2,4,1,1/2;1,1/2,1,5,3,1/2; 
1/4 ‚1/4 ,1/5,1,1/3 ,1/3;1,1,1/3,3,1,1;2,2,2,3,1,1] 
>>[5,Т] =eig(P),Pmax =Т(1,1),С1= (Pmax -6)/(6 - 
1) ,ВТ=1. 24, СВ =СИВГ СВ =0. 0678 | 
>> ВІ = [1,1/4,1/2;4,1,3;2,1/3,1],В2 = [1,1/5 ,1/8;5, 
1,1/2;8,2,1] 
>> ВЗ = [1,1/3,1/5;3,1,1;5,1,1],В4 = [1,1/3,5;3,1,7; 
1/5 ,1/7,1] 
>> В5 = [1,1,6;1,1,7;1/6,1/7,11],В6 = [1,1/7,1;7,1,5; 
1,1/5,1] 
>>[51,Т1] = ев (ВІ), ВІ тах = Т1(1,1),СП = (ВІ тах — 
3)/(3-1),ВП =0. 58,СКІ =СПИВИ 
>> [52,Т2] = eig( B2) ,B2max = T2 (1,1) ,С12 = ( B2max - 
3)/(3-1) ,RD =0. 58 ,CR2 = СІ2/КІ2 
>> [53 ,Т3 | =eig(B3),B3max = ТЗ (1,1), СІЗ = ( B3max - 
3)/(3-1) , КІЗ =0. 58, СВЗ = СІЗ/КІЗ 
>> [54,14] = eig ( B4 ) ,B4max = T4 (1,1), СМ = ( B4max - 
3)/(3-1) КИ =0. 58 ,CR4 = CI4/RI4 
>> [55,15] = eig( B5 ) ,B5max = Т5 (1,1), С15 = ( B5max - 
3)/(3-1) ,RIS =0. 58 ,CR5 = С15/КІ5 
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>>[S6,T6] =eig(B6),B6max=T6(1,1),CI6 = ( B6max – 
3)/(3 -1) ,RI6 =0. 58 ,CR6 = CI6/RI6 
>> ZB =[CR1,CR2,CR3,CR4,CR5,CR6 ] 
ZB =0. 0158 0. 0048 0. 0251 0. 0559 0. 0023 0. 0109 
所 有 值 都 小 于 0. 1 ,通过 一 致 性 检验 。 
4. 总 体 一 致 性 检验 
>> Wl =|[S(:,l)Z/sum(s(:,1))] 
>> ZC = [| CH ,СІ2,СІЗ, СІМ, С15,С16 ] ,ZR = [ВП , КІ2, КІЗ, 
RI4 ,RIS ,RI6 ] 
>>a=ZC * W1,b=ZR *W1,CR =а/Ь 
>> СВ =0. 0142 
由 于 CR 的 值 小 于 0.1, 所 以 总 体 通过 一 致 性 检验 。 
5. 排序 
>> Wl =[S(:,l)/Zsum(s(:,1))] 
>> № =[S1(:,1)/sum(sl(:,1)),S2(:;,1)/sum(s2(:, 
1)),S3(:,1)/sum(s3(:,1)),S4(:,1)/sum(s4(:,1)),S5 
(:,1)/sum(s5(:,1)),S6(:,1)/sum(s6(:,1))] 
>> W = W2 * WI 





W= 0.1694 
0. 5394 
0. 2957 
从 结果 比较 ， 故 应 选 第 二 个 人 。 


Z 训 w 
小 李 是 某 大 学 的 应 届 毕 业 生 ， 参 加 多 家 企业 的 招聘 会 后 ， 有 甲 、 
乙 、 丙 三 个 单位 愿意 录用 他 。 请 你 为 他 提供 建议 ， 以 选择 合适 的 单 


位 。 
























































第 12 章 M 文件 和 编程 


假如 读者 想 灵活 运用 MATLAB 去 解决 实际 问题 ， 想 充 
分 调动 MATLAB 的 科学 技术 资源 ， 想 理解 MATLAB 版 本 升 
级 所 依仗 的 基础 ， 那 么 本 章 内 容 将 十 分 有 用 。 


12.1 编程 入 门 





通过 M 脚本 文件 ， 画 出 下 列 分 段 函数 所 表示 的 曲面。 
0. 5457е 79 523.75 - 1.50) 


xita >] 
p(xi,x,) = {0. 7575е 7-61 


-1 <, +х, <1 
0 5457679755 -3-759 +1. 59 


(1) 编写 М 脚本 文件 的 步 又 


х 二 

















MATLAB Editor/Debugger - [Untitledl] 





4 File Edit View Debug Tools Window Help 














=la] x| 
Dk | х ве | ë аа 
ЕЕ | | Stack: | El 








ГІ 
| [J песи 








Line 1 7:52 РМ 4 
图 12-1 





ө 点 击 MATLAB 指令 窗 工 具 条 上 的 New File Ех O, 
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可 打开 如 图 12-1 所 示 的 MATLAB 文件 编辑 调试 器 MATLAB 
Editor/Debugger。 其 窗口 名 为 untited ， 用 户 可 在 空白 窗口 中 
编写 程序 。 比 如 输入 如 下 一 段 程序 : 
[2х1. m] 
a=2;b =2; 
clf; 
х= -а:0.2;а;у= -b:0.2:b; 
for i = 1 еп (у) 
for j = 1 :Іепећ (х) 
if х()) +у(1) >1 
2(1,)) =0. 5457 x exp( – 0. 75 + у(1)^2 -3.75 *x(j)2 

е 

elseif x(j) +у(1) < = -1 

2(1,)) =0. 5457 x ехр( – 0. 75 + у(1)^2 – 3. 75 *x(j)2 

+1.5жх(])); 

else 2(1,)) =0. 7575 * ехр( —у(1)^2 -6. *х(])^2); 

end 

end 

end 
axis([ —a,a, -—b,b,min(min(z)) ,max(max(z))]); 


colormap ( flipud( winter) ) ;surf( x,y,z); 


e。 点 击 编辑 调试 器 工具 条 图 标 图 ， 在 阐 出 的 Windows 
标准 风格 的 “保存 为 ”对 话 框 中 ， 选 择 保 存 文件 夹 ， 键 入 
新 编 文件 名 (如 zxl1)， 点 动 【 保 存 】 键 ， 就 完成 了 文件 保 
To 

(2) 运行 文件 

e 使 xl.m 所 在 目录 成 为 当前 目录 ,或 让 该 目录 处 在 
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MATLAB 的 搜索 路 径 上 。 
e 然后 运行 以 下 指令 ， 便 可 得 到 图 12-2 所 示 图 形 。 


zx] 





图 12-2 


12.2 MATLAB 控制 流 


1，for 循环 结构 
例 12.1: 一 个 简单 的 for 循环 示例 。 
fori=1:10; % 依次 取 1,2,…,10 
x(i) =i; APEA i 值 ,重复 执行 由 该 指令 构成 的 循环 体 


end; 
x % 要求 显示 运行 后 数组 x 的 值 。 


1 2 3 4 5 6 7 8 9 10 
2. while 循环 结构 
例 12.2: Fibonacci 数组 的 元 素 满足 Fibonacci 规则 : а, 
=G, +G,.15 (k =1,2,:-‹) ; Н. G, = а, = 1. 现 要 求 该 数组 中 
第 一 个 大 于 10000 的 元 素 。 
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а(1) =1;а(2) =1;i=2; 
while a(i) < =10000 
a(i+1) =а(1-1) +а(1); 狗 当 现 有 的 元 素 仍 小 于 
10000 时 , 求解 下 一 个 
元 素 。 





21 
ans = 
10946 
3. if-else-end 分 支 结 构 
例 12. 3 :一 个 简单 的 分 支 结构 。 


cost = 10 ; number = 12; 





if number > 8 
sums = number * 0. 95 x cost; 
end ,Sums 
sums = 
114. 0000 
$1 12.4: 用 for 循环 指令 来 寻求 Fibonace 数组 中 第 一 个 
大 于 10000 的 元 素 。 


п =100;а =опез(1,п); 





for 1=3:n 
a(i) =а(1-1) +а(1-2); 
if a( i) > = 10000 
а(1), 
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break; 多 跳 出 所 在 的 一 级 循环 
end ; 
end ,1 
ans = 
10946 
i= 
21 
4. switch-case 结构 
BJ 12.5: 学 生 的 成 绩 管理 ,用 来 演示 switch 结构 的 应 








用 。 
clear; 

% 划 分 区 域 : 满分 (100), 优秀 (90 -99), 良好 (80 ~ 89), 
及 格 (60 -79), 不 及 格 ( <60)。 

fori=1:10;ali =89 +1;611} =79 +i;c{i} =69 +1;411} = 
59 +1;епа;с= [4,с]; 

Name = !' Jack ',' Marry ',' Peter ',' Rose ', Тот '|; 狗 元 胞 数组 
Mark = {72 ,83 ,56 ,94 ,100} ;Rank = cell(1,5) ; 

P 创建 一 个 含 5 个 元 素 的 构架 数组 S, 它 有 三 个 域 。 

S =struct(' Name ', Маше ,' Marks ', Mark ,' Rank ', Rank ) ; 

和 根据 学 生 的 分 数 , 求 出 相应 的 等 级 。 











for i=1:5 
switch S(i). Marks 
case 100 % 得 分 为 100 时 
5(1). Rank = 满分; 。% 列 为 满分 等 级 
case a % 得 分 在 90 和 99 之 间 


S(i). Rank =; 和 列 为 优秀 等 级 
case b % 得 分 在 80 和 89 之 间 
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S(i). Rank =' 良 好 ; 狗 列 为 良好 等 级 


case c % 得 分 在 60 和 79 Z [j 
5(1). Rank = ЛО; % Я МЕХ 
otherwise % 得 分 低 于 60. 
S(i). Rank = 不 及 格 ; % 列 为 不 及 格 等 级 
end 
end 





% 将 学 生 姓 名 ,得 分 ,登记 等 信息 打印 出 来 。 
disp([ 学 生 姓名 ', 得 分 ,等 级 ] ) ;disp(' ') 
fori=1:5; 

disp ([S (i). Маше, blanks (6), num2str (5 (1). Marks ) , 
blanks(6),S(i). Rank] ); 
end; 
学 生 姓名 ”得 分 ”等 级 

Jack 72 及 格 


Tom 100 满分 
5. try-catch 结构 
例 12. 6: try-catch 结构 应 用 实例 。 
clear,N =4;A =magic(3); % 设置 3 行 3 ЯНА 
try 








А_М=А(М,:), % ХА 的 第 NN 行 元 素 
catch 
А end=4A(end,:), % 如 果 取 A(N,:) 出 错 , 则 改 


ЖА 的 最 后 一 行 
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end 
lasterr % 显示 出 错 原因 
А _end = 
4 9 2 
ans = 


Index exceeds matrix dimensions. 
6. M 文件 的 一 般 结 构 
例 12.7: M 函数 文件 示例 。 
[ circle. m ] 
function sa =circle(r,s) 


% CIRCLE plot a circle of radii r in the line specified by s. 





9% т 指定 半径 的 数值 
%s 指定 线 色 的 字符 串 
% за 圆 面积 


% circle( r) 利用 蓝 实 线 画 半 径 为 > 的 圆周 线 . 

% circle(r,s) 利用 字符 串 s 指定 的 线 色 画 半 径 为 > 的 
圆周 线 . 

% sa=circle(T) 计算 圆 面积 ,并 画 半径 为 > 的 蓝 色 圆 面 . 

% sa=circle(r,s) 计算 圆 面积 ,并 画 半径 为 > 的 * 色 圆 面 . 

if nargin >2 


error( 输 入 宗 量 太 多 。'); 





end; 

if nargin = =1 
s='b'; 

end; 


clf; 


. 120 . MATLAB 基础 与 应 用 





1=0:р1/100:2 * pi; 
x=r*exp(i*t); 
Ш nargout = = 0 
plot(x,s); 
else 
sa=pi*r*r; 
fill(real(x),imag(x),s) 
end 
axis(' square ') 
7. 内 联 函 数 
例 12. 8: 演示 内 联 函 数 的 第 一 种 创建 格式 。 
clear,F1 = inline('sin(rho)/rho') % 第 一 种 格式 创建 内 联 函 数 
ЕІ = 
Inline function: 


ЕІ (rho) = sin(rho)/rho 








fl =F1(2) % 内 联 函 数 的 一 种 使 用 方法 
П = 
0. 4546 
КЕ! = vectorize ( F1 ) % 产生 适 于 “数组 运算 ”的 
内 联 函 数 

xx=|0.5,1,1.5,2];ffl =FF1(xx) 
FFI = 

Inline function: 

ЕЕ! (rho) =sin(rho). /rho 
ffl = 


0.9589 0.8415 0. 6650 0. 4546 
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实 训 题 

1. 鸡 兔 同 笼 ， 共 有 36 个 头 、100 只 脚 ， 问 鸡 兔 各 有 多 少 ? 

2. 公元 5 世纪 我 国 古代 数学 家 张 丘 建 在 《 算 经 》 一 书 中 提出 了 
“ 百 鸡 问题 ” : 鸡 翁 一 值钱 五 ， 鸡 母 一 值钱 三 ， 鸡 雏 三 值钱 一 。 百 钱 买 
AX, MSA, В}, ЗЕЕ ЛИ? 

3. 猴子 分 桃 ， 五 只 猴子 ， 第 一 只 将 桃子 平均 分 成 五 份 ， 多 了 一 个 
扔 掉 ， 并 拿 走 一 份 ， 第 二 只 猴子 ， 把 剩 下 的 平均 分 成 五 份 ， 又 多 了 一 
个 扔 掉 ， 拿 走 一 份 。 第 三 、 第 四 、 第 五 只 猴子 都 这 样 做 ， 问 原来 一 共 
有 多 少 个 桃子 ? 

4. 求解 猴子 吃 桃 问题 。 猴 子 在 第 一 天 摘 下 若干 个 桃子 ， 当 即 就 吃 
了 一 半 ， 又 感觉 不 过 净 ， 于 是 就 多 吃 了 一 个 。 以 后 每 天 如 此 ， 到 第 10 
天 再 想 吃 时 ， 却 发 现 就 只 剩 下 了 一 个 桃子 。 请 编程 计算 第 一 天 猴子 摘 
的 桃子 个 数 。 
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第 1 章 实 训 题 参考 答案 
、1996 ，， 
L 表达 式 8 的 计算 结果 。 


解 : MATLAB 程序 如 下 : 





>> 1996/18 
ans = 
110. 8889 
° 、 2sin0.3 Аа 
2. 请 输入 表达 式 一 一 一 ,并 运行 。 
1+/5 


ff: MATLAB 程序 如 下 : 
>>2 * sin(0.3 *pi)Z/(1 +sqrt(5)) 


ans = 


0. 5000 
1 4 7 
3. 请 输出 矩阵 |2 5 8 
3 6 9 


解 : MATLAB 程序 如 下 : 
>> [1,4,7;2,5 ,8;3,6,9] 


ans = 
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第 2 章 实 训 题 参 考 答案 
1. Е (х) =4x -5 /х+1,Ж КЗ) ШИ, 
解 : MATLAB 程序 如 下 : 
>> syms x 
>>х=3; 
>>f=4 *х3 – 5 ж sqrt(x) +1 
f= 
100. 3397 


2. ВЖЕ (х) = x° -Ix +3. 45sinx – 2соѕх, 3K (3) АЈ 


解 : MATLAB 程序 如 下 : 
>>х=3; 
>> Ё = х4 – х^(1/4) +3. 45 * sin(x) —2 * cos(x) 
f = 
82. 1508 
3. HIH pK у = tanx 的 图 形 。 
解 : MATLAB 程序 如 下 : 
ezplot(' tan(x)') 
图 形 如 图 附 -1 所 示 。 
4. 在 区 间 0<x<8 绘制 y=2e “sin(x) 的 图 形 。 
解 : MATLAB 程序 如 下 : 
>>f='2 жехр( -х). * sin(x) '; 
>> fplot(f, [0 8]); 
>>title(f) ,xlabel(' x ') 
图 形 如 图 附 -2 所 示 。 
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. 124. 
tan(x) 
| | | 
| | 
| 
! | | 
2 / p: РА и 
TH 0). АД у a 
и А и i 
-2 / / [ 
| | | 
-4 | | | 
о 
— — 0 
图 附 -1 
2*exp(-x).*sin(x) 
0.7 ' 




















0.1 i 
0 2 М 
x 
图 附 -2 
х+1,х<1 


解 : MATLAB 程序 如 下 : 
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>>xl = -3:0.1:1;x2 =1:0.1:3; 
> yl =xl +1;y2 =1+1./х2; 
>> plot(x1 ,у1,х2,у2,г-!) 














图 形 如 图 附 -3 所 示 。 
2 = 
z 
15 А ЗЕЕ 
1L 
У 
051 и 
4 
0 Г Р, 
РА 
05| РА 
1р J 
154 
2“ 1 
-3 2 1 0 1 2 3 


图 附 -3 


6. 画 出 y=sin(x) „y =s [+ | +2,y = cos(x) HXT EE 


图 。 
解 : MATLAB 程序 如 下 : 
>> clc ,clear 
>>x= —2 *pi:0.1:2 * pi; 
>>yl =sin(x);y2 =sin(x +pi⁄3) +2;y3 =cos(x); 
>> plot(x,yl,'. =); 
>> hold оп % 图 形 保 持 命令 
>> plot(x,y2,'* —') ;plot(x,y3,'-0'); 
图 形 如 图 附 4 所 示 。 
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第 3 章 实 训 题 参考 
1. КВ: 

(1) lm 

解 : MATLAB 程序 如 下 : 
>> syms x 

>> limit(sin(2 x x)/x,x,0) 


й 


案 





ans = 


2 
(2) 加 | 2 ) 





SE 2 -1 
解 : MATLAB 程序 如 下 : 


>> syms x 
>> limit( (1/ (x-1) -2/(x°2 -1)),x,1) 


ans = 
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1⁄2 


(3) lim 


x—> o 


解 : MATLAB 程序 如 下 : 
>> syms x 
>>limit(((2*x+1)/(2*x-1)) (x+1),x,inf) 
ans = 
exp(1) 
(4) ба 
550-5 
解 : MATLAB 程序 如 下 : 
>> Шии(1/Их,х,0 ей’) 
ans = 


— Inf 





(2 а 


2х -1 


(5) a Es 


x>0 siny 


ff: MATLAB 程序 如 下 : 
>> limit(1/sin(x)) 


ans = 


NaN 
x3 
i 3⁄2 х +1015 . 
2 у= (т) А. 
2X +х +1 х0 


解 : MATLAB 程序 如 下 : 
>> syms x y 


> y=((3*x2-—x+1)/(2*x2 +х+1))^(х^З/ (1 - 
x)); 


>> limit(y,x,0) 
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3. 求 函数 的 导数 : 

(1) y =xe” 

解 : MATLAB 程序 如 下 : 
>> syms x y 

>>y =x ж ехр(х^2); 

>> diff( y) 

ans = 


ехр(х^2) +2 * х^2 * ехр(х^2) 


(2) у= cos yx 

解 : MATLAB 程序 如 下 : 

>> syms x y 

>> у = cos(sqrt(x)); 

>> diff( y) 

ans = 

—sin(x(1⁄2))/(2 ж х^(1/2)) 

(3) y =e *tan2x 

解 : MATLAB 程序 如 下 : 

>> syms x y 

>>у=ехр( —3 *x) *tan(2 жх); 

>> diff( y) 

ans = 

(2 *tan(2 жх)^2 +2) /ехр(3 * x) - (3 ж tan(2 *x))/ 
exp(3 * x) 

lnx 


(4) y= > 


2 
x 
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解 : МАТГАВ 程序 如 下 : 
>> syms x y 

>> у =log(x)/x2; 

>> diff( y ,2 ) 

ans = 

(6 x Іов(х) )/х4 -5/x4 
4. 求 不 定 积分 : 

(1) f esimxds 

解 : MATLAB 程序 如 下 : 
>> syms ху 
>>у=ехр(х) * (sin(x))2; 
>> шКу) 


ans = 


— (ехр(х) * (cos(2 + х) +2 *sin(2 жх) -5))/10 
(2) | (Ve +x) ахах 
解 : MATLAB 程序 如 下 : 


>> syms x y 

>>y= (x0.5 +x) * log(x); 

>> int(y) 

ans = 

(x x (log(x) -1⁄2))⁄/2 + (2 ж х^(3/2) * (log(x) -2/ 
3))/3 

(3 ах р 


1 + cosx 


解 : MATLAB 程序 如 下 : 
>>f=(1 +sin(x)) * ехр(х)/(1 +cos(x)); 
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>> int(f,'x') 

(ехр(х) * sin(x) )/(cos(x) +1) 
5. ЖЕ 

(1) [, (3x – 5 )агссоѕхах 


解 : MATLAB 程序 如 下 : 

>> syms x y 

>>у= (З*х-5) *acos(x); 
>> шКу,0,1) 

ans = 


(3 *pi)/8 -5 


(2) р ТИ 


解 : МАТГАВ 程序 如 下 : 
>>int(sqrt(1 – sin(2 * x) ),'x','0','(pi/2)') 
ans = 

2 x2^(1/2) -2 

6. PERAS) =x e ИН. 
ff: MATLAB 程序 如 下 

>> syms x y 

>>y =x?2 *exp( —x); 

>> dy =diff( y) 

dy = 

(2 xx)/exp(x) – х^2/ехр(х) 
>> solve( ду) 


апѕ = 
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2 
>> ddy = diff( dy) 
ddy = 
2/ехр(х) - (4 жх)/ехр(х) + х^2/ехр(х) 
>> subs( ddy,x,0) 
ans = 
2 
>> subs ( ddy,x,2) 
ans = 
– 0. 2707 
>> ѕирѕ(у,х,2) 
апѕ = 
0. 5413 
>> ѕирѕ(у,х,0) 
апѕ = 
0 
因为 x=2 代入 二 阶 导 的 值 为 负 , 所 以 x=2 是 极 大 值 点 ， 
极 大 值 为 0.5413 ,x =0 代入 二 阶 导 的 值 为 正 ,所 以 x =0 是 极 
小 值 点 , 极 小 值 为 0。 
7. К y =x -2x +5 在 区 间 [ -2,2] 内 的 最 小 值 。 
解 : MATLAB 程序 如 下 : 
>> х = fminbnd ('х^4 -2 ж х^2 +5', -2,2) 


x = 











1. 0000 
>>y=subs('x4 -2 ж х^2 +5',x,1) 
у = 
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4 
所 以 最 小 值 为 4 


第 4 章 实 训 题 参 考 答案 
1. 求 微分 方程 y =1+y 的 通 解 。 
解 : MATLAB 程序 如 下 : 
>> dsolve('Dy =1 + y^2 ') 
ans = 
tan(t+c1) 
2. 求 常 微分 方程 +y + (х – 2у)у =0 的 通 解 。 
解 : MATLAB 程序 如 下 : 
>>a='x2 +у+ (х -2 жу) є Dy =0'; 
>> dsolve(a,'x') 
ans = 
х/2 + ( (4 ж x3)/3 + 
x2 + С2)^(1/2)/2 
x/2 — ( (4 x x'3)/3 + 
x2 + С2)^(1/2)/2 


第 5 章 实 训 题 参 考 答案 
ПЕ: 
1. 输入 矩阵 4=|4 5 
7 8 
解 : MATLAB 程序 如 下 : 
>>а=[123;456;789] 
8 


1 2 3 


3 
6 
9 
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4 6 
7 8 9 
2 1 3 4 25 
2. ФА = ,B= ФО) A+B; 
人 人 4 (5 -3 | 019 
(2) A-B. 


ЙЕ: MATLAB 程序 如 下 : 
>>а=[2 1 3;342]; 
>>Ь= [42 5;2 -3 4]; 


>>a+b 
ans = 
6 3 8 
5 1 6 
>a-b 
ans = 
=s 1 й 
1 7 -2 
>a 
з 设 4= 3 Е, ,).в= 2 ll 求 (1) 24; 
3 1 -2 
-5 0 
(2) AB. 


解 : MATLAB 程序 如 下 : 
>>a=[2 -10 ;31 -2]; 
>>b=[1 -3;21; -5 0]; 
>>2 жа 


апѕ = 
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>>a*b 
ans = 
0 -7 
15 -8 
2 1 0 
4. ЖА=|0 1 -3 |,0(1)А4'; (2)А'; (3)А*; (4) 
1 0 2 


REE A 的 秩 。 
解 : MATLAB 程序 如 下 : 
>>А=[210;01 -3;102]; 
>> шу(А) 





ans = 


43 -3 
-3 1 -9 
4 1 4 
>> rank(A) 


ans = 
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3 
5. WA FIB 是 满足 关系 式 4B =А +2B ЖИ, ЯТА 
4 2 3 

=| 1 1 0 |, 求 矩阵 B。 
-1 2 3 


f: 

分 析 :B = (A -2Е) А, НЕ НЕ 
MATLAB 程序 如 下 : 

>>А=[423 ;110;-123]; 


>> В = іпу(А -2 * еуе(3)) * А 9% еуе(3) Ж = BT 
单位 矩阵 








B = 
3. 0000 – 8. 0000 – 6. 0000 
2. 0000 – 9. 0000 – 6. 0000 
– 2. 0000 12. 0000 9. 0000 


第 6 章 实 训 题 参考 答案 
1. 求 方程 组 的 通 解 
和 二 22 二 
2x, +X, —2x, — 2%, = 0 
XI =% —4x, -3x4 = Ü 
解 : MATLAB 程序 如 下 : 
format rat 
а=[1,2,2,1;2,1,-2,-2;1, -1, -4, -3] 
b = и (а,'г’) % KA EAE 
syms k1 k2 


x=kl1*b(:,1)+k2*b(:,2) % 写 出 方程 组 的 通 解 
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求 得 基础 解 系 为 
(2, -2,1,0)",(5⁄3, -4/3,0,1) 
ВА k (2, -2,1,0)" +k, (5/3, -4/3,0,1)",Ё,,Ё, є 


2. 求 超 定 方程 组 
2x, +4x, = 11 
3x, —5x, =3 
х +25, = 6 
2x, +х =7 
解 : MATLAB 程序 如 下 : 
a=[2,4;3,-5;1,2;2,1]; 
b=[11;3;6;7]; 
solution =а\Ь 
求 得 最 小 二 乘 解 为 x = (3. 0403 ,1. 2418)", 
3. 用 最 小 二 乘 解 法 解 方程 组 


xi +x, =l 








Xi +x; =2 


XI +x, +x, = Ü 


м +2x, —x, = -1 
解 : MATLAB 程序 如 下 : 
format rat 
a=[1,1,0;1,0,1;1,1,1;1,2,-1]; 
b=[1;2;0; -1]; 





xl =a\b % 这 里 \ 和 pinv 是 等 价 的 
x2 =pinv(a) * b 
求 得 最 小 二 乘 解 为 
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17 13 4 
маз Tp” = Е 
4. 求解 方程 组 
XI — X, — X, + X, = Ü 
X, — X, + X, — 3x, = ] 
Xi — x, —2x, +x, = — 1⁄2 


ЙЕ: MATLAB 程序 如 下 : 


format rat 
a=[1,-1,-1,1,0;1,-1,1,-3,1;1,-1,-2,3, 
2]; 


0 00 00 
故 方程 组 有 解 ,并 有 





1 
Xi =X + XÀ г 
1 
ia 2 Е 
因而 方程 组 的 通 解 为 
м 1 1 1/2 
X2 1 0 0 
= К, +k, + 
х3 0 2 1/2 
х4 0 1 0 
5. 求解 方程 组 : 





-1/ 
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25 + = SX, +^. =8 
xı —3x, -6x4 = 9 


2x, —x, +2^. = -5 





xı +4x, — Jx, + 6x, = 0 


解 : MATLAB 程序 如 下 : 
clc ,clear 
format rat 
а=[2,1, -5,1;1, -3,0, -6;0,2, -1,2;1,4, -7,6]; 
b=[8;9; -5;0]; 
В а=гапК(а),В В =тапк([а,Ь]) 
п = size(a,2) 
fR a==R B&R а==п 
fprintf HEHA НЕ) 
x=a\b 
elef R a==R B&R a<n 
fprintf( 方程 组 有 无 穷 多 解 \a ') 
x=a\b 狗 求 非 齐 次 方程 组 的 特 解 
хі =null(a,'r') % 求 齐 次 方程 的 基础 解 系 
else 
fprintf (方程 组 无 解 \ 方程 组 的 最 小 二 乘 解 为 \n ') 
х = а\Ь 
end 
求 得 唯一 解 为 x=(3, -4, -1,1) 。 
6. 求 方程 组 的 解 : 
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ху +х, – 3х; = = 1 
3x; – х, – 3х; +4x, =4 
ху +5х, 9х; — 8х, = 0 
解 : MATLAB 程序 如 下 : 
clc ,clear 
format rat 
а=[1,1, -3, -1;3, -1, -3,4;1,5,-9,-8]; 
Ь = [1;4;0]; 
В а=тапк(а),К В =гапк([а,Ь]) 
п = ѕіхе(а,2) 
КК a==R B&R a==n 
fprintt( 方程 组 有 唯一 解 \n) 
x =а\Ь 
elseif R_a==R B&R a<n 
fprintt( 方程 组 有 无 穷 多 解 \n ') 
x=a\b % 求 非 齐 次 方程 组 的 特 解 
хі =null(a,'r') 狗 求 齐 次 方程 的 基础 解 系 





else 
fprintf 方程 组 无 解 \a 方程 组 的 最 小 二 乘 解 为 \n) 
х =а\Ь 





end 
求 得 方程 组 的 通 解 为 
| 3⁄2 – 3⁄4 0 
X 3⁄2 7/4 0 
=k, +k, + ‚Е, Е R 
х; 1 0 -8/15 


х4 0 1 3/5 
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7. 设 有 线性 方程 组 
(1 +À)x, +x, +2, =Ü 
ж +(1+À)<, +% =3 
Xi +x, +(] +À)x; = À 
H A 取 何 值 时 ,此 方程 组 有 唯一 解 ? 
解 : MATLAB 程序 如 下 : 
clc ,clear 
syms k 
a=[1+k,1,1;1,1+k,1;1,1,1 +k]; 
D =det(a); % 求 系数 矩阵 的 行列 式 
y=factor(D) % 系数 矩阵 的 行列 式 进 行 因 式 分 解 
root = solve( y ) 
b = [0;3;К]; 
АВ = [а,Ь]; 
т = length ( root ) ; 








for n = 1:т 
fprintf('k =) ,disp(root( n) ) 
fprintft( 行 最 简 形 为 \n ) 
sln) =rref(subs( AB ,root(n) ) ) ; 
disp(sin) ) % 显 示 细 胞 数组 的 元 素 
if rank(subs(a,root(n))) ~ =rank( subs( АВ ,root(n))) 
fprintf(' 方程 组 无 解 \n '); 





else 
fprintf (方程 组 有 无 穷 多 解 \m ') ; 
end 
fprintf( ' акаев окоо ок око олова \n ') 


епа 
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ЖА Ч A= -3 H A0 时 ,方程 组 有 唯一 解 ; 当 入 = -3 
时 ,方程 组 有 无 穷 多 解 ; 当 A =0 时 ,方程 组 无 解 。 





第 7 章 实 训 题 参考 答案 
Т. 求解 线性 规划 问题 : 
min z= -2х 一 % +2, 
ху +х, +2х, =6 
х, +4x, — x, S4 
5. L. 
2) -2х, +х. <12 
х 20, x, 20, x, =5 
解 : MATLAB 程序 如 下 : 
> c=[ -2,-1,l];a=[1,4,-1;2,-2,1]; Ь=[4; 
12]; а =[1,1,2]; bl =6; 
>> №=[0;0; -inf]; ub= [| inf; mf; 5]; 
>> [x, z] =linprog(c,a,b,al,bl,1b,ub) 
运行 后 得 到 : 
x= 
4. 6667 
0. 0000 
0. 6667 
z= 
– 8. 6667 
2. 求 下 面 的 规 化 问题 : 


min -5x, -4x, -6x; 





sub. to Xi — X, +x, <20 


Зх, +2x, +4x, <42 
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3x, +2x, <30 
О=лх, ‚0 =х, ‚0 =< x, 
解 : MATLAB 程序 如 下 : 
>>f =| -5; -4; -6]; 
>>A=[1 -11;324;320]; 
> b=|[20; 42; 30]; 
>> 1Ь =zeros(3,1); 
>> [х,а] =linprog(f,A,b,[],L],b) 
结果 为 : 
х= Ф 
0. 0000 
15. 0000 
3. 0000 
fval = % 最 优 值 
– 78. 0000 


第 8 章 实 训 题 参 考 答案 

І. 计算 正 态 分 布 M0,1) 的 随机 变量 X 在 点 0. 6578 的 
密度 函数 值 。 

解 : MATLAB 程序 如 下 : 

>> pdf('norm ',0. 6578 ,0 ,1 ) 
ans = 

0. 3213 

2. 自由 度 为 8 的 卡 方 分 布 ,在 点 2. 18 处 的 密度 函数 值 。 

解 : MATLAB 程序 如 下 : 

>> pdf('chi2 ',2. 18,8) 


ans = 
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0. 0363 

3. 求 标准 正 态 分布 随 机 变量 X 落 在 区 间 ( — % ,0.4) 内 
的 概率 ( 该 值 就 是 概率 统计 教材 中 的 附 表 : 标 准 正 态 数值 
表 ) 。 








ff: MATLAB 程序 如 下 : 
>> cdf('norm ',0. 4,0,1) 
0. 6554 
4. 求 自由 度 为 16 的 卡 方 分 布 随 机 变量 落 在 [0,6.91 ] 内 
的 概率 。 
ff: MATLAB 程序 如 下 : 
>> cdf(' chi2',6. 91,16) 
0. 0250 
5. 设 X~N(3,22) 求 Plj2<X<51,P1-4<X<101,P 
o ` |X>3] 
解 : pl =Pi2 <Х<5} 
р2 =Р! -4<X<10| 
p3 =Р} [Х| >2} =1-Р{ |Х| <2} 
p4=P|X>3} =1-PIX<3| 


MATLAB 程序 如 下 : 
>> pl =normcdf(5,3,2) -normcdf(2,3,2) 
pl = 
0. 5328 
>> p2 =normcdf(10 ,3 ,2 ) -normcdf( – 4,3,2) 
p2 = 


0. 9995 
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>>p3 =1 -normcdf(2,3,2) -normcdf( -2,3,2) 
p3 = 
0. 6853 
>>p4 =1 -normcdf(3,3 ,2) 
p4 = 
0. 5000 
6. 设 随 机 变量 X 的 分 布 律 为 : 








求 E(X)E(X -1)。 

解 : 在 MATLAB 编辑 器 中 建立 M 文件 如 下 : 
X=[-2 -1012]; 
p=[0.30.10.20.10.3]; 

EX =sum(X. * p) 

Y=X. 2-1 

EY =sum(Y. * p) 
运行 后 结果 如 下 : 

EX = 

















EY = 
1. 6000 
7. 设 随机 变量 X 的 分 布 律 为 


X — 2 0 








P 0.4 0.3 
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ЖЕ(Х),Е(ЗХ? +5), 

解 : MATLAB 程序 如 下 : 

调用 累 和 函数 sum 

>>x=| -202]; 

>> pk = [0.4 0.3 0.3];sum(x. * pk) 

ans = – 0. 2000 

>>у=3 жх. ^2 +5 

y =17 5 17 

>> sum(y. * рк) 

апѕ = 13. 4000 

8. 设 随 机 变量 X 的 分 布 密度 函数 为 :f (x) = 
S (0<x<1) 
S. 5 

0, 其 他 

解 : MATLAB 程序 如 下 : 

>> syms x 

>> fx =3/5 +6/5 * x2; 

>> Ех =int(x = fx,0,1) % 采 用 连续 随机 变量 的 数学 期 

ZH PRAXI int(x *fx,a,b) 





, 求 EX 和 DX 


Ех =3/5 % 数学 期 望 
>> Ех2 =int(x2 * fx,0,1) 
Ex2 = 11/25 


>> Dx = Ех2 – Ех^2 

Dx =2/25 % 方 差 

9. 规定 某 型 电子 元 件 的 使 用 寿命 超过 15006 为 一 等 品 ， 
已 知 一 批 样品 20 只 ,一 等 品 率 为 0.2。 中 一 等 品 
元 件 的 期 望 与 方差 为 多 少 ( 即 一 等 品 元 件 的 个 数 的 最 大 值 )? 





. 146 . MATLAB 基础 与 应 用 

















解 : 由 分 析 可 知 电子 元 件 中 一 品 元 件 分 布 为 二 项 分 布 调 
用 binostat(n,p) 函数 

>> [m,v| =binostat(20 ,0. 2) 

m=4 

v =3. 2000 

10. 求 参 数 和 A =6 的 泊 松 分 布 的 期 望 与 方差 。 

解 : MATLAB 程序 如 下 

>> [m,v ] =poisstat(6) 





m =6 

v =6 

可 见 泊 松 分 布 的 期 望 与 方差 与 其 参数 入 相等 。 

11. 随机 抽取 6 个 深 珠 测 得 直径 (单位 :mm) 如下: 
14.70 15.21 14.90 14.91 15.32 15.32 
试 求 样本 平均 值 

解 : MATLAB 程序 如 下 : 

>>Х = [ 14. 70 15.21 14.90 14.91 15. 32 15.32]; 
>> теап(Х) % 计 算 样本 均值 

则 结果 如 下 : 

15. 0600 

12. 求 下 列 样本 的 样本 方差 和 样本 标准 差 .方差 和 标准 




















14.70 15.21 14.90 15.32 15.32 
解 : MATLAB 程序 如 下 
> X = [14.7 15.21 14.9 14.91 15. 32 15.32]; 
>> DX = уаг(Х,1) 2 方差 
DX = 
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0. 0559 
>> sigma = std(X,1) % 标准 差 
sigma = 
0. 2364 
>> DXI = var( X) % 样 本 方差 
DXI = 
0. 0671 
>> sigmal = std( X) % 样 本 标准 差 
sigmal = 
0. 2590 
13. 随机 的 取 8 只 活塞 环 , 测 得 它们 的 直径 (单位 :mm) 
为 : 
74.001 74.005 74.003 74.001 74.000 73.998 74.006 
74. 002 
求 样本 的 均值 方差 值 、 样 本 方差 值 标准 差 值 样本 标准 
差 值 。 
解 : MATLAB 程序 如 下 : 
>> d = [74.001 74.005 74.003 74.001 74.000 73.998 
74. 006 74. 002 ] ; 
>> mean( 1) 
ans = 74. 0020 
>> d = [ 74. 0010 74.0050 74.0030 74. 0010 74. 0000 
73. 9980 74. 0060 74. 0020 | ; 
>> уаг( 4,1) 
ans = 6. 0000е — 006 
>> уаг( 4) 
ans = 6. 8571е – 006 
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>>std(d,1) 
ans = 0. 0024 
>>std(d) 

ans = 0. 0026 


第 9 章 实 训 题 参考 答案 
Т. 假设 某 种 清漆 的 9 个 样品 ,其 干燥 时 间 ( 以 小 时 计 ) 分 
别 为 6.0,5.7,5.8,6.5,7.0,6.3,5.6,6.1,5.0。 设 干燥 时 间 
总 体 服 从 正 态 分 布 WUw,o ) 。 
(1) ио 的 点 估计 值 ; 
(2) Жи То 的 置信 和 度 为 0.95 的 置信 区 间 (o 未 知 ) 。 
解 : MATLAB 程序 如 下 
>>time =[6.0 5.7 5.86.57.06.35.66.15.0]; 
>> [| muhat, sigmahat, muci, sigmaci | = normfit ( time, 
0. 05) 
muhat = 
6 %p 的 估计 值 为 6 
sigmahat = 
0. 5745 %oc 的 估计 值 为 0. 5745 
muci = 
5. 5584 
6. 4416 %u 的 估计 区 间 为 [5. 5584, 6. 4416 | 
sigmaci = 
0. 3880 
1. 1005 Фо 的 估计 区 间 为 [5. 5584, 6. 4416] 
2. 分 别 使 用 金 球 和 铀 球 测定 引力 常数 : 
(1) 用 金 球 测定 观察 值 为 :6. 683 ,6. 681,6. 676 ,6. 678, 
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6. 679 ,0. 672; 


(2) 用 铂 球 测定 观察 值 为 :6. 661 ,6. 661 ,6. 667 ,6. 667, 
6. 664, 

设 测定 值 总 体 为 Wuw,o ) м Мо 为 未 知 ,对 (1) (2) Я 
种 情况 分 别 求 和 o 的 置信 度 为 0.9 的 置信 区 间 。 

解 : MATLAB 程序 如 下 : 

>> j = [6. 683 6. 681 6. 676 6. 678 6. 679 6. 672]; 

>>Ь = [6. 661 6. 661 6. 667 6. 667 6. 664]; 

>> [ muhat ,sigmahat ,muci,sigmaci] = normfit( j ‚0. 1) 





muhat = 
6. 6782 
sigmahat = 
0. 0039 
muci = 
6. 6750 
6. 6813 
sigmaci = 
0. 0026 
0. 0081 
% 金 球 测定 数据 的 置信 和 度 为 0.9 的 人 和 ez 的 置信 区 间 
X : и [ 6. 6750 ,6. 6813 |, с 0. 0026, 0. 0081 ] 
>> [muhat ,sigmahat , тисі ,sigmaci | = normfit( b ‚0. 1) 





muhat = 
6. 6640 
sigmahat = 
0. 0030 


muci = 
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6. 6611 
6. 6669 
sigmaci = 
0.0019 
0.0071 
% 铀 球 测定 数据 的 置信 度 为 0.9 的 人 和 的 置信 区 间 
为 :&[6. 6611 ,6. 6669 | ,с[0. 0019 ,0. 0071 ] 
3. 有 两 组 (每 组 100 个 元 素 ) 正 态 随机 数据 ,其 均值 为 
10, 均 方差 为 2, 求 95% 的 置信 区 间 和 参数 估计 值 。 
解 : MATLAB 程序 如 下 : 
>>r=normrnd(10,2,100,2); % 产生 两 列 正 态 随 机 
数据 
>> | mu ,sigma ,muci ,sigmaci | = normfit( г) 
则 结果 为 
ши = 
10. 1455 10. 0527 ”% 各 列 的 均值 的 估计 值 
sigma = 
1.9072 2.1256 ”% 各 列 的 均 方 差 的 估计 值 
muci = 
9.7652 9.6288 
10. 5258 10. 4766 
sigmaci = 
1.6745 1.8663 
2.2155 2.4693 
说 明 :muci ,sigmaci 中 各 列 分 别 为 原 随 机 数据 各 列 估计 
值 的 置信 区 间 ,置信 度 为 95% 。 
4. МЕР ЖЕЛЕ TEO F , 按 规格 平均 净重 379g ,标准 
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差 11g, 现 抽查 10 个 测 得 数据 (单位 : g):370.74 372. 80 
386. 43 398. 14 369. 21 381. 67 367. 90 371. 93 386. 22 393. 08 。 
试 说 明 平均 净重 是 否 符合 规定 要 求 (a =0. 05 ) 。 

解 : 总 体 的 均值 和 标准 差 已 知 ， 用 u 检验 法 ,由 函数 
ztest(x,mu,sigma,alpha) 实现 。 

假设 НО: р=379;Н1 :рз379 

MATLAB 程序 如 下 : 

>> p = [ 370. 74 372. 80 386. 43 398. 14 369.21 381.67 
367. 90 371. 93 386. 22 393. 08 |; 

>> | h,pweight,ci] =ztest(p,379 ,11) 











h=0 

pweieht = 0. 8154 

ci = 386. 6297 

在 a=0.05 水 平 下 ,不 能 拒绝 原 假设 , 即 平均 净重 符合 
规定 要 求 。 


5. 某 制药 厂 生产 一 种 抗生素 , 已 知 在 正常 情况 下 , 每 瓶 
抗生素 的 某 项 指标 服从 均值 为 22.3 的 正 态 分 布 。 某 天 开工 
后 , 测 得 10 瓶 的 数据 为 : 22.3 21.5 21.7 23.4 21.8 21.4 
23.4 19.8 24. 4 21. 2 ,该 天 生产 抗生素 每 瓶 的 该 项 指标 的 均 
值 是 否 正 常 ? 

解 : MATLAB 程序 如 下 : 

>>z = [22.3 21.5 21.7 23.4 21.8 21.4 23.4 19.8 24.4 
21.2]; 














>> std( z) 
апѕ = 1.3295 ”% 用 样本 标准 差 代替 总 体 标 准 差 
标准 差 方差 已 知 时 , 单个 正 态 总 体 的 样本 均值 是 否 等 于 


22. 3 的 检验 ,用 函数 ztest 实现 
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> [h,sig,ci] =ztest(z,22.3,(1.33)22) 

>> [h,sig,ci] =ztest(z,22.3,(1.33)22) 

h=0 

sig = 0. 7073 

cl = 

20.9936 23. 1864 % ДЕ а = 0. 05 水 平 下 ,不 能 拒 

绝 原 假设 , НП m =22. 3 时 该 天 生产 抗生素 每 瓶 的 该 项 指标 的 
均值 正常 。 

6. 某 加 油 站 有 2008 年 1 月 和 2 月 二 组 汽油 价格 数据 样 
本 :119 117 115 116 112 121 115 122 116 118 109 112 119 112 
117 113 114 109 109 118 和 118 115 115 122 118 121 120 122 
120 113 120 123 121 109 117 117 120 116 118 125 , 设 该 州 汽 
油价 格 的 标准 差 为 每 加 仑 4 美 分 。 

(1) 判断 1 月份 每 加 仓 汽油 平均 价格 是 否 为 1.15 美元 ; 
(2) 设 价格 标准 差 未 知 ,比较 1.2 月 的 二 个 样本 均值 。 

解 :(1) 方差 已 知 时 , 单个 正 态 总 体 均 值 的 检验 用 u № 
ре, АРА ztest(x,mu ,sig,alpha) 

假设 HO: u =m =1. 15 

MATLAB 程序 如 下 : 

>> р! = [119 117 115 116 112 121 115 122 116 118 109 
112 119 112 117 113 114 109 109 118]; 

>> р2 = [118 115 115 122 118 121 120 122 120 113 120 
123 121 109 117 117 120 116 118 125]; 

>> [В ,руаше,с1] =ztest(p1Z100 ,1. 15 ,0. 04) 

В =0 

pvalue =0. 8668 

сі = 1. 1340 1. 1690 % а = 0.05 水 平 下 ,不 能 拒 
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绝 原 假设 , 即 1 月 份 每 加 仑 汽油 平均 价格 为 1. 15 美元 。 

>> [В ,руаше,с1] =ztest(p2/100,1. 15 ,0. 04) 

в =1 

pvalue =9. 1112е — 005 

ci=1.1675 1.2025 Фа =0.05 水 平 下 ,拒绝 原 
假设 , 即 2 月 份 每 加 仑 汽油 平均 价格 与 1. 15 美元 有 显著 差 
异 。 

(2) 具有 相同 方差 的 2 个 正 态 总 体 均 值 差 用 t 检验 法 
检 , 用 函数 ttest2 (x,y,alpha) 实 现 

假设 HO:p1 =p2,H1 :р1 7р2, 

MATLAB 程序 如 下 

>> [h,sig,ci] = ttest2 ( р1/100 ,p2Z100 ) 

h=1 

sig = 0. 0083 

сі = - 0. 0578 -0.0092 % ДЕ a =0. 05 水 平 下 , 拒 
绝 原 假设 , 即 2 月 份 每 加 仑 汽油 平均 价格 与 1.15 美元 有 显 
Е о 

7. 从 一 批 化 工 产品 中 取 8 份 样品 ,每 份 一 分 为 二 ,分 别 
用 甲 、 乙 两 种 方法 检测 某 成 分 得 数据 (单位 :mg) :57 56 61 60 
47 49 63 61 和 65 69 54 60 52 62 57 60。 设 含量 测 值 NV (u, 
о?) , 问 两 法 测定 的 平均 值 是 否 有 显著 差异 ? (а =0. 05) 

解法 一 : 

构造 随机 变量 z=pl -Pp2， 即 两 法 测定 值 的 差 z, 仍 服 从 
正 态 分 布 。 

MATLAB 程序 如 下 

>>pl = [57 56 61 60 47 49 63 61]; 

>> р2 = [65 69 54 60 52 62 57 60]; 
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>>z=pl - р2 % z 两 法 测定 值 的 差 

z=-8 -13 7 0 -5 -13 6 1 

>> mean (z) 

ans = -3. 1250 % 两 法 测定 值 的 差 的 均值 
>> std( z) 

ans =7. 8819 % 两 法 测定 值 的 差 的 标准 差 





两 法 测定 值 的 差 z 的 均值 和 标准 差 已 知 , 用 u 检验 法 ， 
Н РА ztest(x,mu sigma „alpha ) 实现 。 

假设 H:z=0 %pl = p2 

>>z=[ -8 -13 7 0 -5 -13 6 1]; 

>> [h,sig,ci] = ztest(z, – 3. 1250, (7. 8819)^2) % $R 








认 值 =0. 05 
h=0 
sig = 1 


сі = – 46. 1742 39. 9242 % ДЕ а = 0.05 水 平 下 ,不 
能 拒绝 原 假设 , В pl = р2 

两 法 测定 的 平均 值 无 显著 差异 。 

解法 二 : 

用 t 检 验 法 检验 样本 均值 与 一 常数 比较 , z=pl -p2 与 0 
比较 

Б НО: = ш1 -u2 =0,Н1: = =ul -u2 #0 

>> [h,sig,ci] = ttest(z) 

h =0 

sig =0. 2991 

ci = -9.7145 3. 4645 % ДЕ а =0. 05 水 平 下 ,不 
能 拒绝 原 假设 , В рі = р2 

两 法 测定 的 平均 值 无 显著 差异 。 
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解法 三 : 

用 t 检验 法 检验 具有 相同 方差 的 2 个 正 态 总 体 均 值 差 ， 
用 函数 ttes2(x,y,alpha) 实 现 

假设 НО :ul -u2 =0,H1 :ul -u2 =0 

用 t 检验 法 检验 具有 相同 方差 的 2 个 正 态 总 体 均 值 差 ， 
用 函数 ttes2(x,y,alpha) 实 现 

>> [h,sig,ci] =ttest2(pl,p2) 

h=0 

sig = 0. 2942 

ci = – 9. 2746 3. 0246 % ДЕ a=0.05 水 平 下 ,不 
能 拒绝 原 假设 , B pl = p2 

两 法 测定 的 平均 值 无 显著 差异 。 

8. 近期 在 某 城市 一 中 学 收集 到 24 名 17 岁 男 性 中 学 生 
身高 数据 如 下 : 

170.1 179.0 171.5 173.1 174.1 177.2 170.3 176.2 

163.7 175.4 163.3 179. 0 176. 5 178.4 165.1 179. 4 

176.3 179.0 173.9 173.7 173.2 172.3 169.3 172. 8 

又 查 到 20 年 前 同一 所 学 校 同龄 男生 的 平均 身高 为 
168cm ,根据 上 面 的 数据 回答 ,20 年 来 城市 17 岁 男 性 中 学 生 
的 身高 是 否 发 生 了 变化 ? 

解法 一 :根据 题 意 : 原 假设 :有 :p=168 ЧБ: А, и 
#168 

MATLAB 程序 如 下 : 

>> x = [170.1 179.0 171.5 173.1 174.1 177.2 170.3 
176.2 163.7 175.4 163.3 179.0 176.5 178.4 165.1 179.4 
176.3 179.0 173.9 173.7 173.2 172.3 169.3 172.8]; 

>> [H,sig] =ttest(x, 168, 0. 05, 0) 
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H= 


Slg = 
8. 7275е – 006 
% 返 回 值 日 =1, 即 在 显著 水 平 为 0.05 的 情况 下 ,拒绝 原 
假设 。 即 认为 学 生 刁 高 的 确 发 生 了 变化 。 
解法 二 : 
>> x = [170.1 179.0 171.5 173.1 174.1 177.2 170.3 
176.2 163.7 175.4 163.3 179.0 176.5 178.4 165.1 179. 4 
176.3 179.0 173.9 173.7 173.2 172.3 169.3 172.8]; 
>> [H,sig,ci] =ttest(x, 168, 0. 05, 0) 
Н = 
1 
sig = 
8. 7275е – 006 
cl = 
171. 4660 175. 4340 % 置信 区 间 
% 返 回 值 日 =1, 即 在 显著 水 平 为 0.05 的 情况 下 ,拒绝 原 
假设 。 即 认为 学 生 身 高 的 确 发 生 了 变化 , 且 平 均 身 高 有 95% 
的 可 能 性 在 [171. 4660 ,175. 4340 ] 。 
9. 为 了 人 研究 某 一 化 学 反应 过 程 中 ,温度 XX 对 产品 生成 率 
了 的 影响 , 测 得 数据 如 下 : 
温 JÆ X 100 110 120 130 140 150 160 170 180 190 
生成 率 了 45 51 54 61 66 70 74 78 85 89 
ЕВЕ y =a + bx 型 回归 。 
解 : MATLAB 程序 如 下 : 
>>x=[100 110 120 130 140 150 160 170 180 190]; % 
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温度 
>>у = [45 51 54 61 66 70 74 78 85 89]; Фр 
>> plot(x,y,'*'); % 作 数据 散 点 图 
图 形 如 图 附 -5 所 示 。 
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图 附 -5 


>>p=polyfit(x,y,1) — % (解法 一 ) 计 算 直 线 拟 合 系 
数 

p =0. 4830 -2.7394 % 结 果 表 明 回 归 直 线 为 y = 
0. 4830х —2. 7394 


>> [р,з] =polyfit(x,y,1) %( 解 法 二 ) 计 算 直 线 拟 合 系 





数 
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p =0. 4830 -2.7394 % 结果 表明 回归 直线 为 y = 
0. 4830x —2. 7394 
== К; [ 2х2 double | 





4:8 % s 为 中 间 参 数 ,用 于 polyval 
函数 和 ployconf 函数 进行 预测 Y 的 误差 估计 。 
normr :2. 6878 


>>[Y,D] =polyval(p,x,s) % 给 出 预测 值 Y+D, 设 
ployfit 函数 的 输入 数据 Y 的 误差 服从 方差 为 常 值 的 正 态 分 
fi, 则 ployval 产生 的 Y+D 至 少 包 含 预 测 值 的 50% 。 

Y =45. 5636 50. 3939 55. 2242 60. 0545 64. 8848 

69. 7152 74. 5455 79. 3758 84. 2061 89. 0364 

D =1. 1023 1. 0618 1. 0304 1. 0089 0. 9980 

0. 9980 1. 0089 1. 0304 1. 0618 1. 1023 

>> [У, Оейа | = polyconf(p,x,s,0.05) % 给 出 100(1 
— alpha) % 置信 区 间 У = Рена, DZ ployfit 函数 的 输入 数据 Y 
的 误差 服从 方差 为 常 值 的 正 态 分 布 。 

Y =45. 5636 50. 3939 55. 2242 60. 0545 64. 8848 

69. 7152 74. 5455 79. 3758 84. 2061 89. 0364 

Delta =2. 5418 2. 4485 2. 3761 2. 3266 2. 3015 

2. 3015 2. 3266 2. 3761 2. 4485 2. 5418 

>> lsline % 画 最 小 二 乘法 拟 合 直 线 

图 形 如 图 附 -6 所 示 。 

10. 用 ployfit 函数 拟 合 y =a; + ax + a, x2 型 回归 
































ЙЕ: MATLAB 程序 如 下 : 
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图 附 -6 


>>x=[0.5 1.01.52.02. 53.0]; 

>>у= [1.75 2. 45 3. 81 4. 80 7. 00 8. 60]; 

>> plot(x,y,'*') 

图 形 如 图 附 -7 所 示 。 

>>p=polyfit(x,y,2) % (解法 一 ) 计 算 曲 线 拟 
合 系 数 

p=0.5614 0.8287 1.1560 %% 表 明 回归 直线 为 y = 
0. 5614x +0. 8287x +1. 1560 

>> [р,з] = рой (х,у,2) % (解法 二 ) 计 算 曲 线 拟 
合 系数 
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图 附 -7 


р=0. 5614 0.8287 1.1560 % 表明 回归 直线 为 》 
0. 5614х? +0. 8287x +1. 1560 

з=В: [ 3х3 double ] 

df:3 

normr :0. 4220 

>>xl = [0. 5:0. 05:3.0]; 

> yl = р(1) + х1. ^2 + р(2) + х1 + р(3); % у 
0. 5614x? +0. 8287x +1. 1560 二 次 曲线 回归 

>> plot( x1 ,yl,'—r') % '-r' 表示 红色 

图 形 如 图 附 -8 所 示 。 
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图 附 -8 


第 10 章 实 训 题 参考 答案 
L 用 最 小 二 乘法 求 一 个 形 如 y =a + bx 的 经 验 公式 ,使 
它 与 下 表 所 示 数 据 拟 合 。 


x 19 25 31 38 44 





y 19.0 32:3 49 73.3 97. 8 


№: MATLAB 程序 如 下 : 

>>х= [19 25 31 38 44|; 

>>y = [19.0 32.3 49.0 73.3 97. 8 ]'; 
>>r= [опеѕ(5 ,1) „х. 2]; 

>> ab =т\у 

>> х0 = 19:0. 1:44; 
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>>y0 =ab(1) +ab(2) * х0. ^2; 
Plot(x,y,'0',x0,y0,'r') 

求 得 的 经 验 公式 为 y=0.9726 +0. 05x 
2. 某 乡镇 企业 1990—1996 年 的 生产 利润 如 表 : 


年 份 1990 1991 1992 1993 1994 1995 1996 








利润 /万 元 70 122 144 152 173 


试 预测 1997 年 和 1998 年 的 利润 。 

解 : 作 已 知 数据 的 散 点 图 ,有 

>> x0 = [ 1990 1991 1992 1993 1994 1995 1996 | ; 

>> y0 = [70 122 144 152 173 196 202]; 

>> plot( xÜ ,y0 ,' * ') 

发 现 该 乡镇 企业 的 年 生产 利润 几乎 直线 上 升 ,如 图 附 -9 
所 示 , 因 此 可 以 用 y=aix + oo 作为 拟 合 函数 来 预测 该 乡镇 企 
业 未 来 的 年 利润 。 编 写 程序 如 下 : 

>> x0 = | 1990 1991 1992 1993 1994 1995 1996 | ; 

>> y0 = [70 122 144 152 174 196 202]; 

220 т 


196 202 


























2001 


* 


1807 
1607 











1990 1997 1992 1903 1994 1995 1996 
图 附 -9 
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>> а = ployfit (x0 ,y0,1); 

>> Ү97 = роууа[( а 1997) ; 

>> Ү98 = ployval( a ,1998 ) ; 

得 到 1997 年 的 利润 233. 4286, 1998 年 的 利润 253. 9286 


Жп 章 实 训 题 参考 答案 

小 李 是 某 一 大 学 的 应 届 毕 业 生 ,参加 多 家 企业 的 招聘 会 
后 ,有 甲乙 、 丙 三 个 单位 愿意 录用 他 。 请 你 为 他 提供 建议 ,以 
选择 合适 的 单位 。 

解 : 根据 往届 同学 的 就 业经 验 , 可 以 考虑 以 下 几 个 因素 : 
研究 能 力 发 展 前 途 、 待 遇 、 同 事情 况 ` 地 理 位 置 . 单 位 名 气 。 
甲乙 两 三 个 单位 为 供 选 择 的 工作 地 。 

(1) 建立 目标 为 工作 满意 程度 的 层次 结构 模型 。 如 图 
附 -10 所 示 。 























目标 层 0 


工作 满意 度 
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(2) 构造 两 两 判断 比较 矩阵 
比较 6 个 因素 对 工作 满意 度 的 影响 大 小 ,可 以 构造 比较 
НА: 











3 
А e P: 
5 
人 
2 2 
A= 
ЕС Е 1 
54 4 3 7 
1 1 a 
4 4 3 3 6 
А а ея 
2 4 





类 似 于 和 矩阵 4 的 构造 方法 , 男 外 构造 方案 层 3 个 因素 对 
准则 层 6 个 因素 中 的 每 一 个 因素 的 判断 矩阵 4;(i=1,2,3,4， 
5,6) 











1 1 1 1 т: S 
аа 12325 3 4 
4 2 4 5 
A= 4 1 3|4=|，， а fa =| 2 1 1| 
二 
= 5 2 1 ти 
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1 
A,=| 3 
1 


5 


1 1 7 9 
— 5 1 1 7 
3 l 
1 1 7 КА, d 
1 714, = „Ав=| 7 
1 1 
l — — 1 l 
= j 7 7 — 1 1 
7 9 


(3) 单项 一 致 性 检验 (用 MATLAB 程序 检验 ) 

>>р=[1,1,1,5,4,2/3;1,1,2,4,4,4/5;1,1/2,1,4,3, 
1/2;1/5 ,1/4,1/4,1,3/4,1/7;1/4,1/4,1/3,4/3,1,146;3/2, 
4/5,2,7,6,1] 


>>[S,T 


] =eig(P),Pmax =T(1,1),CI= (Pmax -6)/(6 


-1),RI=1.24,CR =CI/RI СК =0. 06781 


>> ВІ = 
5;4,1,1/2;5, 
>> B3 = 
5;3,1,7;1/5, 
>> В5 = 


[1,1/4,1/2;4,1,3;2,1/3,1],В2 = [1,1/4,1/ 
2,1] 

[1,2/3 ,3/4;2/3 ,1,1;3/4,1,1],В4 = [1,1/3, 
1⁄7,1] 


[1,1,7;1,1,7;1/7,1/7,1],B6 =[1, 7,9;1⁄/7, 





1,1;1/9,1,1] 


>> [S1,T1 ] =е1=( ВТ) В] тах =Т1(1,1),СП = (ВІ пах 


-3)/(3-1) 


‚ВП =0.58,CR1 = CH /RII 


>> | S2 ,T2 | =eig(B2),B2max =Т2(1,1),С12 = ( B2max 


-3)/(3-1) 


, RI2 =0. 58 ,CR2 = CI2/RI2 


>> [53 ,T3 | = eig( B3 ) ,B3max = Т3(1,1), СІЗ = ( B3max 


-3)/(3-1) 


,RI3 =0. 58 ,CR3 = СІЗ/КІЗ 


>> [S4,T4] =eig(B4),B4max =Т4(1,1),С14 = ( B4max 


-3)/(3-1) 





‚КИ =0. 58 ,CR4 = CI4/RI4 





>> [85 ,T5 | =eig(BS),BSmax =T5(1,1) ,С15 = (В5шах 


-3)/(3-1) 


,RIS =0. 58 ,CR5 = CI5/RI5 
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>> [S6,T6] = eig( B6 ) ,B6max = Т6(1,1),С16 =(B6max 
-3)/(3 -1) ,RI6 =0. 58 ,CR6 = CI6/RI6 

>> ZB = [ СВТ, СВ2 ,CR3 ,CR4 ,CRS ,CR6 ] 

ZB =0. 0158 0. 0212 0. 0462 0. 0559 0 0. 0061 

所 有 值 都 小 于 0. 1 ,通过 一 致 性 检验 

(4) 总 体 一 致 性 检验 

> Wl =[S(;:,1)/sum(s(:,1))] 

>> ZC = [ CH ,CD СІЗ, C ,CIS ,CI6 ] ,ZR = [ RII , RI2, 
RI3 ,RI4 ,RIS , RI6 | 

>>a=ZC x W1,b=ZR *=W1,CR =а/Ь 

>> CR =0. 0342 

由 于 СВ 的 值 小 于 0. 1 ,所 以 总 体 通 过 一 致 性 检验 

(5) 排序 

>> Wl =|[S(:,1)/sum(s(:,1))| 

>> W2 =[|S1(:,1)/sum(sl(:,1)),S2(;,1)/sum (s2 
(:,1)),S3(:,1)/sum(s3(:,1)),S4(:,1)/sum(s4(:,1)), 
S5(:,1)/sum(s5(:,1)),S6(:,1)/sum(s6(:,1))] 

>> W= W2 x Wl 

W =0. 3843 

0. 3512 
0. 2645 
从 结果 比较 , 故 应 选 第 一 个 单位 。 
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1. 鸡 兔 同 笼 ,共有 36 1525 100 只 脚 , 问 鸡 兔 各 有 多 少 ? 
解 : 程序 如 下 : 

for i=1:36 
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for j =1:36 
ifi+j==36 &2 *1+4 *] == 100 
а= [1,}] 
end 
end 
end 
2. 公元 5 世纪 我 国 古代 数学 家 张 丘 建 在 4 算 经 》 一 书 中 
提出 了 ” 百 鸡 问题 ”: 鸡 翁 一 值钱 五 , 鸡 母 一 值钱 三 , 鸡 雏 三 值 
钱 一 。 百 钱 买 百 鸡 , 问 鸡 翁 В} Е ЛЬ? 
解 : 设 x: 鸡 俩 数 , 则 x 的 范围 :0 ~19 
у: ВЕРЯ, Ш у 的 范围 :0 ~33 
2: ASHER, IJ z 的 范围 :0 ~ 100 
则 : 
x+y+z=100 
5х +3y +z/3 =100 
这 是 一 个 不 定 方程 。 
for x =0.19 
for y =0:33 
for z = 0.100 
Ш(х+у +z==100)&(5 * x +3 жу + z/3 



































end 


3. 猴子 分 桃 , 五 只 猴子 ,第 一 只 将 桃子 平均 分 成 五 份 ,多 








. 168 . MATLAB 基础 与 应 用 








了 一 个 扔 掉 ,并 拿 走 一 份 ,第 二 只 猴子 ,把 剩 下 的 平均 分 成 五 
份 ,又 多 了 一 个 扔 掉 , 拿 走 一 份 。 第 三 ,第 四 ,第 五 只 猴子 都 这 











样 做 的 , 问 原 来 一 共有 多 少 个 桃子 ? 
s=1; 
а=[0,0,0,0,0]; 
while s >0 





c =s; 

for i=1:5 

a(i) =rem(c,5); 
с= (с-1) * (4/5); 
end 
На==[1,1,1,1,1] 
break; 

end 

s=s+l; 

end 


s 





4. 求解 猴子 吃 桃 问题 。 猴 子 在 第 一 天 摘 下 若干 个 桃子 ， 
当即 就 吃 了 一 半 , 又 感觉 不 过 瘾 ,于 是 就 多 吃 了 一 个 。 以 后 每 
天 如 此 ,到 第 10 天 再 想 吃 时 , 却 发 现 就 只 剩 下 了 一 个 桃子 。 














请 编程 计算 第 一 天 猴子 摘 的 桃子 个 数 。 
5=1; 
while s >0 
c =s; 
fori=1.:.9 
c=c/2-1; 


end 
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if c ==1 break 
end 
s=s+l; 


end 


s 
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